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STATEMENT  OF  MISSION 


The  mission  of  the  Office  of  Telecommunications  in  the  Department 
of  Commerce  is  to  assist  the  Department  in  fostering,  serving,  and 
promoting  the  nation's  economic  development  and  technological 
advancement  by  improving  man's  comprehension  of  telecommuni- 
cation science  and  by  assuring  effective  use  and  growth  of  the 
nation's  telecommunication  resources. 


In  carrying  out  this  mission,  the  Office 


• Conducts  research  needed  in  the  evaluation  and  development 
of  policy  as  required  by  the  Department  of  Commerce 

• Assists  other  government  agencies  in  the  use  of  telecommuni- 
cations 

• Conducts  research,  engineering,  and  analysis  in  the  general 
field  of  telecommunication  science  to  meet  government  needs 

• Acquires,  analyzes,  synthesizes,  and  disseminates  information 
for  the  efficient  use  of  the  nation's  telecommunication  re- 
sources. 

• Performs  analysis,  engineering,  and  related  administrative 
functions  responsive  to  the  needs  of  the  Director  of  the  Office 
of  Telecommunications  Policy,  Executive  Office  of  the  Presi- 
dent, in  the  performance  of  his  responsibilities  for  the  manage- 
ment of  the  radio  spectrum 

• Conducts  research  needed  in  the  evaluation  and  development 
of  telecommunication  policy  as  required  by  the  Office  of  Tele- 
communications Policy,  pursuant  to  Executive  Order  11556 
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COMPUTER  SOFTWARE  FOR  EWCS  PERFORMANCE  PREDICTION 


* 

E . J . Dutton 


A computer  program  in  FORTRAN  IV  language 
entitled  PROGRAM  PREDIC  has  been  prepared.  This 
calculates  performance  prediction  and  its  5 and 
95  percent  confidence  levels  for  microwave 
terrestrial  links  (8  to  30  GHz)  operating  in  the 
European  Wideband  Communication  System  (EWCS) . 
This  program  predicts  atmospheric  attenuation, 
principally  due  to  rain,  at  the  indicated  fre- 
quencies . 

Key  Words:  Microwave  Terrestrial  Links,  Europe, 
Path  Loss  Prediction,  Prediction 
Variability,  Rainfall. 


1.  INTRODUCTION 

Given  minimal  input  data  on  a certain  location  of  its  zone, 
latitude,  longitude,  and  the  elevation  of  the  station  to  within 
"standard"  map  accuracy  (assumed  as  100  ft.  or  35  m) , it  is 
intended  to  predict  link  performance  and  its  5 and  95  percent 
confidence  limits.  Minimal  link  information  of  carrier  fre- 
quency, f,  and  path  length,  L,  is  also  required.  The  prediction 
methodology  is  the  heart  of  the  conputer  software,  but  in  order 
to  achieve  prediction  via  so  few  input  parameters,  an  inter- 
polation procedure,  called  SUBROUTINE  IDBVIP  (Akima,  1975),  is 
also  required. 

This  subroutine  is  a highly  sophisticated  interpolation 
procedure,  without  which  historical  data  of  the  other  input 
parameters  of  average  annual  pressure,  temperature,  relative 
humidity,  thunderstorm  ratio,  8,  and  precipitation  M,  would  be 
required  for  each  location  of  interest.  Of  course,  as  can  be 
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seen  from  the  listing  of  the  Appendix,  an  input  option  is  avail- 
able to  the  user  if  any  of  these  parameters  are  available  at 
any  given  location.  therwise,  these  data  are  interpolated 
from  the  data  of  the  249  European  data-station  sample  (Dutton 
et  al. , 1974) . 


There  are  essentially  two  predictable  uncertainties 

associated  with  link  performance  prediction.  First  there  are 
2 2 

and  Sg,  the  predicted  spatial  and  temporal  variation  of  M 
and  within  a zone.  These  are  values  that  would  exist  even 
if  the  interpolation  procedure  were  perfect,  as  outlined  in 
Dutton  et  al.  (1974) . The  fact  that  the  interpolation  pro- 
cedure is  most  likely  not  perfect,  however,  introduces  the 

2 2 

second  uncertainties  and  Sg^.  The  interpolation  uncertainty 
is  of  consequence  for  the  M and  B parzuneters,  which  bear  directly 
on  rain  rate  prediction.  However,  it  is  not  of  much  consequence 
for  the  other  parameters,  whose  bearing  upon  the  prediction  pro- 
cess is  decidedly  minor.  Hence,  anytime  IDBVIP  is  used  to 
obtain  M and  6,  the  second  uncertainty  must  be  considered  in 
the  confidence  level  prediction  procedure.  Since  historical 
data  are  often  not  available  for  arbitrary  locations,  there  is 
usually  no  alternative  to  this  interpolation. 

We  have  estimated  this  prediction  uncertainty  by  zones, 
choosing  a "location"  to  be  one  of  the  stations  of  the  zonal 
data  sample,  and  using  other  nearby  stations  to  estimate  M and 
6 from  IDBVIP.  When  the  interpolated  values  of  M and  6 are 
compared  with  actual  station  values  for  all  stations  in  a zone, 

the  zonal  rms  uncertainty  results.  It  is  then  not  difficult 

2 2 2 2 
to  show  that  and  Sg^  can  be  combined  with  and  Sg,  the 

predicted  spatial-temporal  variances  of  M and  6,  as 


and 


2 2 
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M Me 
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(1) 

(2) 
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2 2 

to  produce  new  variances,  a.,  and  a„  of  M and  g to  be  used  in  pre- 

” ^2  2 

dieting  rain-rate  variance  instead  of  S„  and  S . Note  that,  in 

2 2 ” ° 2 2 

(1)  and  (2)  , and  are  assumed  independent  of  and  , 
Mg  Me  ge 

respectively.  The  estimated  values  of  Sj^,  S^,  and  are 

given  in  Table  1 for  those  zones  pertinent  to  the  EWCS  net. 

The  geographical  distribution  of  zones,  of  the  average 
annual  precipitation,  M,  and  of  the  thunderstorm  ratio,  g,  were 
obtained  previously  by  Dutton  et  al.  (1974)  . They  are  shown, 
respectively,  in  figures  1,  2,  and  3. 

There  are  certain  limitations  of  which  the  user  should 
be  aware  in  this  computer  software.  First,  it  is  not  intended 
for  use  outside  of  the  EWCS  network  grid  in  Europe,  or  approxi- 
mately an  area  enclosed  by  the  latitudes  52® O'  on  the  north  and 
42®0'  on  the  south,  and  by  the  longitudes  1®20'  on  the  west  and 
25® O'  on  the  east.  Second,  although  the  contribution  of  the 
gaseous  atmosphere  was  included  in  Methods  1 and  2,  it  was  not 
included  in  Methods  3 and  4.  However,  in  the  8 to  30  GHz 
region,  this  is  not  a severe  limitation  for  the  percentages  of 
an  average  year  (0.01  to  1 percent)  considered  by  the  software 
package . 

2.  FORTRAN  PROGRAM  PREDIC 

The  actual  software  package  for  the  path  loss  prediction 
and  its  variability  is  a FORTRAN  program  entitled  PREDIC. 

This  program  and  its  concommitant  subprograms  are  described 
in  detail  by  means  of  the  listing,  and  comment  cards  in  the 
Appendix.  A brief  flow  diagram  of  PROGRAM  PREDIC  is  also 
given  at  the  end  of  the  Appendix.  Only  a overview  of  the 
program  is  given  here. 

The  location  input  variables  discussed  in  section  1 are 

used  to  obtain  a point  distribution  of  rain  rates,  R,  and 
2 

its  variance  S_.  They  are  obtained  from  SUBROUTINE  DELTT. 

2 2 

This  subroutine  requires  the  zonal  and  Sg,  which  are  stored 
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in  arrays  in  the  main  program.  A SUBROUTINE  TABLES  is  used  as 
an  intermediary  to  achieve  the  interpolations  discussed  in 
section  1.  It  stores  data  from  each  of  the  249  European  data 
stations;  to  be  used  as  needed  by  SUBROUTINE  IDBVIP. 

There  are  four  procedures  for  estimating  mean  (50  percent 
confidence  level)  path  loss  due  to  rain  noted  in  comment  cards 
in  the  Appendix,  that  are  then  incorporated  into  PREDIC,  with 
the  aid  of  FUNCTIONS  GAMMA  and  REDCO,  and  SUBROUTINES  RAINRT, 
ATCOS,  and  PROMO.  Briefly,  FUNCTION  GAMMA  obtains  the  attenua- 
tion coefficient  per  unit  rain  rate  used  for  the  Barsis  et  al. 
(1973)  and  the  Battesti  et  al.  (1971)  methods.  SUBROUTINE  RAINRT 
performs  the  modification  of  rain  attenuation  prescribed  by 
Barsis  et  al  (1973)  for  terrestrial  links.  FUNCTION  REDCO  obtains 
the  necessary  parameter  for  modification  of  terrestrial  link 
rain  attenuation  prescribed  by  Battesti  et  al.  (1971). 

SUBROUTINES  ATCOS  and  PROMO  obtain  and  modify  terrestrial  link 
rain  attenuation  in  accordance  with  two  new  procedures. 

The  user,  through  an  integer  parameter  called  METHOD,  is 
able  to  select  whichever  of  the  procedures  he  wishes  to  run, 
from  whence  he  will  acquire  a mean  path  loss  due  to  rain  and 
atmospheric  gases  for  eleven,  critically  placed,  percentages 
(PCT)  of  an  average  year.  The  variances  and  confidence  levels 
accounting  for  year-to-year  (temporal)  variation  and  location- 
to- location  variation  within  the  specified  zone  (spatial) 
variation  are  then  computed  for  the  same  values  of  PCT,  using 
the  results  of  SUBROUTINE  DELTT.  The  resulting  5 and  95 
percent  confidence  levels  of  the  path  loss  prediction  for  the 
designated  values  of  PCT  are  then  printed,  along  with  the  pre- 
dicted mean  path  loss. 
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3 . RECOMMENDATIONS 

Some  attenuation  data,  taken  in  the  Federal  Republic  of 
Germany  (F.R.D.)  are  presented  in  figures  and  6 (CCIR,  1970- 
1974) . These  data  are  compared  with  four  possible  prediction 
methods,  designated  in  the  Appendix,  as:  Methods  1 and  2, 
corresponding  to  two  extrapolations  of  the  earth/space  pro- 
bability modification  factor  (Dutton  and  Dougherty,  1973)  to 
terrestrial  link  application;  Method  3,  the  method  of  Barsis  et 
al.  (1973);  and  Method  4,  the  method  of  Battesti  et  al.  (1971). 

Based  on  figures  4 and  6,  then,  it  is  the  personal 
recommendation  of  the  author  that  Method  4 , the  method  of 
Battesti  et  al.  (1971),  be  tentatively  given  preference  in 
path  performance  predictions  via  PROGRAM  PREDIC.  This  is 
because  it  fits  the  data  of  figures  4 and  6 reasonably  well. 

It  appears  from  figure  4 that  the  Method  3 the  method  Barsis 
et  al.  (1973)  is  best.  In  figures  5 and  6 it  appears  that 
Methods  1 and  2 are  best.  Based  on  these  results,  and  the  fact 
that  they  are  his  methods,  the  author  would  naturally  have  a 
predisposition  to  select  either  Method  1 or  2.  However,  the 
limited  size  of  the  data  sample  analyzed  is  still  justification 
for  caution  in  this  regard.  Furthermore,  Method  4 was  developed 
in  France,  and  thus,  presumably,  has  applicability  to  that  part 
of  Europe.  For  the  present,  the  author  recommends  Method  4;  as 
additional  experimental  data  become  available,  this  recommenda- 
tion should  be  reconsidered. 
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Figure  5.  Comparison  of  prediction  methods  with  observed  data  on  a 69.0  km  link  at  9.3  GHz  between 
Hambach  and  Darmstadt,  F.R.D.,  Dec.  1968  to  Aug.  1971  (except  Jan.  1969).  See  Figure  4 
for  explanation  of  insets  (a),  (b) , (c) , (d) , and  curves  (1),  (2),  (3),  (4). 
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PilOGRAM  PREDIC( INPUT,  OUTPUT, TAP?5  = INPUT, TAPE6  = 0UTPUTJ 


this  program  ESTIHATES  atmospheric  attenuation  and  its  5 ANO  95 
PERCENT  CONFIDENCE  LIMITS  ON  MICRONAVE  <8  TO  30  GHZ)  TERRESTRIAL 
LINKS  FOR  THE  EUROPEAN  NIOEBANO  COMMUNICATION  SYSTEM  (ENCS) , SY 
MEANS  OF  FOUR  DIFFERENT  PROCEDURES.  PRIMARY  EMPHASIS  IS  GIVEN  TO 
RAINFALL-CAUSED  ATTENUATION. 


INPUT  — 

TNO  CAPOS  READ. 


FIRST  CAPO  — 

STATIC  - ALPHANUMERIC  ARRAY  FOR  INPUT  OF  ANY  IDENTIFYING  HEADER 

DR  COMMENTS  (COLS  1-80,  8A10), 

SECOND  CARO  -- 

XLAT  - DEGREE-MINUTES  (DO. MM)  LATITUDE  OF  DESIRED  LOCATION  (COLS 
1-8,  F8.Q). 

XLON  - CEGREE-MINUTES  (DO. MM)  LONGITUDE  OF  OESIREO  LOCATION  (COLS 
9-16,  F8.0), 

ELEV  - ELEVATION  IN  METERS  AT  DESIRED  LOCATION,  TO  MAP  ACCURACY, 
ASSUMED  NEAREST  100  FT  OP  35  M (COLS  17-24,  F8.0). 

F - CARRIER  FREQUENCY  IN  GHZ  OF  TRANSMISSION  LINK  (COLS  25-32, 
FS.O). 

D - DISTANCE  ALONG  TRANSMISSION  PATH  (COLS  33-40,  F8.0). 

IZONE  - METEOROLOGICAL  ZONP  APPLYING  TO  DESIRED  LOCATION  (COLS  43- 
44,  12). 

PI  - average  annual  surface  pressure  IN  MILLIBARS,  IF  KNOWN 
(COLS  47-54,  E8.r». 

RHl  - AVERAGE  ANNUAL  SURFACE  RELATIVE  HUMIDITY  AS  A DECIMAL  FRAC- 
TION, IF  KNOWN  (COLS  f3-70,  F8,0). 

Ml  - AVFRAGE  ANNUAL  PRECIPITATION  IN  MILLIMETERS,  IF  KNOWN  (COLS 
71-79,  FJ.O). 

NOTE  — IF  ANY  OR  ALL  OF  THE  APOVE  METEOROLOGICAL  DATA  (PI, RHl, Ml) 
FOR  THE  CESIREO  LOCATION  IS  UNKNOWN,  LEAVE  THE  FIELD  FOR 
THE  UNKNOWN  PARAMETER  PLANK. 


C0MM0N/RRATE/RP(12) ,VRP(12) ,PCT(12) 

COMMON/FROM/TAUDBT  (1")  .RPVTAUd'’)  ,RELI(  15)  ,HT0P(12) 

DIMKNSICN  NUMSTA (10 ) , ST ATIO ( 8 ) , T AU5 ( 12 , 4 ) ,TAU95(12,4) ,TM00(12,4) , 
IVARBETdC)  ,VARM(10)  ,VTAU(12  (4)  , AT  ( 12 ) , RHS0ET  (1 C ) , RMSM  (1 0 ) 

PEAL  N,M1 

DATA  NUMSTA (1)  .NUMSTA (2) .NUMSTA (3)  .NUMSTA (4) .NUMSTA (5) .NUMSTA (6)  , 

1 NUMSTA (7) .NUMSTA (8) , NUMSTA (9 ) ,NUMS TA d 0 ) /3 , 54, 14 , 58 , 2 0, 10 , 1 ,69, 9 , 
213/ 

DATA  RMSeETd) .RMSBET (?) .PMSEET (3) .PMSBET ( 4 ) , RMS9ET (5 ) , RMSBET (6 ) , 
1RMS9ET(7) ,RMSBET(«) .RMSBET (9) .RMSBET ( 1 0 ) /O . , 0 . 06 , 0 . , . 04 , 0 . , . 06 , 0 . , 

2. 05. . 05. .05/ 

TATA  RMSM (1), RMSM (2) ,PMSM(3) ,RMSM(4) , RMSM ( 5 ) , RMSM ( 6 ) , PPSM ( 7 ) , RMSM ( 
19) , RMSM (9), RMSM( 10 )/0. 0.149. 3, 0.0, 297.4, 0.0,835. 6, 0.0,359. 3, 306.2, 
2372.2/ 

TATA  VARB'^Td)  , VARPET  (2)  ,VARPET(3)  , V ARBET  ( 4 ) , V ARBET  (5  ) , VARBET  ( 6 ) , 
IVARPET  (7)  .VARBET  (A  ) , VARBET  (9)  , VARPET  (1 0 )/.  0 99''2,  . 0 0 7762  ,.  0 09643 , 

2. 00 7106..  007744. .009092.0.0. .004706..007957..0 09467/ 

-<ATA  VAPMd)  ,VARM(2)  ,VAPM(3)  ,VARM(4)  , V ARM( 5 ) , V ARM  (6 ) , VARH (7  ) , 

IV ARM (8) ,VARM(9)  ,VARM(10) /7549.e,64669. , 493 87. , 5 21 32. , 186 76, , 396765 
2.  ,0.0, 87599.,  18165 5. ,196  97./ 


o o o o o r)  o o nnn  nnn  nan  nnn 


.....<{EAO  INPUT  DATA. 

READ (5*10001 ISTATIOITt 

REAO(5,1300)XLAT,XLON,ELEV.P,0,IZONE.P1*T1.RH1,N1 

. . .OBTAIN  INTERPOLATED  lETEOROLOGICAL  DATA  (KNOMN  DATA  AT  DESIRED 
. . .LOCATION  IS  GIVEN  PRIORITY*.  USER  IS  AGAIN  CAUTIONED  THAT.  IP 
. . .DATA  IS  UNKNOMN,  TO  LEAVE  THE  SPACES  POR  THE  DATA  ON  THE  INPUT 
. . .CARO  BLANK.  THIS  IS  THE  INPUT  OPTION  HENTIONEO  IN  THE  MAIN  TEXT. 

CALL  TAeLES(XLAT,XLON,ELEV,P,T,RH,BFT,H) 

IP  (Pl.GT.P.O)  P*P1 
IP(Tl.GT.O.O)  T=T1 
IP(RHl.GT.O.O)  RHsRHl 
IP(Ml.GT.O.O)  M=H1 

. . .OBTAIN  VARIANCES  OF  H AND  BETA. 

ENaNUHSTAdZONFJ 

VBFTA  = VARBFT(IZONE) ♦RHSBFT  aZONE»*PMSBET(IZONE» 

VEM=VARH(TZONEI  ♦RHSM ( IZONFJ •RMSH ( I ZONE ) 

GAM*(1.14-0.07*(  IF-Z.  ***(l./3.)  »)•  (!.♦!). 085* «F-3.5»*EXP(-.0  06*F*F» 
1) 

IF(GAH.LT.l.O)  GAH^l.O 
CAVsGAHHA(F) 

HAVa29.97B3/F 

CALL  OFLTT(BFT,H,VEM,VBETA* 

00  105  1=1,11 
105  RELT(I)=PCT(I) 

RELI (12)=100. 

IF  METHOD  = 1,  THE  PROBABILITY  MODIFICATION  FACTOR,  PTl,  IS  USED. 

.....IF  METHOD  = 2,  THE  PROBABILITY  MODIFICATION  FACTOR,  PT2,  IS  USED. 

IF  METHOD  = 3,  THE  METHOD  OF  0ARSIS  ET  AL.  (1973)  IS  USED. 

IF  HETHOO  = 4,  THE  METHOD  OF  BATTESTI  ET  AL.  (1971)  IS  USED. 

.....CALCULATIONS  OF  ATTENUATION  FOP  METHODS  3 AND  4. 

DO  100  1=1,11 
TF(RR(I).EQ.P.C)GOTO  220 

GOTO (165, 170, 175, ISO, 185, 190, 195, 200 ,2  05,210) , I ZONE 
165  HTOP(I) *0.344144*rr(I) ♦1.14796 
GOTO  215 

170  HTOP (I ) =9.46601758* (RR (I) ••.182178) 

GOTO  215 

175  HTOP(I)=10.011935574*(RR(I)** 0.171 134) 

GOTO  215 

180  MTOP( I) *10.44433786 3* ( PR (1 )••. 1419197 ) 

GOTO  215 

185  MTOP (I) *14.91071378 5* (RP(I) ••.0765 891) 

I GOTO  215 

) 190  MTOP(I)*10.63765431*(RR(I)**. 100104) 

GOTO  215 

j 195  HTOP(I)*14.0348C057*(RR(I)**C.C690994) 

■ GOTO  215 


200  HT0P(I)*5.686528747*(9R(I)**.213556) 
GOTO  215 


210  HTOP(I)=12.19«»316  71*  (RRII)**.100871) 

215  IF»HTOP«I».LT.9.0)  HTOPJI)»9.0 
GOTO  225 
220  HTOP<II«9.0 
225  HTOPI12I*10.8 
C 

C.... .CALCULATE  METHOD  3. 

C 

METH00*3 

CALL  RAlNRT(RRa)  .D.RM) 

OMsO 

rF(OM.GT.22.0»  DM=22.0 

TMOO <1, METHOD) «C AY* (RM**GAH)*DM 

VTAUa.METHOO)  = CTMOO  (I,  METHOD)  •TM00(I.METH00)*GAM*GAM/CRR«I)*RR(I) 
1) )*VRRII) 

C 

C CALCULATE  METHOD  <». 

C 

method*** 

RF»REOCO«n,PCT(I) ) 

TMOC(I,METHOO)=CAr*( (PF’RRII) )**GAM)*0 

VTAUd.METHOD)*  (TMOD(I,METHOP)*TMOOtI,METHOD)  ’GAM’D  AM/ ( RR  ( 1 ) •RR  ( I ) 
1) )*WPR)I) 

CALL  ATCOS(F,T,P,RH,HTOPIl) , PEL H I ) » BET  ,RR ( I) ,AT(I) ,MAV) 

100  TAUDBT(I) *AT(1) *0 
C 

CALCULATIONS  OF  ATTENUATION  FCR  METHODS  1 AND  2. 

C 

CALL  ATCOS«F,T,P,RH,10.0,RELI(12)  . BET , 1 . E-E  , AT ( 12)  .WAV) 

TAU0eT(12)=AT(12)*C 

DO  160  METHOD*!. 2 

CALL  PROMOCMAV.O, METHOD) 

00  160  1=1,11 

TM00( I, METHOD) =RE VT AU ( I ) ’TAUt  BT(12) 

160  VTAU(I,METHOO)=tTMOO(I,METHOO)*TMOO(I,METHOO)*GAM*GAM/(RR(I)*RR(I) 
1) )*WRRII) 

C 

C CALCULATE  «=  AND  95  PERCENT  CONFIDENCE  LIMITS  OF  ATTENUATION 

C DISTRIBUTION  FOR  ALL  METHODS  Cl-i.), 

C 

00  150  METHOn  = l,«, 

on  150  1=1,11 

X = 1.645*SQRT(  VTAUd  ,MFTHOn)  ) 

TAU95  « I,  METHOD)  *TMOnd,  METHOD  )-X 
IF  IT AU9E d, METHOD) .LT.O.O)  T AU95 ( I , MET  HOD) =0. 0 
I'O  T AUEd,  METHOD)  * TMOD  d,  METHOD)  ’X 

c 

C.  . .OUTPUT  HEADERS  AND  RESULTS. 

C 

WRITE  <6, 1400)  (ST ATI Od)  ,1  = 1,8) 

00  300  J=l,4 
IF(J,E0.1)  HPITE(6,1450) 

TF(J.EQ.2)  WRITE (6,1451) 

IF(J.E0.3)  WRITE (6,1452) 

IFIJ.EQ.4)  MRITF(6,1453( 

WRITE  (6, 1500)  (PCTd)  ,1*1,11) 

WRITE (6,2000) (PR (I) ,1*1,11) 

WRITE  (6,2500)  (VRR ( I ) , 1*1 , 11) 

WRITE (6,3000) I TMOD ( I, J) ,1*1,11) 

WRITE (6, 3500) ( VT AU ( I , J) , 1=1 , 1 1 ) 


• I 

I ! 


« i 


MVITE(6«l>000l  ITAU95 
300  NRITE(6«<»500)  (TAU5  ( I « J)  « I>1 1 11) 

STOP 

1000  FORMATIOAIOI 

1300  FORHATISF8.0,2XI2.2X4F9.0) 

1400  FORNATIlHltSAlO///) 

1450  FORHAT(51H  THE  PROBABILITY  HOOIFICATION  FACTOR,  PTl,  IS  USED.,//) 

1451  FORHATI51H  THE  PROBABILITY  HOOIFICATION  FACTOR,  PT2,  IS  USED.,//) 

1452  F0RHATI44H  THE  HETHOO  OF  BARSIS  ET  AL.  11973)  IS  USED.,//) 

1453  F0RNAT(46H  THE  HETHOO  OF  BATTESTI  ET  AL.  (1971)  IS  USED.,//) 

1500  F0RNAT(1X,12HPCT  3 ,11(F9.3,1X)/) 

2000  F0RHATI1X,12HR(HM/HR)  s ,11(F9.3,1X)/) 

2500  F0RHAT(1X,12HVAR(R)  > , 11 (F9. 3 , IX ) / ) 

3000  F0RHAT(1X,12HATTEN(09)  > ,11 (F9.3,1X) /) 

3500  F0RNAT(1X,12HVAR(ATT)  s ,11(F9.3,1X)/I 
4000  F0RNAT(1X,12HATT. (95)  : ,11 (F9. 3, IX) /) 

4500  F0RMAT(1X,12HATT. (5)  = , 11 ( F9. 3 , IX ) /// » 

END 


r>  o o 


SLOCK  DATA 

THIS  SUBROUTINE  INITIALIZES  ARRAY  PCT  IN  COMMON/RRATE/. 
C0MN0N/RRATE/RR(12) tVRR«12»»FCTI12) 

OAT A PCT(ll«PCTf2l tPCT(3l«PCT(A),PCT(5) ,PCTC6» ,PCT(7) ,PCT(6> ,PCT(9 
II  fPCTdO)  ,PCT(11)  ,PCT(12)/.01«.01S,.02f03f  .05t.08,.lf.2..St.S,l.  . 

e.aol/ 

END 
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SUBROUTINE  OELTI I BET A tEN« VN* VBET) 

THIS  SUBROUTINE  USES  THE  NETHOO  OF  DUTTON,  DOUGHERTY  AND  MARTIN, 
(197%},  PREDICTION  OF  EUROPEAN  RAINFALL  AND  LINK  PERFORNANCE 
COEFFICIENTS  AT  8 TO  30  GHZ,  SECTION  5,  PLUS  SOME  NEM  VARIANCE 
PREQICTION  PROCEDURES,  TO  OBTAIN  VARIATION  OF  T^l  NIN.  RAIN  RATES 
IN  TERMS  OF  ESTIHATEO  STANDARD  DEVIATIONS,  BASED  ON  CURRENTLY 
AVAILABLE  YEAR-TO-YEAR  PRECIPITATION  DATA. 

INPUT 

BETAbTHUNOERSTORH  ratio  at  EACH  STATION. 

E'1«AVERAGE  ANNUAL  PRECIPITATION  IN  MM,  AT  EACH  STATION, 

VM«VARIANCE  OF  YEAR-TO-YEAR  ANNUAL  PRECIPITATION.  IT  NAS 
NOT  PRfS«'NTEn  IN  THE  ABOVE  REPORT  AND  IS  THUS  REQUIRED  AS 
AOniTIONAL  INFORMATION. 

VBET»VARIANCE  OF  BETA 

OUTPUT  (IN  COMMON/RRATE/J 

RR  =T-MINUTE  RAIN  RATE  FROM  THE  MODIFIED  R-M  MODEL, 

VRR  sESTIMATEO  VARIANCE  OF  YEAR-TO-YEAR  EXPECTED 

RCT*EXPECTANCY  PER  YEAR,  IN  PERCENT,  OF  TT-  MINUTE  RAIN  RATE  (IN 
OATA  STATEMENT), 

DIMENSION  Al(6) ,A5(6> ,AE(F)  ,A7(6) ,Ae(6) ,B1(6) 

1,  82  (6),  B3  (6),  B4  (6),  BS  (El,  B6  (6),  TT  (6), 

2 TTR(i2),Sl(GI,S2(E ',S3(6) ,XX(12) 

COMMON/TEST/VRC  (12>,VRO(12) 

COMMON/RRATE/RR(12) ,VRR(12) ,PCT(12) 

. . .DATA  STATEMENTS  CONTAINING  COEFFICIENTS  FOR  MODIFIED  RICE- 
. . .HOLMBERG  (RH)  MODEL  RARAMFTEFS,  FOUND  IN  APPENDIX  C DF  DUTTON  ET 

. . .AL.  (1974). 

DATA  TT(1),TT(2) ,TT(3) ,TT(4) ,TT(5) ,TT(6)/1.,5. ,30, ,60,,360.,1440./ 
DATA  A1  (1),A1(2) ,A1(3) ,A1(4) ,A1 (5 ) , A1 ( 6 ) /. 0 091 33, . 0 12548, . 021934 , . 
AO  27642,. 10841,.  38794/ 

DATA  A5(1),AC(2)  ,A5(3)  ,A‘:(4)  ,A5(5)  , AS  (6) /O  . 0 ,-.  0157  82  ,-.  0084347  , 

AO  053345, -.0  01 8022, 0.0/ 

DATA  AE(1),A6(?) ,A6 (3),A6(4) ,A6 (5) , A6 ( 6) /O . 0 • 14. 31 3 , 7.4519 , 4.6 086 , 
Al. 7116, 0.0/ 

DATA  A7(1),A7 (2) , A7 (31 ,A7(4I ,A7(5) , A7 ( 6 1 /O . 0 , . 199  S3  , . 10638 , . 06705 3 

1..  022986. 0.0/ 

DATA  A8(1),A8(2I  ,A8(3),A8(4) ,A6(5> ,A9(6)/34.13?9,18.278,12.276,9.2 
1079,2.1793,1.00038/ 

DATA  B1(1),B1(2) ,B1 (3I,B1(4) ,B1 (9) ,B1 (6) /. 30045, . 30 049, . 32633, .491 
ASS, 1.3139, 1.8209/ 

DATA  62(1),B2(2) ,B2(3) ,B2(4) ,B2(5) ,B2 ( 6 ) /2 07. 05 , 2 07 . 0 5 , 224. 89 , 33 8 . 
884,904.10,1254.9/ 

DATA  B3(1),B3(2),B3(3),B3(4) ,B3(5) ,B3 (61 /3. 5329E-4 ,-6.6457E-4, 
l-,0013234,-4,9893E-4,-1.7799E-4,-6.262E-5/ 

DATA  B4(1),B4(2)  ,84 (3) .B4(4) ,B4(5) ,B4 ( 6) /. 24476 , 1 . 4 071, 2. 31 83 , 1. 14 

101. 0. 49868. .3022/ 

DATA  B5 (1),B5 (2) ,S5 (31 ,B5(4I ,B5(5) ,B9 (6 >/. 0 033902 ,. 01 6705 ,. 025179 , 

1.01194..  0045351.. 00 231 86/ 

DATA  86(11,86(21  ,86(31,86(41 ,B6 (9) ,B6 (6 ) /1 . 2807, . 44634, -. 17309, . 1 9 
A847,. 083186,. 1152/ 

DATA  SI (1),S1(2) ,S1 (3) ,S1 (4) ,S1 (5) .SI (6)/.2823,.0683,  .0411..0289, . 
A0060,0.C/ 


DATA  S2(1),S2{2) ,S2(3),r2(4) ,S2(5) ,S2(6)/25.0338,25.0338,27.1901,4 


oooo  onooo  ooon 


f fl3. 9671*, 109. 391, 151. 7197/ 

OAT  A S3(1»,S3(?)  ,S3<3),S3(«»),S3(5»  ,S3  ( 6) /.  0915  ,.  0915  ,.  0884  ,.  060U  , . 
A0ai,0,.0166/ 

NS=1 
NT  = 1 

IP  (NS  .LT.  1)  GO  TO  16'^ 

00  160  I = 1,  NS 

. . .05T5RMINE  PARAMETERS  OF  MODIFIFP  RH  MODEL  FROM  THE  ABOVE  COEF- 
. . .FTCIENTS. 

OE  = 0.C76561  • EM  - 93.632  ♦ BETA  ♦ 62.523 
VDE=5.8ei4E-3*VM*6,9943E3*V9ET+34.44*34.44 
IF  (NT  .LT.  1)  GO  TO  17P 
00  155  J = 1,  NT 
T2T  = B1  (J)  • EM  ♦ 92  (J) 

RIT  = Af  (J)  • EM  ♦ A6  (J)  • BETA  ♦ A7  (J)  * OE  ♦ A9  (JJ 

TIT  = BETA  ♦ EM  / RIT 

RPT  = 64  (J)  • BETA  ♦ 93  (Jl  • EM  ♦ B5  (J)  • OE  ♦ BE  (J) 

RC  = (RIT  • RPT  • ALOG  ((TIT  ♦ T2T)  / TITH  / (RIT  - RPT) 

OETERMINE  VARIANCES  OF  TH^  PARAMETERS  IN  THE  MODIFIED  RH  MODEL, 

AND  COMBINE  WITH  CERTAIN  PARTIAL  DERIVATIVES  (SEE  APPENDIX  0, 

DUTTON  ET  AL.  (1974)). 

VRP  = VM* (93(J) ♦B5( J) *0 . 076561 ) •*2*VRE T* (94 ( J ) -95 ( J ) * 8 3 .632 ) **2 
1«-B5  (J)*B5  (J)*34.44*34.4I.*S3  (J)*S3(  J) 

VT?=B1( J) •BK J)*VM  + S2( J)*S2(  J) 

VR1  = VM*  (A5(J)  *8  7 (J)  *0.07656 1)**2«’VBET*  ( A6  ( J ) - A7  ( J ) *93 .632 ) **2 
1*A7(J)*A7(J)*34.4<,*34.444S1  (J)*S1  ( J) 

VT1=( (BET A/RIT) **2. )*VM+ ( (EM/ PI T) **2 . ) • V9E T ♦ ( ( BET A* EM/ ( RIT’RIT ) )♦• 
12.) ’VPl 

RST  = 0.845  / (((30.  / RIT)  - (5.  / RPT))  ♦ ALOG  (1.  ♦ (T2T  / TIT) 
1)  ) 

X = (30.  • RPT  * RPT)  / (ALOG  (1.  ♦ T2T  / TIT)  • ((30.  • RPT  - 5. 

1*  RIT)  • (30.  • PPT  - 5.  • RIT))) 

Y = (5.C  • RPT  • RPT)  / (ALOG  (1.  ♦ T2T  / TIT)  ♦ ((30.  • RPT  - 5. 

1*  RIT)  * (30.  • PPT  - 5.  • PIT))) 

7 = (RIT  • RPT  • T2T)  / ((30.  * RPT  - 5.  • PIT)  • (TIT  • (TIT  ♦ T2 

IT)  • (ALOG  (1.  ♦ T2T  / TIT)  • • 2.0))) 

H = Z ♦ TIT  / T2T 

VRS= (X*X*VR1»Y*Y*VRP*Z*Z*VT1*H*W*VT2) •0.945*0. «45 
0 = EXP  ((2.340347319  / RST)  - (30.  / RIT)) 

TST  = TIT  * 0 

U = 30.  * TIT  • Q / (RIT  * RIT) 

V = 2.340347319  * TIT  * Q / (RST  * RST) 

VTS=Q*Q*VT1^U*U*VR1 ♦V*V*VPS 

OBTAIN  RAIN  RATES  AND  THRIP  VARIANCES  FDR  12  SELECTED  PERCENTAGES 

OF  AN  AVERAGE  YEAR. 

DO  150  K = 1,  12 
TTR  (K)  = 97. 6E  * PCT  (K) 

AL=AL06(TST/TTR(K) ) 

AL1=AL0G( (T1T*T2T) /TIT) 

AL2=ALOG( (T1T*T2T) /TT9(K)) 

9MR=BETA*EN/(R1T*R1T) 

“RR1T=AL0G(T1T/TTR(K) ) 

PRT1T=P1T/T1T 
T1=PRT1T*BETA/P1T 
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o o r>  rj  o o n o r>  n n rt  n 


i 


«0.967«*,  109. 361,151.71 97/ 

0«T«  S3(1),S3(?)  ,S3(3),S3(<»)  ,S3(5)  ,S3  1 6) /.  0915  , . 0915,  .0884  , . 0604  , . 
A0340,.0166/ 

NS»1 
NT  = 1 

IF  (NS  .LT.  1)  GO  TO  165 
00  160  I = 1,  NS 

. . .DETERMINE  PARAMETERS  OF  MODIFIED  RH  MODEL  FROM  THE  ABOVE  COEF- 
. . .FICIENTS. 

OE  = 0.C76561  • EM  - 83.632  * BETA  ♦ 62.523 

VOE =5. 8t 14E-3*VM46. 9943E3* VRET4 34. 44*34.44 

IF  (NT  .LT.  IJ  GO  TO  17C 

DO  155  J = 1,  NT 

T2T  * B1  (J)  • Em  ♦ 92  (J) 

RIT  » AE  (J»  • EM  ♦ A6  (J)  • BETA  ♦ A7  (J)  • OE  ♦ A8  (JJ 

TIT  = BETA  • EM  / RIT 

RPT  = B4  (J)  • BETA  ♦ B3  (J)  • EM  ♦ B5  (J)  • OE  ♦ BE  (J) 

RC  = (RIT  • RPT  • ALOG  ((TIT  ♦ T2TJ  / TIT))  / (RIT  - RPT) 


.BETERMINE  VARIANCES  OF  IN'-  PARAMETERS  IN  THE  MODIFIED  RH  MODEL, 
•AND  COMBINE  WITH  CERTAIN  PARTIAL  DERIVATIVES  (SEE  APPENDIX  0, 
.DUTTON  ET  AL.  (1974) ) . 

VRPcVM* (93 (J)  *95 (J) *0. 076561 )**2*VBET* (94 ( J ) -95 ( J ) • 8 3 .632) **2 
1*05  (J)*e5(J)*34.44*34.44*S3 (J(*S3( J) 

VT2=B1(J)*B1(J)*VM*S2(J)*S2(J) 

VRlsVM* (A5(J) *A7 ( J) *0 . 076561) *• 2* VBET* ( A6 ( J ) - A7 ( J ) • 83 .632) **2 
1*A7 (J)*A7 (J)*34.44*34.444S1 (J)*S1 ( J) 

VT1«( (BETA/R1T)**2. )*VM» ( (EM/ PIT) ••2. ) • V9ET* ( ( BET A* EM/ ( R1T*R1T ) )•• 
12.)*VR1 

RST  = 0.845  / (((30.  / RIT)  - (5.  / RPT))  ♦ ALOG  (1.  ♦ (T2T  / TIT) 
D) 


X 

3 (30. 

* RPT  * RPT) 

/ 

(ALDG  (1.  ♦ 

T2T 

/ TIT)  * ( (30. 

RPT  - 5. 

1* 

RIT)  * 

(30.  * PPT  - 

5. 

* RIT)  ) ) 

Y 

= (5.0 

* RPT  • RPT) 

/ 

(ALOG  (1.  ♦ 

T2T 

/ TIT)  * ((30. 

RPT  - 5. 

1* 

RIT)  * 

(30.  * PPT  - 

5. 

* RIT)  ) ) 

Z 

= (RIT 

* RPT  * T2T) 

/ 

((30.  * RPT 

- 5. 

• RIT)  * (TIT 

(TIT  ♦ T2 

IT)  • (ALOG  (1.  ♦ T2T  / TIT)  • • 2.0))) 

M = Z • TIT  / T2T 

VRS»(X*X*VRlTY*r*VRP»Z*7*VTl*M*H*VT2) *0.845*0. «45 
D = EXP  ((2.340347319  / RST)  - (30.  / RIT)) 

TST  = TIT  • 0 

U = 30.  * TIT  * (J  / (RIT  * RIT) 

V = 2.340347319  * TIT  * 0 / (RST  * RST) 

VTS  = Q*(J*VT1*U*U*VR1*V*V*VRS 


OBTAIN  RAIN  RATES  AND  THRIP  VARIANCES  FOR  12  SELECTED  PERCENTAGES 
OF  AN  AVERAGE  YEAR, 

DO  150  K = 1,  12 

TTP  (K)  = 87.66  * PCT  (X) 

AL=AL0G(TST/TTR(K) ) 

ALl=ALOG( (T1T»T2T) /TIT) 

AL2=ALOG( (T1T*T2T) /TT9(K)) 

9MR*BETA*EN/(R1T*R1T) 

“RR1T=AL0G(T1T/TTR(K( ) 

PRT1T*R1T/T1T 

T1=PRT1T*BETA/P1T 
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T2»PRT1T*EH/R1T 

T3=PRR1T-PRT1T*BMR 

PRRSTa«*.0*CRST**3.0)*<*L**l»,C) 

PRTST=PRRST*RST/«TST*#L) 

PRSTRP*-0.8i*5*5.0*RlT*RlT/( ( J30.*RPT-5.*R1T)**2,0) •ALl) 

HixPRSTRP* (93 U» 485  076561) 

PRSTT1*0.8<i5*T2T/(T1T*(T1T4T2T)*»L1*»L1*(30./R1T-5./RPT)) 

PRSTT2*-PRSTT1*T1T/T2T 

T6=PRSTT1*BCT»/R1T 

T7*PRSTT2*B1(J) 

PTSTTlxQ 

T8=PTSTT1*BETA/R1T 

T10=PRSTT1*BFTA/R1T 

TllxPRSTT2*Bl(J) 

T9=PRSTRP*(B3(J)495 ( J ) *0 .07656 1 ) 

“TSTRS=-W 

T12=PRSTRP*(91»(J)  -B5(  J)  *83.632) 

T13*PRSTT1*EM/R1T 
T1«.  = PTSTT1*EM/R1T 
Tl5rPRSTPP*{94(J) -B5( J) *83.632) 

T16xPRSTTl*EM/RlT 

Tl7*PRSTPP*9‘:(J) 

PRSTR1=-6.0*(RFT/R1T)*(RPT/R1T)*PRSTRP 

ri8*PR6TTl*BBR 

TiqsPTSTTl*B1R 

720x118 

“TSTRlxU 

PRRPTXAL2 

T21=PRRPT*(93(J) *95 ( J) *0.076561) 
oRT2TxPPT/(T1T4T2T) 

T22*PRT2T*BETA/R1T 

T23=PRT2T*B1(J) 

T24xPRRPT*(B4(J)-B5 ( J)*83.632> 

T25»PRT2T*EM/R1T 

730  X 717  * EXP  ( - 30.  / R17) 

76  X (717  ♦ 727)  * FXP  (-<=./  RP7) 

7C  = 717  * EXP  ( - RC  / R17) 

RS  = - 0.844996  / (ALOG  (730  / 75 > ) 

7S  X 730  * EXP  ((30.  * * 0.25)  / RS) 

RI  X ( - RS  * ALOG  (77R  (K)  / 7S) ) * * 4.0 
IF  ((RI  .L7.  30.0)  .AND.  (PI  .G7.  5.0))  GO  70  100 

9G  X - RiT  * ALOG  (77R  (K)  / 717) 

9L  X - RPT  * Al  OG  (7TR  (X)  / (717  ♦ 727)) 

GO  70  106 
100  RLxRl 
RGxRI 

105  IF  (77R  (K)  .G7.  70  GO  70  110 
RR  (K)  X RG 
GO  70  115 
110  RR  (X)  X RL 
115  C0N71NUE 

IF  (77  (J)  .G7.  60.0)  GO  70  125 
IF  (RR  (K)  .GE.  5.0)  GO  70  120 
GO  70  130 

120  IF  (RR  <K)  .07.  30.0)  GO  70  135 
GO  70  140 

125  IF  (RR  (K)  .G7.  RC)  GO  70  135 

130  VRC  (K)  X ((V71  ♦ V7’)  * (RP7  / (717  ♦ 7 27 ) ) **2 . ) ♦ ( VRP* ( ALOG ( ( 71 7* 
1 727)/77»(X)))**2.) 


1 ♦ < IPRRPT*B5  « J)  I ••2.0»*3«».4«.*34.<*4*»PRT2T*BMR 

2 •PRT2T*eHR*SHJ»*SHJ»»PRT2T*PRT2T*S2(J>*S2<J) 

3 ♦PRRPT*PRRPT*S1CJ»*S3fJ) 

GO  TO  lli9 

135  VRC  «K>*WT1*HR1T/T1TI  ••2.)  ♦ CVRl  *(*100  (TIT  / TTR  (K»n**2.) 

tfRO  (K»«T1*T1*VM^T2*T2*VBET^T3*T3*S1(J)*S1(J» 

GO  TO  I'tS 

140  F*C»«*.0*(  (RST*AL0G(TST/TTR(K)  I )**4.0» 

VRC(K)*I(FAC/RST)**2.>*VRS*( (F»C/(TST*»L0G(TST/TTR(K) n)**2.)*VTS 
WRO  (K)  = ((PRRST*(T(»*T6*T7  » ♦PRTST*  (T8*PTSTRS*  ( T9*T  1 0 ♦Tlin  ) 

1 •♦2.0)*VM4(  (PRP<;t*(T12*T13)*PRTST*(T14 

2 ♦PTSTRS*(T15>T16» » »**2,0)*VBeT*( IPRRST*T17 

3 ♦PRTST*PTST9S*T17)**2.0)*34.44*34.i.4 

4 ♦(  (PRRST*(PRSTR1-T1«) ♦PRTST* ( -TlB^PTSTRl 

A ♦PTSTRS*(PRSTR1-T20))»»*2.0I*S1(J)*S1(JK(  (PRRST* 

5 PRSTT2^PRTST*PTSTRS*PRSTT2)**2.0»*S2(J) *$2 ( JJ 

6 ♦ ((PRRST*PRSTRP*PRT5T»PTSTRS*PRSTRP) ••2. 0 ) *33 ( Jt • S3 (J) 

145  CONTINUE 

150  CONTINUE 
155  CONTINUE 
170  CONTINUE 
16C  CONTINUE 
165  continue 
C 

C.  . .FIT  THE  RAIN  RATE  VARIANCE  RESULTS  MITH  A SMOOTH  CURVE  FOR  PREOIC- 
C.  . .TION  PURPOSES.  AND  PERFORM  THE  PREDICTION  FOR  THE  12  PERCENTS  OF 

c.  . .an  average  YEAP. 

C 

00  200  Isl.12 

XX(II*AL0G(AL06( (TlT*0.5*T?TI/TrR(r>  > » 

VRO(I)*ALOG(VRO(I)  ) 

200  CONTINUE 

CALL  LINREG(12«XXtVR0«XBAR.YeAR,8t A, VARX.VARY.COV.SOX.SOr.SER.Z.R) 
AE=EXP(A) 

00  275  I»l,12 

VRRtI)  = AE*(  (EXP  (XX  ( in  )**8I 
275  CONTINUE 
RETURN 
END 
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SUBROUTINE  LlNREC  IN,  X,  Y,  XBAR,  YBAR,  B,  A,  VARX,  VARY,  COV,  SOX 
t,  SOY,  SER,  Z,  R> 

perforns  a standard  linear  regression  iy*bx«ai  fitted  through  the 

MEAN. 

INPUT 

N - NUHBER  OF  DATA  POINTS 
X • INOEPENOeNT  VARIABLES. 

Y - DEFENOENT  VARIABLES. 

OUTPUT 

XBAR  * ESTIMATED  MEAN  OF  X. 

YBAR  - ESTIMATED  MEAN  OF  Y. 

B - SLOPE  OE  FITTED  LINE. 

A - INTERCEPT  OF  FITTED  LIME. 

VARX  - ESTIMATED  VARIANCE  OF  X. 

VARY  - ESTIMATED  VARIANCE  OF  Y. 

COV  - ESTIMATED  COVARIANCE  OF  X AND  Y. 

3DX  - ESTIMATED  STANDARD  DEVIATION  OF  X. 

SOY  - ESTIMATED  STANDARD  DEVIATION  OF  Y. 

SER  - ESTIMATED  STANDARD  ERROR  OF  ESTIMATE. 

Z - ESTIMATED  FISHERS  Z. 

R - ESTIMATED  CORRELATION  COEFEICIENT. 


1 o: 

105 


DIMENSION  X 
SUMX=0. 
3UMY=0. 
SOMXXsC. 


(ICOG)  , V (lOOw) 


SUMVYsO. 

SUMXYsO. 

:n  = N 

IF  IN  .LT.  U 

GO 

TP  1U5 

DO  IjO  1=1, 

N 

SUMX  = SUMX  ♦ 

X 

(I) 

SUMY  = SURV  ♦ 

Y 

(I) 

SUMXV  * SUMXY 

* 

X m • 

Y 

(I) 

SUMXX  = SI'MXX 

♦ 

X (I)  • 

2 

SUMYY  = SCMYY 
CONTINUE 

♦ 

Y 1 1)  • 

2 

X3AR  = SURX  / 

EN 

Y9AR  = SUt-Y  / 

IN 

VARX  = ISIMXX 

- 

X3AR  • 

SUMX)  / 

(EN 

VARY  = ISliMYY 

- 

Y3AR  • 

SUMY)  / 

(EN 

COV  = (SUMXY  - 

X8AR  • SUMY)  / 

lEN 

1.) 


3 = COV  / VAi-X 
A = YBAR  - B • X8AR 
SOX  = Srj?T  (VARX) 

SOY  = SORT  IVARrT 
R = COV  / (SDX  • .PY) 

Z = .5  * ALSO  1(1.  ♦ R) 
S'R  = SORT  KEN  - l.DI  ‘ 
RETURN 
c TO 


/II. 

VARY 


- R)  ) 

• (1.0 


R)  / (EN 


Z.i)) 
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c 

c 
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FUNCTION  6«MH«(FREQ» 

FREQ  IS  FREQUENCY  IN  GHZ  BETWEEN  7 GHZ  ANO  100  GHZ. 

OINENSION  F(Z7I ,G(Z7) 

THESE  OATA  ARE  TAKEN  FROM  THE  CCIR  CURVES  OF  ATTENUATION  PER  KM. 


GAMMA  IS  RAIN  ATTENUATION  COEFFICIENT  AT  FREQ  IN  DB/KM/MH/HR. 

GAMHAsTERPCZS.FREQ.F.G) 

RETURN 

END 


ooonoonoo 


FUNCTION  RrOCO(D<T) 


THIS  FUNCTION  OBTAINS  THE  REDUCTION  FACTOR  FOR  POINT  RAIN  RATES 
USED  BY  BATTFSTI  ET  AL.  11971), 


I 

1^ 


INPUT 

0 - PATH  LENGTH  IN  (CM. 

T - PERCENT  OF  TIME  OF  OPERATION  OF  CONCERN, 

OIHENSION  X(12) .Y(12) ,9(3) 

DATA  X(1),XI2),X(3),X(4),X(S),X(6I,X(7),XI8),X(9I,X(10),X(11),XI12 
A) /O • ,2,5,9, ,7,5 ,10 , ,15, ,20, ,25, , 30, ,35, ,40, ,50,/ 

BATA  Yd)  ,Y<2),Yt3)  ,Y(4),VI5),Y(6>  , Y (7)  , Y(0>  , Y (9)  , Y dO)  , V (11)  , Y d2 
A) /I, ,,91, ,86, ,81, ,78,, 72, ,67, ,62,, 58, ,57, ,54, ,5/ 

20  Rd)»l 

R«2)»10,**«-, 0036*0) 

IF(D,GT,25,)  R(2)»10.**(-,003*C-,05) 

R(3)sTERPd2,0,X,Y) 

IFd.LT.,!)  REOCOsI  (T-,01)/.09)*(R(2)-R(3))*R(3) 

IF«T,GE,,1)  PEDCO=( (T-,1 ) /, 9) • (P d )-R ( 2) ) ^R  (2) 

RETURN 

END 


<;uaROUTINr  RAINRT(ROtRSS«RBI 


r 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  SUBROUTINE  MODIFIES  THE  RAIN  RATE  IN  ACCORDANCE  HITH  THE 
PROCEDURE  OF  BARSIS  ET  AL.  <19731. 

INPUT 

RO  • UNHOOIFIEO  IPOINT)  RAIN  RATE  IN  HN/HR. 

RSS  - PATH  LENUTH  IN  KH. 

OUTPUT 

RB  - MOOTFIEO  IPATH)  RAIN  RATE  IN  MH/HR. 

OINENSION  RATI3)tPL(3l 
IFfRO.LE.lO.O)  229.230 

229  RTB«1.0 
GO  TO  290 

230  RAT  (!)»-. 0907679«»672*ALOGIR0)*J. 209 
PL(ll>S.O 

RAT (21 «>.18B91B0996*ALO6<R0> 41.439 
PL(2>>10.0 

IFCRQ.GT. 20. 0)239. 243 

239  RAT (3 )>>.13O7074S21«ALOG(R0 I 41. 0 36771423 
GO  TO  249 

240  RAT <3 )«-. 3999940639 *4106(90) 4 1,911717924 
249  PL(3)s22.0 

RSSSsRSS 

IF(RSSS.GT.22.0)  RSSS*22.0 
RTB*TERP(3,RSSS.PL,RAT) 

290  RB*RTB*R0 
RETURN 
END 


oooooooooon 


SUBROUTINE  ATCOS<FRtT,P,RHtHTOf>,PROBtBETA«R«ATEtMAV) 

INPUT 

PROB  - PERCENT  OF  AN  AVERAGE  YEAR  FOR  NHICH  RESULTS  ARE  DESIRED 
RATIO  OF  THUNDERSTORM  RAIN  TO  NON-THUNDERSTORM  RAIN. 

R - RAIN  RATE  CORRESPONDING  TO  PROB  AT  THE  SURFACE  OF  THE  EARTH 

OUTPUT 

TOTAL  ATMOSPHERIC  ATTENUATION  PER  UNIT  LENGTH  IN  OB/KM, 


PI-3.1A16 

TK*T*273.16 

CALL  RATTCO(NAV,TK,P,RH«HTOP,0.2.R.AT,PROB,BETA,SCAT,PHIRR> 
SRMOa2lE.68*ESUBSITK)*RH/TI{ 

CALL  OXYGENITKtP.SRHO.MAV.RMAV.OXAT) 

CALL  MATER (TKfP,SRHO.MAV,MN«NATl«NAT2.PHM) 

IFIFR.LE.AS.OI  GO  TO  200 

CALL  TOPOXYfTK.P.MAV.l.OXAT.SRHOtGAMAUP) 

200  ATE=2.0*AT*OXAT»MAT1*MAT2 

CALL  REFRAC(<»t3t6tTtP«RH,ENS«O0  .MO  ) 

«»HIRD=PHIRR^PHN 

PHIT*l2,lFtS*PI/MAV>*(1.0»1.0E-6*ENS)^PHIRO 

RETURN 

END 
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SUBROUTINE  CHPLXN  IHAVEt  T«  CSO«  HIHHKt  CFPT,  CTPT,  Ct  0> 

THIS  SUBROUTINE  COMPUTES  THE  EFFECT  OF  THE  DIELECTRIC  CONSTANT  OF 
hater  in  CLOUO  attenuation  COMPUTATIONS. 

INPUT 


HAVC>NAVELENGTH  IN  CENTIMETERS. 
TaTEMPERATURE  IN  DECREES  CENTIGRADE. 

OUTPUT 


CSQaOIELECTRIC  INFLUENCE  OF  HATER  ON  SCAT  IN  ATTCOE. 
HIMMKxOIELECTRIC  influence  OF  MATER  ON  ABS  IN  ATTCOE. 
CFPT  NOT  USED  IN  TEST30. 

CTPT  NOT  USED  IN  TEST30. 

CaREAL  PART  OF  THE  DIELECTRIC  COEFFICIENT  OF  HATER. 
OalNAGINARV  PART  OF  THE  DIELECTRIC  COEFFICIENT  OF  HATER. 


DIMENSION  EDO  (6).  OOLAM  16),  TTAB  (6) 

TTAB  (1)  a 0.0 
EDO  (1)  a 83.0 
OOLAM  ID  a 3.59 
TTAB  12)  * 10.0 
EOO  (2)  a OA.O 
DOLAN  12)  a 2,2k 
TTAB  IS)  a 10.0 
EOO  13)  a 81.0 
DOLAN  13)  a 1.66 
TTAB  (A)  a 20.0 
EOO  (A)  a 80.0 
OOLAM  lA)  a i,5J 
TTAB  (5)  a 30.0 
EOO  15)  a 76. A 
OOLAM  15)  a 1.12 
TTAB  (6)  a 40.0 
EOO  (6)  a 73.0 
OOLAM  16)  a 0.059 
EIN  a 5.5 

EO  a TERP  (6,  T,  TTAB,  EOO) 

OLAM  a TERP  (6,  T,  TTAB,  OOLAM) 

3 a OLAM  / have 

A a lEO  - EIN)  / (1.0  ♦ IB  • • 2)) 

C a EIN  - 1.0 
0 a EIN  ♦ 2.0 
E a A ♦ C 
H a A ♦ 0 
G a A • B 

HIMMK  a IG  ♦ CH  - E))  / I (H  • • 2 ) ♦ I G • • 2)) 

CSQ  a <(CE  • H ♦ (G  ••  2))  / (IH  ••  2)  ♦ (G  ••  2)))  ••  2)  ♦ HIMMK 
1 ••  2 
C a A ♦ EIN 
0 a - A • B 
U a 0 / 45. 

V a C(20.  • C • D ♦ 13.  • 0)  • (15.  a (C  • ♦ 2)  - 15.  • (0  • • 2) 

!♦  60.  a C ♦ 60.)  ♦ (30.  a c a c a 60.  a q)  a ( - i q . a (c  a a g)  « 

2 10.  a (0  a • 2)  - 10.  a c ♦ 20.))  / ((15.  a (c  a ♦ 2)  - 15.  a <o 

3a  a 2)  ♦ 60.  a c ♦ 63.)  a • g ♦ (30.  • C a 0 ♦ 60.  a Q)  a • g) 


u u uu 


3*  21  « 60.  • C ♦ 60.)  ••24  130.  • C • 0 ♦ 60.  • 0)  • • 21 
T » <0  • (30.  • C * 65.)  - 30.  • D • (C  - 1.))  / ((30.  *04  65.) 

• 2 4 (30.  • 0)  • 4 2) 

CFPT  >3.  *04  3.  41(4  5.  •Y 

CTPT  « 3.  4 V 

RETURN 

END 


FUNCTION  ESUBS  (T) 

THIS  SUBROUTINE  CALCUCRTES  SATURATION  YAPOR  PRESSURE  OF  HATER  IN 
AIR  AT  TEMPERATURE  T IN  QEGREES  KELVIN. 

X * .05  4 (T  - 263.) 

y * 26.661775  - X • (.336222  - (X  - 1.)  • (.95e9E-2  - (X  - 2.)  • r 
1.166666E-3  4 (X  - 3.)  4 .229E-6))) 

ESUBS  * EXP  (Y  - 659'.. 6076  / T) 

IF  (T  - 283.)1C5,  105.  100 
100  ESUBS  » ESUBS  4 (100.  * (T  - 293.)  • • 2)  • 8.E-6 
105  RETURN 
END 
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SUBROUTINE  REFRAC  (KT«  KP.  KHt  TN«  PR«  HUt  FN*  Ot  W) 

SPECIFIC  LOCATIONS  (HEICHTS)  IN  THE  ATMOSPHERE  FOR  VARIOUS  KINDS 
THIS  SUBROUTINE  CALCULATES  THE  REFRACTIVITY  (REFRACTIVE  INDEX)  AT 
OF  METEOROLOGICAL  INPUT  PARAMETERS. 

INPUT 

KT3I»«  means  TEMPERATURE  INPUT  IN  DEGREES  CENTIGRADE. 

KP«3*  means  pressure  INPUT  IN  MILLIBARS. 

KH=6«  means  RELATIVE  HUMIDITY  INPUT  AS  A DECIMAL  FRACTION. 
TN«TEMPERATURE  IN  DEGREES  CENTIGRADE. 

“R»PRESSURE  IN  MILLIBARS. 

HUsRELATIVE  HUMIDITY  AS  A DECIMAL  FRACTION. 

OUTPUT 

EN»REFRACTIVITY  in  N-UNITS  (PARTS  PER  MILLION  OF  REFRACTIVE  INOEXi 
OsORY  TERM  OF  REFRACTIVITY  IN  N-UNITS. 

N*NET  TFRM  OF  REFRACTIVITY  IN  N-UNITS. 


T a TM 

P a PR 

H a HU 

60  TO  (100,  105,  110, 

115,  130),  KT 

100 

T a T - 459.69 

105 

T a .555555555  • T - 
GO  TO  115 

17,7777777 

110 

T a 1,25  a T 

115 

T a T ♦ 373, 

130 

GO  TO  (135,  130,  135) 

, KP 

135 

o a 25,4  • P 

130 

P a 1,333334  • P 

135 

GO  TO  (140,  140,  175, 

150,  185,  190, 

195),  KH 

140 

GO  TO  (145,  150,  155, 

160,  165),  KT 

145 

H a H - 459.69 

150 

H a .555555555  • H - 
GO  TO  160 

17.7777777 

155 

H a 1.35  • H 

160 

H a H ♦ 373. 

165 

E a ESUeS  (H) 

GO  TO  (315,  170) , KH 

170 

P a E - (,66E-3  a (i, 
GO  TO  315 

♦ .115E-3  • (H 

- 373,))) 

175 

H a H / (1.  - H a i.F 

-3) 

180 

H a H • l.E-3 

X = T - 373. 

FN  = l.OOORL  - X • (.33F-fc  ♦ .175E-6  • X)  ♦ (.39E-5  ♦ X ♦ (.45E-9 
1*  X - .3B5E-7))  • P 
E « H • P / (FH  • (H  ♦ .E3197)) 

GO  TO  315 

185  E » 4.61«:0136  • H • T 
GO  TO  315 

190  E » M • ESU9S  (T) 

GO  TO  315 

195  GO  TO  (300,  305,  310),  KP 
300  H = 35,4  • H 

305  H s 1.333334  • H 

310  E = M 

315  N « 373356.  • E / T • • 3 

0 « 77.6  ♦ P / T 

EN  a M ♦ 0 

RETURN 

END 
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SUBROUTINE  OXYGEN  (Tt  R«  RHOt  PLANDA*  RMAV«  6ANAUI 

c this  subroutine  calculates  atnospheric  absorption  per  unit  length 

C OUE  TO  ATNOSPHERIC  OXYGEN  FOR  FREQUENCIES  ROUGHLY  IN  THE  RANGE 

C GHZ.  USES  VAN  VLECK  FORNULATION  (SEE  RADIO  HETEOROLOGY.  PAGE 

C Z72>. 

C 

C INPUT 

C 

C TsTEHPERATURE  IN  DEGREES  KELVIN,  AT  LOCATION  ON  TRANSHISSION  PATH 

C OF  INTEREST. 

C RsPRESSURE  IN  HILLIBARS  AT  LOCATION  ON  TRANSMISSION  PATH  OF 

C INTEREST. 

C RHO*ATMOSPHERIC  MATER  VAPOR  DENSITY  IN  G/M(3),  AT  LOCATION  ON 

C TRANSMISSION  PATH  OF  INTEREST. 

C PLAHOA>NAVELENGTH  IN  CENTIMETERS. 

C 

C OUTPUT 

C 

C RMAVsRECIPROCAL  OF  WAVELENGTH. 

C GAHAU»OXYGEN  ABSORPTION  IN  OP/KH. 

C 

PM  = (<f.615E-3)  * RHO  ♦ T 
P = P - PM 

D = O.OAS  • (P  / 1013. 2E)  • ((300.  / T)  • • .75> 

RMAV  si./  PLAHOA 
RMAV2  s RMAV  • RMAV 
02  = 0 • 0 

FI  = D / (RMAV2  ♦ 02J 

F2  = 0 / ((RMAV  - 2.)  • • 2 ♦ 02) 

F3  = 0 / ((RMAV  *2.)  • • 2 ♦ C2) 

GAMAU  = .34  * RMAV2  * ((293.  / T)  • ♦ 2)  ♦ (FI  ♦ F2  ♦ F3)  • (P  / 1 
1013.25) 

RETURN 

END 
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SUBROUTINE  MATER  (T«  Pt  RHO*  MAVE*  MNt  CAHARt  GAMANR,  PHM) 

This  subroutine  calculates  atmospheric  absorption  per  unit  length 

OUE  TO  ATMOSPHERIC  MATER  VAPOR  FOR  FREQUENCIES  ROUGHLYIN  THE  RANGE 
5-70  GHZ.  USES  VAN  VLECK  FORMULATION  (SEE  RADIO  METEOROLOGY,  PG. 
2721. 

INPUT 

TsTEMPERATURE  IN  DEGREES  KELVIN,  AT  LOCATION  ON  TRANSMISSION 
PATH  OF  INTEREST, 

P=PRESSURE  IN  MILLIBARS,  AT  LOCATION  ON  TRANSMISSION  PATH  OF 
INTEREST. 

RHO=ATMOSPHERIC  MATER  VAPOR  TENSITY  IN  G/M(3>,  AT  LOCATION  ON 
TRANSMISSION  PATH  OF  INTEREST. 

MAVE*MAVELENGTH  IN  CENTIMETERS. 

OUTPUT 

MNsRECIPROCAL  OF  WAVELENGTH. 

GAMARsCONTRIBUTION  OF  MATER  RESONANCE  LINE. 

GAMANR«NON-RESONANT  CONTRIBUTION  TO  MATER  VAPOR  ABSORPTION  IN 
OB/KM. 

PHM=PHASE  DISPERSION  DUE  TO  HATER  VAPOR  IN  RADIANS/KM. 

F = 2B.R793  / HAVE 
FO  = 22.23515 
OF  = FO  - F 
Cl  - 0.00361 
C?  = 0.06089 

A = 0.08478  • (P  / 1013.25) 

B = 0.00708 
WN  = 1.  / MAVE 
WNR  = 0.7417 
X = 318.  / T 

D = A • (1.  ♦ B • RHD)  • (X  • • 0.625) 

FORMM  = (MN  - MNP)  •♦2*0**2 
FORMM  = 1.  / FORMM 

FORMP  = (MN  ♦ MNP)  **2+0**2 
FORMP  X 1.  / FORMP 


FORM  X C • (FORMM  ♦ FOPMP) 

TT  X (X  • ♦ 2.'^)  • EXP  ( - » (1,  / T - 1.  / 315.)) 

YP  X Cl  ♦ TT  • FORM  • (MN  • * 2) 

YNR  X C2  • 0 • X • (MN  • • 2) 

GAMAR  X YR  • RHO 
GAMANR  X YNR  • PHO 

GGxO. 0186823*RHO>0. 00 28129* (P-PH0*T/216.68)*( (300. 0/T)* *0,63) 
PHMx (4.19168E-2)*F*0. 110132*DF/(OF*0F+GG*GG) 

RETURN 

END 
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SUBROUTINE  RATTCO ( N A VE , TRtPRE , REL , HITE t ZHH » RO , AT , RE LV |BET4 , SC ATP , 
IPHIR» 

C THIS  SUBROUTINE  CALCULATES  RAIN  ATTENUATION,  PROVIDED  AN 

C AIR/GROUND  PATH  IS  USED  FROM  SURFACE  METEOROLOGICAL  DATA  AND 

C RELIABILITY  REQUIREMENTS. 

C 

C INPUT 

C 

C HAVE*NAVELENGTH  IN  C'=’NTI METERS. 

C TA=SURFACE  TEMPERATURE  IN  DEGREES  KELVIN. 

C PREaSURFACE  PRESSURE  IN  MILLIBARS. 

C REL=SURFACE  RELATIVE  HUMIDITY  AS  A DECIMAL  FRACTION. 

C HITEsSTORM  TOP  HEIGHT  IN  KILOMETERS. 

C ZHHaHEIGHT  ABOVE  EARTHS  SURFACE  IN  KILOMETERS. 

C RO=SURFACE  RAINFALL  RATE  IN  MILLIMETERS  PER  HOUR. 

C RELYsPERCENT  OF  AN  AVERAGE  YEAR. 

C BETAaRATIO  OF  THUNDERSTORM  RAIN  TO  NON-THUNDERSTORM  RAIN  (SEE 

C EARLIER  REFERENCE  IN  TEST30). 

C 

C OUTPUT 

C 

C SCATRsRAINFALL  REFLECTIVITY  IN  KM(-l). 

C PHIRaPHASF  DELAY  PER  UNIT  LENGTH  DUE  TO  RAIN,  IN  RAOIANS/KM. 

C ATaATTENUATION  COEFFICIENT  DUE  TO  RAIN  IN  DB/KH. 

C 

DIMENSION  PF  (3),  FACT  (3) 

DATA  PF(1»,PF(2» ,PF (3)/.01,.l,l./ 

FR  r 29.9793  / HAVE 

CALL  SFCG  (TA,  PEL,  PRE,  HITE,  ZHH,  RO , EMZl,  EMZ2,  BETA,  PELY , UN 
IN) 

CALL  CRANE  (FR,  A,  8,  Al,  Bl) 

TAC  a TA  - 273,16 

CALL  CHPLXN  (NAVE,  T»C,  CSQ,  HIMMK,  CFPT,  CTPT,  CC , DO) 

Fo  a 0.P25970  • FR  - 0.r013S 
IF  (FR  .LE.  19.31)  FP  = 0.0 
FQ  a 0.01A8889e89  • FR  - D.402 
TF  (FR  .LE.  27.0)  FP  = 0.0 
FACT  (1)  a 1.0  ♦ FP 

fact  (2)  a 1.0  ♦ FQ 

FACT  (3)  a 1.0 

FftC  a TERP  (3,  RELY,  PF,  FACT) 

TF  (RELY  .GT.  1.0)  FAC  = 1.0 

U a (0.3  * 3.1415927  / NAVE)  • (FAC  • ( (CC  - 1.)  • (CC  ♦ 2.)  ♦ DO 
1*  DO))  / ((CC  ♦ 2.)  • (CC  ♦ ?,)  ♦ DO  • OD) 

OHIR  a U • EMZl 


AT 

a A • (EMZl  • 

• B) 

IF 

((ABS  (RELY  - 

O.D)  .LT. 

0.01) 

AT 

= 

BETA  • A • (EMZl  • • 9)  ♦ 

1(1. 

- BETA)  • A 

• (EMZ2  • 

* 9) 

IF 

((ABS  (RELY  - 

0.1))  .LT. 

C.Cl) 

PHIR 

a beta  • U • EMZl  ♦ (1.  - 9 

l^’TA)  • U * EMZ2 

IF 

((RELY  - 0.1) 

.GP.  0.01) 

AT  a 

A ♦ 

(EMZ2  • • B) 

TF 

((RELY  - 0.1) 

.GE.  0.01) 

PHIR 

a U 

« 

EMZ2 

C a 

(2.E-7  / 3.) 

* (306.01812)  • 

CSQ 

/ 

(NAVE  • • 4) 

X a C • Al  • (EMZl  • • 91) 
Y a C • Al  * (FMZ2  • • PI) 
E a 0.64  • 81 
SCATR  a Y 
10?  RETURN 
END 
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SUBROMiNE  SFCG  IT»,  REH,  PR,  HOT,  ZH,  RO,  EMZl,  EMZ2,  BETA,  RELY, 
1 Tl 

THIS  SUBROUTINE  CALCULATES  THE  LIQUID  MATER  CONTENT  IN  A RAIN> 
STORN.  DETAILS  ARE  IN  A METEOROLOGICAL  MODEL  FOR  USE  IN  THE  STUDY 
OF  RAINFALL  EFFECTS  ON  ATMOSPHERIC  RADIO  TELECOMMUNICATIONS,  BY 
E.J.  OUTTON,  OFFICE  OF  TELECOMMUNICATIONS  REPORT  OT/TRER  2A, 
PPCEhBER,  1971. 

INPUT 

TA=SURFACE  TEMPERATURE  IN  DEGREES  KELVIN, 

reh»surface  relative  humidity  as  a decimal  fraction. 

PR=SURFACE  PRESSURE  IN  MILLIBARS. 

HGT=STORM  TOP  HEIGHT  IN  KILOMETERS. 

zh=height  above  earths  surface  in  kilometers. 

RiJsSURFACE  RAINFALL  RATE  IN  MM/HR, 

9ETA=RATI0  of  thunderstorm  rain  TO  NON-THUNDERSTORM  RAIN  (SEE 
“RFVIOUS  REFERENCE). 

RELYsPERCENT  OF  AN  AVERAGE  YEAR. 

OUTPUT 

EMZl*CONTRIBUTION  IN  G7M(3)  OF  CONVECTIVE  RAINSTORMS  TO  LIQUID 
MATER  CONTENT, 

EMZ?=CONTRIBUTION  IN  G/Mt3)  OF  STRATIFORM  RAINSTORMS  TO  LIQUID 
HATER  CONTENT. 

REAL  L 

Z = 1000.  • ZH 
HIT  = 1000.  • HGT 
Cl  = 1.9031 
C?  = 1.5625 

TO  = TA  / (1.  - 1.8594E-‘*  • TA  • ALOG  (REH)  ) 

L = (123.  ♦ .227  * (TO  - 273.16))  • (TA  - TO) 

HTE  = ((RO  / Cl)  ♦ ♦ (1.  / C2))  ♦ (2.  • L / 1852.0) 

IF  ((Z  / 1852.0)  .LT.  HTE)110,  105 
105  EMZl  = C.O 
GO  TO  115 

no  RR  s Cl  • (((HTE  - (L  / 1852.0))  - ABS  ((Z  - L)  / 1852.0))  * • C2) 
9Z  = 8.2  • (PR  • • ( - .21)) 

■^MZl  = 64,  • 3.1415927  • (BZ  • • ( - 4)» 

115  9 = 8.2  ♦ (RO  • • ( - .21)) 

EML  s 64.  * 3.1415927  ♦ (B  • * ( - 4)) 

IF  (Z  .LE.  DlOO,  120 
100  EMZ2  = EML 

T s TA  - Z • 9.8E-3 
GO  TO  125 

120  TL  = TA  - L • “.8E-3 
E = ESUeS  (TL) 

D - PR  • FXP  ((980.62  / 28704.)  ♦ ( - LI  • (2.  / (TA  ♦ TL))) 

M = (E  * ,622)  / P 

TS  = 9.8E-3  • ((1.  ♦ (M  • 597.3  / (TL  • 6.8557E-2)))  / (1.  ♦ ( (M  • 
1 22.191E4)  / (TL  • TL  • 1.64537E-2 ) ) I ) 

T = TA  - L • 9.8E-3  - (Z  - L)  • TS 
EXZ  = EXP  ( - .064  • TS  • (Z  - D) 

=^XH  = EXP  ( - .064  • TS  • (HIT  - D) 

T9Z  = TL  - .5  • TS  • (Z  - L) 

TBH  s TL  - ,5  ♦ TS  • (HIT  - L) 

G s (.622  / (T  • 2.8704E-3))  • I - .064  • TS  • E • EXZ  ♦ (980.62  • 
1 F • EXZ  / (T  • 2.8704E4))) 


36 


R * 2.e7QliE-3 
R1  « 2.870i»E4 


GAH2  s 

(1.622 

E 

/ 

(R 

• TB2n  • (EX2 

- 1.1 

- (.622 

» 980.62  • 

E 

/ 

1 I.06N 

• TS  • 

R 

R1 

• 

TB2  • TB2n  • 

(EX2  - 

l.)l 

GANN  « 

((.622 

E 

/ 

(R 

• TBH»)  • (EXM 

- l.J 

- (.622 

• 980.62  • 

E 

/ 

1 (.06<i 

• TS  * 

R 

R1 

TBH  • TBH))  • 

(EXH  - 

i.n 

PH22  » EML  • (1.  - (GAH2  / GAMHI) 
IF(EH22.LT.0.0)  EH2230.0 
126  CONTINUE 
RETURN 
END 


i 
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SUBROUTINE  CRANE  (F«  A«  St  Ct  D> 

THIS  SUBROUTINE  GETS  THE  COEFFICIENTS  A ANO  B IN  A«(H»*ei  AND 
C ANO  D IN  C*(N**OI  FOR  (INPUT)  FREQUENCY*  Ft  BELOM  94  GHZt 
FOLLOWING  CRANE*  R.K.t  HICROHAVE  SCATTERING  PARAMETERS  FOR  NEM 
ENGLAND  RAIN*  HIT  LINCOLN  LABORATORIES  REPORT  R TR  426*  AO  647798* 
OCT.  1966. 

THE  OUTPUTS  ARE  A ANO  B*  COEFFICIENTS  IN  COMPUTING  ATTENUATION 
COEFFICIENT*  C ANO  0*  COEFFICIENTS  USED  IN  COMPUTING  REFLECTIVITY 
PER  UNIT  VOLUME.  M IS  LIQUID  MATER  CONTENT  IN  G/M(3)*  BUT  IS  NOT 
OUTPUT  OF  THIS  ROUTINE.  THE  SUBROUTINE  USES  THE  APPROXIMATE  F 
SQUARED  ABSORPTION  OEPENO^NCY. 


100 


12E 


130 


135 


140 


145 


DATA  : 

X1*X2 

*X3 

*X4/.l 

*1. 

* . 1 * 

IF 

(F 

.LT. 

1. 

29) 

GO 

TO 

160 

IF 

(F 

.LT. 

2. 

8) 

GO 

TO 

100 

IF 

(F 

.LT. 

8. 

0) 

GO 

TO 

125 

IF 

(F 

.LT. 

9. 

35) 

GO 

TO 

130 

IF 

(F 

.LT. 

15 

.5) 

GO 

TO 

135 

IF 

(F 

.LT. 

35 

.0) 

GO 

TO 

140 

IF 

(F 

.LT. 

70 

.0) 

GO 

TO 

145 

IF 

(F 

.LT. 

94 

.0) 

GO 

TO 

150 

WRITE 

(6*1502) 

Y1 

s 

1.5E-4 

Y2 

s 

1.8E-3 

Y3 

s 

9.1E-4 

Y4 

s 

1.05E-2 

FI 

= 

1.29 

F2 

r 

?.0 

GO 

TC 

1 155 

ri 

s 

9.1E-4 

Y2 

= 

1.05E-2 

Y3 

s 

1.3E-2 

Y4 

s 

0.18 

FI 

s 

2.8 

F2 

r 

8.0 

GO 

TO 

155 

Y1 

= 

1.3E-2 

Y2 

s 

0.18 

Y3 

r 

2.0E-2 

Y4 

s 

0.32 

FI 

s 

8.0 

F2 

= 

9.35 

GO 

TC 

1 155 

Y1 

= 

2.0E-2 

Y? 

s 

0.32 

Y3 

s 

6.1E-2 

Y4 

= 

1.3 

FI 

s 

9.35 

F? 

r 

15.5 

GO 

TO 

i 155 

Y1 

= 

6.1E-2 

Y? 

s 

1.3 

Y3 

r 

0.41 

Y4 

= 

5.8 

FI 

r 

15.5 

F2 

s 

35.0 

GO 

TO 

155 

Y1 

s 

0.41 

Y2 

3 

5.8 

Y3 

S 

1.5 
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V4  « 10.1 
FI  « 35.0 
F2  « 70.0 
GO  TO  155 
150  Y1  « 1.50 

T2  s 10.1 

Y3  « 1.40 

Y4  » 12.0 

FI  s 70.0 

F2  » 94.0 

155  31  s »LOG  «Y2  / Yll  / •LOG  (X2  / Xl» 

•1  = Y1  / (XI  • * Bll 

32  « •LOG  (Y4  / Y3»  / •LOG  (X4  / X3) 

•2  = Y3  / (X3  • • B2» 

• s (((F  • F - FI  • Fl>  * (^2  - «1)>  / (F2  • F2  - FI  • FD)  ♦ •! 

3 = ((F2  - F»  / (F2  - Fin  • B1  ♦ ((F  - FI)  f (F2  - FI)  ) • B2 

GO  TO  1B5 

160  • » 0.000973  • F • F 

9 3 1.0 

165  IF  (F  .LT.  1.29)  GO  TO  210 
IF  (F  .LT.  2.8)  GO  TO  170 
IF  CF  .LT.  8.0)  GO  TO  175 
IF  (F  .LT.  9.35)  GO  TO  180 

IF  (F  .LT.  15.5)  60  TO  185 

IF  (F  .LT.  35.0)  GO  TO  190 

IF  (F  .LT.  70.0)  GO  TO  195  . . ' 

IF  (F  .LT.  94.0)  GO  TO  20n  ^ .v- 

170  Y1  = 530.0 
Y2  = 2.1E+4 
Y3  = 600.0 
Y4  = 2.35*4 
FI  = 1.29 
F2  = 2.80 
60  TO  205 
IT*^  Y1  s 60C.C 
Y2  = 2.3E*4 
Y3  = 690.0 
Y4  = 2.55*4 
FI  = 2.8 
F2  = 8.0 
GO  TO  205 
180  Y1  = 690.0 
Y2  « 2.55*4 
Y3  = 61C.0 
Y4  = 2.15*4 
FI  = 8.0 
F?  s 9.35 
GO  TO  205 
185  Y1  = 610.0 
Y2  = 2.15*4 
Y3  = lOCO.O 
Y4  = 3.35*4 
FI  = 9.35 
F2  = 15.5 
GO  TO  205 
193  Y1  = 1000.0 
Y2  = 3.35*4 
Y3  = 890.0 
Y4  = 1.25*4 
FI  = 15.5 
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F2  = 35.0 
GO  TO  205 
195  Y1  » 890.0 
Y2  = i.eE*** 

Y3  = 170.0 
Y<i  > 910.0 
FI  » 35.0 
F2  » 70.0 
GO  TO  205 
200  Y1  » 170.0 
Y2  s 910.0 
Y3  s 51.0 
Y4  3 260.0 
FI  a 70.0 
F2  * 94.0 

205  01  = ALOG  <Y2  / Tl)  / ALOG  <X2  / XI) 

Cl  » Y1  / (XI  • • 01) 

02  s ALOG  (Y4  t Y3)  / ALOG  (X4  / X3) 

C2  a Y3  / (X3  * ♦ 02) 

C _ • F - FI  * FI)  • (C2  - CD)  / (F2  • F2  - FI  • FD)  ♦ Cl 

0 = l(F2  - F)  / (F2  - FD)  • 01  ♦ ((F  - FD  / (F2  - FD)  • 02 
CO  TO  220 
210  WRITE  16,1500) 

1500  F0RMAT(48H  FREQUENCY  TOO  LON  FCR  REFLECTIVITY  IN  SUB  CRANE) 

1502  FORHAT(32M  FREQUENCY  TOO  HIGH  IN  SUE  CRANE) 

220  RETURN 
ENO 


FUNCTIOK  TERP  (N,  P,  X,  Y) 

THIS  FUNCTION  DOES  LINEAR  INTERPOLATION  ANO  EXTRAPOLATION. 
NaNUMBER  OF  DATA  POINTS  IX, Y)  TO  BE  USED  IN  INTERPOLATION. 
P=  X VALUE  THAT  PRODUCES  INTERPOLATED  Y VALUE. 

OIHENSION  X(75),Y(75) 

IF  (N  .LT.  2)  GO  TO  240 

DO  225  I a 2,  N 

IF  (P  - X (I))230,  2T5,  225 


225 

240 

CONTINUE 

CONTINUE 

I a H 

230 

TFRP  a V 
1-  X (I  - 
RETURN 

(I  - 1)  ♦ (Y  (I)  -Yd-  D)  • (P  - X (I  - D)  / (X  (I) 

D) 

235 

TPRP  a y 

RETURN 

ENO 

(I) 

• 
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SUBROUTINF  TOPOXY IT .P,PLAHO*«L«GAMAUt ARRON.GAMAUP) 

THIS  SUBROUTINF  CALCULATFS  ATHOSPHFRIC  ABSORPTION  PER  UMT  LENGTH 
DUE  TO  ATMOSPHERIC  OXYGEN  FOP  FREQUENCIES  ROUGHLY  IN  THE  RANGE 
4S<>70  GHI.  SEE,  AS  AN  EXAMPLE,  FALCONE,  V.J*,  ATMOSPHERIC 
ATTENUATION  OF  HICRONAVE  PONER,  J.  MICRONAVE  POMER  ICANAOAI,  VOL. 5 
NO.  4,  DEC.  1970,  PP.  269-278, 

INPUT 

t*tempera7ure  in  OEGREES  KFLVIN. 

PsPRESSURE  IN  MILLIBARS. 

PLAHOAsHAVELENGTH  IN  CENTIMETERS. 

L*ACCOUNTEO  for  by  TEST30. 

ARRON*ATMOSPHERIC  HATER  VAPOR  DENSITY  ON  TRANSMISSION  PATH  OF 
INTEREST. 

OUTPUT 

GAMAUsOXYGEN  ABSORPTION  COEFFICIENT  IN  OB/KM. 

GAMAUP*PHASE  DISPERSION  OUE  TO  OXYGEN  IN  RAOIANS/KM. 


■ 


1 

t 


COMMON  /9L0CX2  /PMUPL  (491,  PMUH  (49),  PHUNOT  (49),  RSRLNl  (49),  R 
1SRLN2  (49) 

REAL  L3 
C 

VP  = T ♦ ARROM  / 216.68 
GOTOdOO.lOE)  ,L 
IOC  XI  = .021333 
X2  s .04523 
X3  = .36748 
X4  » .027351 
X5  a 1. 

L3  a ALOG  (X3) 

OXOLOG  = (X2  - XI)  / (ALOG  (X4)  - L3) 

105  RLMOA  a 1.0  / PLAHDA 
B9  a 2,068666098  / T 
SUM  a o.n 
SUMP  a 0.0 
FEE  a 0.0 

NN  a 1,9 

X a (P  ♦ VP)  / 1013,25 
IF  (X  .GT.  X3)  GO  TO  110 

IF  (X  .LT,  XL)  GO  TO  11' 

OLT2aX140XOLOG*(ALOG(X)-L3) 

OLTlaOLTZ 
GO  TO  120 
110  OLTBaXl 
OLTlaXl 
GO  TO  120 
115  0LTlaX2 
0LT2aX2 
GO  TO  120 

12C  niTl  a OLTl  • ((300.  / T)  • • X5>  • X 

0LT2  a CLT2  • ((300.  / T)  • » X')  • X 

9LTNUA  a 0,5  • (CLTl  ♦ rLT2) 

G0T0(125,135) ,L 
125  IF  (NN  .LT.  1)  GO  TO  145 
90  130  K a 1,  nn,  2 

OK  a K 

“NUPL  (K)  a PK  a j?,  • PK  ♦ 3.)  / (PK  ♦ 1.) 
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PHOH  (K)  ■ tPK  ♦ t*)  • <2.  • PK  - 1.)  / PK 

PHUNOT  IKI  » IPK  • PK  ♦ PK  ♦ 1.0)  / <PK  • <PK  ♦ 1.01) 

ISO  PMONOT  (K)  » PHUNOT  (K)  • «2.  • PK  ♦ 1.)  • 2. 

I<i5  CONTINUE 
CALL  FREQ 

135  IF  (NN  .LT.  1)  GO  TO  150 
00  lAO  K ■ 1.  NN.  ? 

PK  « K 

FAC  ■ EXP  I - BB  • PK  • (PK  ♦ 1.)) 

CALL  FARM  (OLTI.  RSRLNl  (K) « RLMOA.  AA.  AAP) 

CALL  FARM  (OLT2.  RSRLN2  (Kl * RLMOA*  AB.  A8P) 

CALL  FARM  (OLTNUA.  FFE*  RLMOA.  AC*  ACP) 

AC  « AC  • 0.5 

TERM  * (AA  • PMUPL  (K)  ♦ AB  • FMUM  (K)  ♦ AC  • PMUNOT  (K))  • FAC 
TERMP  « (AAP  • PMUPL  (K)  ♦ ABP  • PMUM  (K)  ♦ ACP  • PMUNOT  (K))  • FA 
1C 

SUM  s SUM  « TERM 
SUMP  « SUMP  ♦ TERMP 
1*»0  CONTINUE 
150  CONTINUE 

GAMAU  * SUN  • 59.4681826  • P / (T  ♦ T • T) 

GAMAU  > GAMAU  / (PLAMOA  • PLANDA) 

GAMAUP  « SUMP  • 6.846527564  • F / (T  • T • T) 

GAMAUP  * GAMAUP  / (PLAMOA  • PLAMOA) 

RETURN 

END 


SUBROUTINE  FREQ 

THIS  SUBROUTINE  IS  ASSOCIATED  KITH  SUBROUTINE  TOPOXV,  BUT  IS  NOT 
DIRECTLY  ASSOCIATED  KITH  ANY  OF  THE  PARAHETERS  USED  IN  PREOIC. 


COMMON  /8LOCK2  /PMUPL  (49).  FMUM  (49).  PMUNOT  (49).  RSRLMl  (49).  R 
1SRLM2  (49) 

PMU  « - 252.72 
100  DO  125  K « 1*  49.  2 
PK  » K 

IF  (K  - 1)105.  110.  115 
105  STOP 
110  FUDGE  * - 1.0 
GO  TO  120 
115  fudge  « 1.0 
120  CONTINUE 

PLAMOA  * 59501.6  ♦ .057*  • PK  • (PK  ♦ 1.0) 

9 » 43101.6  - .14  • PK  • (PK  ♦ 1.0) 

RSRLMl  (K)  « RSLMOl  (PK*  PLAMOA.  B*  PMU) 

RSRLH2  (K)  « RSLND2  (PK*  PLAMOA*  8*  PMU*  FUDGE) 

125  CONTINUE 
RETURN 
END 
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FUNCTION  RSLNOl  (PKt  PLAHO««  B«  PHU) 


THIS  SUBPOUTINF  IS  ASSOCIATED  NITH  SUBROUTINE  TOPOXY,  BUT  IS  NOT 
DIRECTLY  associated  NITH  ANY  OF  THE  PARAMETERS  USED  IN  PREOIC. 

Y » B - PHU  / 2. 

Y1  » IPX  ♦ PK  ♦ 3.»  • Y 

X * PLAHOA  - PHU  • IPX  ♦ 1.)  - Y1  ♦ SORT  (PLAHDA  • PLAMOA  - (PLANO 
lA  ♦ PLAMOA)  • Y ♦ Y1  • Yl) 

RSLMOl  s X / (2.9B793E»A) 

RETURN 

ENO 


FUNCTION  RSLM02  IPX*  PLAMOA,  B,  PHU,  FUDGE) 

THIS  SUBROUTINE  IS  ASSOCIATED  NITH  SUBROUTIIC  TOPOXY,  BUT  IS  NOT 
directly  ASSOCIATED  NITH  ANY  OF  THE  PARAMETERS  USED  IN  PREOIC, 


Y * B - PMU  / 2. 

Yl  s IPX  ♦ PX  - 1.)  • Y 

X « PLAMOA  ♦ PMU  • PX  ♦ Yl  - FUDGE  • SORT  (PLAMOA  • PLAMOA 
lOA  ♦ PLAMOA)  • Y ♦ Yl  • Yl) 

RSLH02  » X / l2.997q3E»u) 

RETURN 

ENO 
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subroutine  farm  (oltn,  resnav,  actnav,  forma,  FORMP) 


this  subroutine  is  associated  NITH  SUBROUTINE  TOPOXY,  BUT  IS  NOT 
DIRECTLY  ASSOCIATED  NITH  ANY  OF  THE  PARAMETERS  USED  IN  PREOIC, 

X « RESNAV  - actnav 

X2  » X • X 

Y « RESNAV  ♦ actnav 

Y2  * Y * Y 

02  » oltn  • OLTN 

FORMl  * X?  ♦ 02 

FORMl  * 1,  / FORMl 

F0RM2  » Y2  ♦ 02 

F0RM2  ■ 1,  / F0RM2 

FORMA  * OLTN  • (FORMl  ♦ F0PM2) 

F0RMP*(02/ACTNAV4IRESNAV/ACTNAV)*X)*F0RM1*(02/ACTNAV*(RESNAV/ACTNA 

1V)*Y)*F0RM2 

RETURN 

ENO 
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SUBROUTINr  PROMO  (MAVE,0,N) 

. . .THIS  SUBROUTINE  PERFORMS  THE  MODIFICATION  OF  ATTENUATION  OF  NICRO- 
MAVE  TERRESTRIAL  LINKS  FOR  METHODS  1 AND  2.  DESCRIBED  IN  THE  MAIN 
OROGRAH. 

INPUT 

HAVE  - CARRIER  WAVELENGTH  IN  CM. 

0 - PATH  LENGTH  IN  KM. 

N - INTEGER  CORRESPONDING  TO  METHOD  1 OR  2. 

DIMENSION  M(6)«CA(6)  .BS(6)  tRELISdSI  . RELIT  (15>  .RELY  (15)  .ETA  (IE) 
DIMENSION  EMTAU(IE) 

COMHON/FROM/TAUDBT  (IS)  .REVTAUdSi  ,RELI  (15)  .HT0P(12) 

Fs2q.97B3/MAVE 
RFSQ=F*F/225.0 
00  10  1=1.12 
RELY(I)=RELI(I) 

IF(N.EO.l)  EMTAU(I)=(RFSQ*0.987/TAUOBT(I))*(SQRT(16994.70663*HTOP( 
lI))/0) 

IF(N.EQ.2)  EMTAU(I)=(RFSQ*0.987/TAUOBT(I))*(SQRT(1699i».70663*HTOP( 
II) )/D)«2.0 

RELIT(I)=RFLY(I)*EMTAU(I) 

10  CONTINUE 

DO  20  1=1,12 

IF(RELIT(I).GE.RELY(I)I  6,8 
6 EMTAU(I)=1.0 

REWTAU(1)=TAUDBT (I) 

GO  TO  20 

8 REVTAU(I)=EXTERP(12,RELY(I) .RELIT .TAUOBT) 

20  CONTINUE 
9ETURN 
ENO 
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FUNCTION  EXTERP(N»PfX,Y» 


C 

c this  function  pfrfopns  fxponential  interpolation  and  extrapolation 
C ON  Y VALUFSt  GIVEN  X VALUES. 

C 

C INPUT 

C 

C N - NUHBER  OF  DATA  POINTS. 

C P - POINT  ON  X-COOROINATE  NEEOING  A CORRESPONDING  Y INTERPOLATION. 

C X - ARRAY  OF  N X-COOROINATE  DATA  POINTS. 

C Y - ARRAY  OF  N Y-COOROINATE  OATA  POINTS. 

C 

OINENSION  XIEOl.YISOI 

00  1 I>2.N 

IF  (X(I»-P)  1,?,2 

1 CONTINUE 

1 « N 

2 IF  IYai*Y(I-in  3,3,4 

3 cxTERP  = (Y(I>-Y(I-1))*<P-X(I-1))/(X(I)-X(I-1))  ♦ Y(I-l) 

RETURN 

4 EXTERP  = Y(I-1»*EXP  ( ALOG(Y (I ) /Y ( I-l ) » • (P-X (I-l ) ) / (X ( I ) -X (I-l)  ) I 
RETURN 

ENC 


ooooonooooooo(^nnoonoooonr>r>nooooor>noooooooonoonooonooooooooo 


SUBR0U7INF  TABLES«XLAT«XLON,ELEV,PREStTMP,RLHUM,BET,H) 


j 


TABLES  FINOS  THE  CLOSEST  20  OF  THE  Z'jB  SAMPLE  OATA  STATIONS  | 

(OUTTON  ET  AL.t  197A)  TO  AN  INTENOEC  INTERPOLATION  POINT.  TABLES 
THEN  CALLS  lOBVIP  (AKINA,1975)  TO  OBTAIN  METEOROLOGICAL  CONOITIONS 
AT  THE  INTERPOLATION  POINT.  FOR  BEST  RESULTS,  THE  INTERPOLATION 
POINT  SHOULO  NOT  BE  OUTSIDE  OF  THESE  LATITUDE  AND  LONGITUDE 
BOUNDS  100. MM)  — 92.00N  TO  ItQ.OQN  BY  1.20E  TO  25.00E  ••  WHICH 
CORRESPOND  ROUGHLY  TO  THE  LIMITS  OF  THE  ENCS  ZONE.  IF  THE  INTER- 
POLATION POINT  IS  OUTSIDE  THESE  BOUNDS,  A WARNING  WILL  BE  PRINTED 
BUT  PROCESSING  WILL  CONTINUE. 


INPUT  

ELEV  - ELEVATION  OF  THE  INTERPOLATION  POINT  IN  METERS. 

XLAT  - DEGREE-MINUTES  (00. MM)  LATITUDE  OF  INTERPOLATION  POINT. 
XLON  - OEGREE-MINUTES  (DO. MM)  LONGITUDE  OF  INTERPOLATION  POINT. 


OUTPUT  

SET  - INTERPOLATED  VALUE  OF  BETA  (DIMENSIONLESS). 

M - INTERPOLATED  TOTAL  ANNUAL  PRECIPITATION  (MILLIMETERS). 
PRES  - INTERPOLATED  SURFACE  PRESSURE  (MILLIBARS). 

RLMUM  - INTERPOLATED  RELATIVE  HUMIDITY  (DECIMAL  FRACTION). 
TMP  - interpolated  TEMPFRATURE  (DEGREES  CENTIGRADE). 


DESCRIPTION  OF  MAJOR  VARIABLES  

BETA  - ARRAY  CONTAINING  THE  BETAS  OF  THE  20  CLOSEST  STATIONS. 
DATAPT  - ARRAY  CONTAINING  METEOROLOGICAL  OATA  FOR  ALL  2L9  STATIONS 
ARRAINGED  3Y  COLUMNS  - 1=PRESSURE,  2=TEMPERATURE , 

RELATIVE  HUMIDITY,  LsBETA,  5*1. 

OCXLAT  - DECIMAL  DEGREES  OF  XLAT. 

OCXLOK  - DECIMAL  DEGREES  OF  XLCN, 

DEC-'.JIT  - ARRAY  CONTIANING  DECIMAL  LATITUDES  OF  21.9  DATA  STAIONS. 
OECLCM  - ARRAY  CONTAINING  DECIMAL  LCNGITUOES  OF  249  OATA  STATIONS. 
ELEV’  - ARRAY  CONTAINING  ELEVATIONS  OF  ALL  249  STATIONS. 

ELV  - VALUE  OF  ELEV  ROUNDED  TO  THE  NEAREST  CONTOUR  INTERVAL. 

ISU8  - ARRAY  CONTAINING  THE  ARRAY  SUBSCRIPTS  OF  THE  20  CLOSEST 
STATIONS. 

TWK  - ARRAY  INTERNAL  TO  lOBVIP  WHICH  MUST  BE  DIMENSIONED  IN  THE 
CALLING  PROGRAM  (AKIMA,  1975). 

LAT  - ARRAY  CONTAINING  THE  OEGREE-MINUTES  (00. MM)  LATTITUDE  OF  249 
STATIONS. 

LON  - ARRAY  CONTAINING  THE  OEGREE-MINUTES  (00. MM)  LONGITUDE  OF  249 
STATIONS. 

LATT  - ARRAY  CONTAINING  THE  DECIMAL  LATITUDES  OF  THE  20  CLOSEST 
STATIONS. 

LONG  - ARRAY  CONTAINING  THP  DECIMAL  LONGITUDES  OF  THE  20  CLOSEST 
STATIONS. 

MM  - ARRAY  containing  H VALUfS  OF  THE  20  CLOSEST  STATIONS. 

N - NUMBER  OF  STATIONS. 

NCP  - NUMBER  OF  ADDITIONAL  POINTS  FCR  IDBVIP  TO  USE  IN  CALCULATION 
(INTERNAL  TO  lOBVIP-POlNTS  NOT  SUPPLIED  SEPARATELY  BY  TABLES 
) . 

NOP  - NUMBER  OF  CLOSEST  STATIONS  USED  IN  THE  INTERPOLATION  OF 


1 
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C PRESSURE,  TEHPERRTURE,  RELATIVE  HUMIDITY  AND  BETA. 

C NORM  - NUMBER  OE  CLOSEST  STATIONS  USED  IN  INTERPOLATING  H VALUES. 

C NIP  - NUMBER  OF  POINTS  TO  INTERPOLATE  AT, 

C PRESS  - ARRAY  CONTAINING  PRESSURES  OF  THE  20  CLOSEST  STATIONS. 

C RELHUH  - ARRAY  CONTAING  REL.  HUMIDITY  OF  THE  20  CLOSEST  STATIONS. 

C TEMP  - array  CONTAINING  THE  TEMPERATURE  OF  THE  20  CLOSEST  STATIONS 
C NR  - ARRAY  INTERNAL  TO  lOBVIP  WHICH  MUST  BE  DIMENSIONED  IN  THE 
C CALLING  PROGRAM  (AKIMA,  1975). 

C 

c 

C 

DIMENSION  ELEV1I2<»9),PRESSI20) ,TEMP(20) ,RELHUM(20) ,OETA(20) , 
10ECLAT(249) tDECL0N(249) ,HK(160> ,OATAPT ( 249 , 5) , ISUB ( 20 ) ,PRS (249 > , 
2IHK(641) 


REAL  LAT(249),LON(249),MM(20),LATT(20).LONG(20),M 
C 

C COMMON  REQUIRED  BY  lOBVIP  (AKIMA,  1975). 

C 

t COMMON/IOCM/NCP 

C 

C DATA  STATEMENTS  CONTAINING  POSITION, ELEVATION  AND  METEOROLOGICAL 

I C DATA  ON  249  STATIONS. 


OATA  LAT(  l),LON(  D.ELEVK  1)/  68.98,  33.03,  46./ 

DATA  OATAPK  1 , 1 ) , OATAPT ( l,2),OATAPT(  1,3),0ATAPT(  1,4), 

lOATAPK  1,5)/1004.4,  .2,  .78,  .039,  384.0/ 

OATA  LAT(  2),L0N(  2),ELEV1(  2)/  68.39,  43.18,  0./ 

OATA  OATAPT(  2,l),OATAPT<  2,2),OATAPT(  2,3),0ATAPT(  2,4), 

10ATAPT(  2,9)/1011.2,  -1.9,  .81,  .036,  329.0/ 

DATA  LAT(  3),LON(  3),ELEV1(  3)/  70.22,  31.06,  15./ 

OATA  OATAPT ( 3,1),0ATAPT(  3,2),0ATAPT(  3,3),OATAPT(  3,4), 

inATAPT(  3,9)/1007.4,  1.6,  .89,  .021,  949.0/ 

OATA  LAT(  4),L0N(  4),ELEV1(  4)/  60.08,  -1.11,  82./ 

OATA  OATaPT(  4,1),DATAPT(  4,2),0aTAPT(  4,3),DATAPT(  4,4), 

10ATAPT(  4,5)/1000.0,  7.2,  .93,  .034,1129.0/ 

OATA  LAT(  9),LON(  9),ELEV1(  9)/  98.13,  -6.20,  3./ 

DATA  OATAPT(  9,l),OATAPT(  5,2),DATAPT(  9,3),OATAPT(  9,4), 

10ATAPT(  9,9)/1010.0,  8.6,  .89,  .091,1041.0/ 

OATA  LAT(  6),L0N(  6),ELEV1(  6)/  96.30,  -6.93,  10./ 

OATA  OATAPT(  6,1),0ATAPT(  6,2),0ATAPT(  6,3),OATAPT(  6,4), 

10ATAPT(  6,9)/1010.0,  9.8,  .84,  .030,1128.0/ 

OATA  LAT(  7),L0N(  7),fLrvi(  7)/  97.12,  -2.12,  59./ 

OATA  OATAPT(  7,1), OATAPK  7,2),DATAPT(  7,3),OATAPT(  7,4), 

10ATAPT(  7,5)/1004.6,  7.9,  .87,  .038,  837.0/ 

OATA  LAT(  8),LON(  8),ELEV1(  8)/  55.97,  -3.21,  242./ 

OATA  OATAPT(  8,l),OATAPT(  e,2),DATAPT(  8,3),OATAPT(  8,4), 

inATAPT(  8,9)/1008.2,  7.7,  .88,  .040,  696.0/ 

DATA  LAT(  9),L0N(  9),ELEV1(  9)/  59.19,  -3.12,  32./ 

OATA  0ATAPT(  9,1),DATAPT(  9,2),DATAPT(  9,3),0ATAPT(  9,4), 

10ATAPT(  9,9)/  983.7,  7.1,  .89,  .049,1927.0/ 

OATA  LAT(  10),LON(  10),ELEV1(  10)/  94.18,  -1.32,  68./ 

OATA  OATAPT(  10 , 1) , 0 A T APT ( 10 ,2) , OATAPT ( 10 ,3) , OATAPT ( 10,4), 
10ATAPT(  10,5)/1009.9,  8.0,  .89,  .070,  649.0/ 

OATA  LAT(  11),L0N(  11),FLEV1(  11)/  94.39,  -6.13,  39./ 

OATA  OATAPT(  11 , 1 ) , D A TAPT ( 1 1 , 2) , OATAPT ( 1 1 , 3)  .OATAPT ( 11,4), 


[ 
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DATA  LATI  13).L0N(  13),ELEV1(  13)/  53.31.  *2.16.  75./ 

OATA  OATAPT(  13 « 1 ) « 0 ATAPTt  13 , 2) , 0 ATAPT ( 13 .3) »OAT APT ( 
lOATAPTI  13.5)/100<».8,  9.5.  .85.  .079.  799.0/ 

OATA  LATI  ID.LONI  14).CLeVl(  141/  52.27.  >1.45.  97./ 

OATA  OATAPTI  14.1) .OATAPT I 14.2) .OATAPT  I 14.3) .OATAPT  I 
lOATAPTI  14.5)/1003.0.  9.7.  .84.  .070.  727.0/ 

OATA  LATI  15).L0NI  15),ELrVll  15)/  53.10.  -.31.  68./ 

OATA  OATAPTI  15 . 1 1 , OA TAPT I 15 .2) .OATAPT I 15 .3) . OAT APT  I 
lOATAPTI  15.5)/1006. 0.  9.4.  .87.  .045.  597.0/ 

OATA  LATI  16).L0N{  16).ELEV1{  16)/  52.35.  1.43«  4./ 

OATA  OATAPTI  16 . 1) . OAT APT  I 16 . 2) . OAT APT  I 16 . 3) . OAT APT  I 
lOATAPTI  16.5)/1014. 3.  9.9.  .84.  .060.  583.0/ 

OATA  LATI  17),L0NI  17).ELEW1I  17)/  51.28,  -.19.  5./ 

OATA  OATAPTI  17 , 1 ) . OATAPT I 17 ,2) . OATAPT I 17,3) .OATAPT I 
lOATAPTI  17.5)/1014.7,  10.5.  .81.  .069.  594.0/ 

OATA  LATI  18).LONI  ID.ELEVK  18)/  51.09,  -.11.  60./ 

OATA  OATAPTI  18 , 1 ) . 0 AT  APT  I 18,2) , OATAPT  I 18,3) .OATAPT  I 
lOATAPTI  18.5)/1008.1.  10.3,  .81.  .070,  683.0/ 

OATA  LATI  19),L0NI  19),ELEV1I  19)/  51.24,  -3.21,  67./ 

OATA  OATAPTI  19 . 1 ) . D ATAPT I 19 .2 ). OATAPT I 19 .3) . OAT APT  I 
lOATAPTI  19.5)/1006.7,  10.2.  .81,  .060,1100.0/ 

OATA  LATI  20),LONI  20),ELEV1I  20)/  50.47,  -1.50,  11./ 

OATA  OATAPTI  20 , 1 ) , 0 ATAPT I 20 , 2 ), OAT APT  I 20 ,3) .OATAPT  I 
lOATAPTI  20.5)/1014, 3,  9.9,  .81,  .070,  800.0/ 

OATA  LATI  21).L0NI  21).ELEV1I  21)/  50.21,  -4. 07,  27./ 

OATA  OATAPTI  21 , 1 ) , OATAPT I 21 ,2 ), OATAPT I 21 ,3) , OATAPT  ( 
lOATAPTI  21,5)/1012.2.  10.8,  .87.  .068,  990.0/ 

OATA  LATI  22),LONI  22),ELEV1(  22)/  55.22,  -7.20,  25./ 

OATA  OATAPTI  22 . 1 ) , 0 AT  APT  I 22 . 2 ), OAT APT  I 22 . 3) , OAT  APT  I 
lOATAPTI  22.5)/1008.9,  9.5,  .83.  .025,  980.0/ 

OATA  LATI  23),L0NI  23),FLEVH  23)/  54.14,-10.0  0,  10./ 

OATA  OATAPTI  23 , 1 ) . 0 AT APT  I 23 ,2) , OATAPT I 2 3 .3) , OAT APT  ( 
lOATAPTI  23.5)/lQ11.2.  10.0,  .84.  .040,1133.0/ 

OATA  LATI  24).L0NI  24).ELEV1I  24)/  53.26,  -6.15,  81./ 

OATA  OATAPTI  24 , 1) , OATAPT I 24 , 2) , OAT APT  I 24,3) .OATAPT ( 
lOATAPTI  24.5)/1003.5,  9.6,  .83,  .035.  758.0/ 

OATA  LATI  25),L0NI  25),fLEV1(  25)/  52.41,  -8.55,  7./ 

OATA  OATAPTI  25 , 1 ) , 0 A T APT  I 25 , 2 ), OATAPT  I 25 ,3) .OATAPT  I 
lOATAPTI  2*^,5) /1012. 8,  10.2,  .86,  . 062,  930.0/ 

OATA  LATI  26),L0NI  26).ELEV1I  26)/  51.51,  -8,30,  162./ 
OATA  OATAPTI  26 , 1 ) , OATAPT I 26 , 2 ), DAT APT  I 26 ,3) .OATAPT  ( 
lOATAPTI  26.5)/  994.4,  10.2,  .83,  .055,1049.0/ 

OATA  LATI  27),L0NI  2r),FLFVll  27)/  51.56,-10.15,  14./ 

OATA  OATAPTI  27 . 1 ) , 0 AT APT  I 27 , 2 ), DAT APT  I 27,3) .OATAPT  I 
lOATAPTI  27,5) /1012. 5,  10.8,  .83,  .076,1398.0/ 

OATA  LATI  2e),LONI  28),ELEV1I  28)/  57,06,  9.52,  3./ 

OATA  OATAPTI  28 , 1 ) , 0 A TAPT I 28 ,2 ), OATAPT  I 28 ,3) .OATAPT ( 
lOATAPTI  28,5) /1012. 3,  7.6,  .83,  .036,  576,0/ 

OATA  LATI  29),L0NI  29),ELEV1I  29)/  55.41,  12.33,  *=./ 

OATA  OATAPTI  29, 1 ), OAT APT  I 29,2) .OATAPT  I 29 , 3) , OAT APT  I 
lOATAPTI  29,5) /1013. 1,  8,5,  .79,  .048,  602.0/ 

OATA  LATI  30),LONI  30),ELEV1I  30)/  55.00,  15.05,  6./ 

DATA  OATAPTI  30 , 1 ) , OA TAPT I 3C ,2) .OATAPT  I 30 ,3) .OATAPT  I 
lOATAPTI  30,5)/1013,4,  8.0,  .84,  ,033,  553,0/ 

OATA  LATI  31),LONI  31),ELEV1(  31)/  63.42,  9.37,  7./ 

OATA  OATAPTI  31 , 1) , 0 A T APT  I 31 , 2) , OATAPT  I 31 , 3) ,OAT APT  I 
lOATAPTI  31,5) /1009. 3,  5.9,  .82,  .055,  999.0/ 

OATA  LATI  32),LONI  32),ELEV1I  32)/  63.25,  10.26,  115./ 
OATA  OATAPTI  32 , 1 ) , OA TAPT I 32, 2 ), OATAPT I 32,3) .OATAPT I 
lOATAPTI  32,5)/  996.9,  4.9,  .78,  .045,  857.0/ 


13.4)  , 

14.4) , 

15.4) , 

16.4)  , 

17.4) , 

18.4) , 

19.4)  , 

20.4) , 

21.4)  , 

22.4) , 

23.4)  , 

24.4)  , 

25.4) , 

26.4)  , 

27.4) , 

28.4)  , 

29.4) , 

30.4)  , 

31.4)  , 

32.4)  , 
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OAT«  LATt  331 
DATA  OATAPTI 
lOATAPU  33t5» 
DATA  LATI  3AI 
OATA  OATAPH 
l']ATAPT(  34.5) 
OATA  LATI  35) 
OATA  OATAPTI 
lOATAPTi  35t5> 
OATA  LATI  36) 
OATA  OATAPTI 
lOATAPTI  36«5I 
OATA  LATI  37) 
OATA  OATAPTI 
lOATAPTi  37,5) 
OATA  LATI  36) 
OATA  OATAPTI 
lOATAPTi  38,5) 
OATA  LATI  39) 
OATA  OATAPTI 
lOATAPTi  39,5) 
OATA  LATI  90) 
OATA  OATAPTI 
lOATAPTi  40,5) 
OATA  LATI  41) 
OATA  OATAPTI 
lOATAPTi  41,5) 
OATA  LATI  47) 
OATA  OATAPTI 
lOATAPTi  47,5) 
OATA  LATI  43) 
OATA  OATAPTI 
lOATAPTi  43,5) 
OATA  LATI  44) 
OATA  OATAPTI 
lOATAPTi  44,5) 
OATA  LATI  45) 
OATA  OATAPTI 
lOATAPTi  45,5) 
OATA  LATI  46) 
OATA  OATAPTI 
lOATAPTi  46,5) 
OATA  LATI  47) 
OATA  OATAPTI 
lOATAPTi  47,5) 
OATA  LATI  48) 
OATA  OATAPTI 
lOATAPTi  48,5) 
OATA  LATI  49) 
OATA  OATAPTI 
lOATAPTi  49,5) 
OATA  LATI  50) 
OATA  OATAPTI 
lOATAPTi  50,5) 
OATA  LATI  51) 
OATA  OATAPTI 
lOATAPTi  51,5) 
OATA  LATI  57) 
OATA  OATAPTI 
lOATAPTi  57,5) 


,LONI 

33.1) 
/1006 
,LONI 

34.1) 
/lOll 
,LONI 

35.1) 
/ 999 
,LONI 

36.1) 
/1003 
,LONI 

37.1) 
/1014 
,LONI 

38.1) 
/1008 
,LONI 

39.1) 
/ 990 
,LONI 

40.1) 
/ 998 
,LONI 

41.1) 
/ 990 
,LONI 

47.1) 
/ 996 
,LONI 
43«1) 
/ 983 
,LONI 

44.1) 
/lOll 
,LONI 

45.1) 
/ 999 
,LONI 

46.1) 
/1003 
,LONI 

47.1) 
/1015 
,LONI 

48.1) 
/ 975 
,LONI 

49.1) 
71007 
,LONI 

50.1) 
71017 
,LONI 

51.1) 
71014 
,LONI 

57.1) 
71P05 


33) ,CLrVl{  33) 
OATAPTI  33,7), 

4,  7 ,8,  , 87 , 

34) «rLEVll  34) 

OATAPTI  34,7), 
1 , 7,4 f , 83 , 

35) ,rLEVll  35) 
OATAPTI  35,7), 
6 , 11*3,  ,87 , 

36) ,ELEV1I  36) 
OATAPTI  36,7), 
8,  10,8,  , 86 , 

37) ,ELEV1I  37) 
OATAPTI  37,7), 
1,  11.7,  .53, 
3S),ELEV1I  38) 
OATAPTI  38,7), 

6,  10.9,  .78, 
39I,ELEV1I  39) 
OATAPTI  39,7), 

5,  9.5,  .87, 

40) ,ELEV1I  40) 

OATAPTI  40,7), 
7 , 9.7 , . 5C  , 

41) ,ELEV1I  41) 
OATAPTI  41,7), 
1,  10.5,  . 77 , 
4?),FLEV1I  47) 
OATAPTI  47,7), 

7,  11.1,  .78, 

43) ,ELEV1I  43) 
OATAPTI  43,7), 
3,  10.6,  .78, 

44) ,ELEV1I  44) 
OATAPTI  44,7), 
7,  17.3,  .51, 

45) ,FLEV1I  45) 
OATAPTI  45,7) , 
1,  17.5,  .76, 

46) ,ELEV1I  46) 

OATAPTI  46,7) , 
0,  9.9,  .83, 

47) ,ELEV1I  47) 

OATAPTI  47,7) , 
3 , 9.3,  . 83 , 

48) ,5LEV1I  48) 
OATAPTI  46,7), 

5,  *.8,  .61, 
4')),FLEV1I  49) 
OATAPTI  49,7), 
9 , C . 7,  . 55  , 
51),FLEV1I  50) 
OATAPTI  50,7), 

, 8.5,  . 84  , 

51),ELEW1I  51) 
OATAPTI  51,7), 

6 , 5.8,  . 86  , 
57) ,ELEV1I  57) 
OATAPTI  57,7), 

7 , 5.9,  . 5? , 


/ 60.74,  5.19,  44./ 
OATAPTI  33,3) .OATAPTI 
.134,1958.0/ 

/ 58.53,  5.38,  6./ 
OATAPTI  34,3) .OATAPTI 
.140,1016.0/ 

/ 49.39,  -1.78,  138./ 
OATAPTI  35,3) .OATAPT  I 
.078,  931.0/ 

/ 48.77,  -4.75,  103./ 
OATAPTI  36,3) .OATAPT  I 
.065,1176.0/ 

/ 47.10,  -1.37,  77./ 
OATAPTI  37,3) .OATAPTI 
.068,  757.0/ 

/ 46.55,  7.77,  65./ 
OATAPTI  38,3) .OATAPTI 
.097,  505.0/ 

/ 48.41,  6.13,  717./ 
OATAPTI  39,3) .OATAPT  I 
.091,  717.0/ 

/ 46.33,  7.38,  154./ 
OATAPTI  40,3) .OATAPTI 
.098,  607.0/ 

/ 47.16,  5.05,  777./ 
OATAPTI  41,3) .OATAPTI 
.099,  777.0/ 

/ 47.04,  7.77,  166./ 
OATAPTI  47,3) .OATAPTI 
.089,  671.0/ 

/ 45.49,  1,17,  396./ 
OATAPTI  43,3) , OATAPTI 
.047,  937,0/ 

/ 44.50,  -.47,  51./ 
OATAPT  I 44,3) , OATAPT  I 
.114,  901.0/ 

/ 43.35,  1.77,  153./ 
OATAPTI  45,3) .OATAPTI 
.140,  659.0/ 

/ 50.45,  4.71,  104./ 
OATAPTI  46,3) .OATAPT  I 
.060,  785,0/ 

/ 57.06,  5.11,  -0./ 
OATAPTI  47,3) .OATAPT  I 
.173,  767.0/ 

/ 49.37,  6.03,  378,/ 
OATAPT  I 45,3) .OATAPT  I 
.070,  ^40.0/ 

/ 54.37,  9.33,  40./ 
OATAPTI  49,3) .OATAPTI 
.057,  575,0/ 

/ 53.38,  10.00,  16./ 
OATAPTI  50,3) .OATAPT  I 
.136,  716.0/ 

/ 53.77,  7.13,  17./ 
OATAPTI  51,3) .OATAPT  I 
.097,  767.0/ 

/ 57.77,  9.43,  54,/ 
OATAPTI  57,3) .OATAPTI 
.091,  665.0/ 


33.4)  , 

34.4) , 

35.4) , 

36.4)  , 

37.4)  , 

38.4) , 

39.4) , 

40.4)  , 

41.4)  , 

47.4) , 

43.4)  , 

44.4) , 

45.4)  , 

46.4) , 

47.4)  , 

45.4) , 

49.4)  , 

50.4) , 

51.4) , 

52.4) , 
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OAT*  LAT(  53),L0N(  5T),rLEWl(  53)/  51.19,  9.29,  163./ 

OATA  DATAPTt  53 , 1 ) , OATAPT ( 53 , 2 ) , 0 AT APT ( 53 ,3) , OAT APT ( 53,6), 
lOATAPTt  53,5)/  996.3,  <».0,  .76,  .106,  628.0/ 

OATA  LATI  56),L0N(  54),ELEV1(  56)/  51.26,  6.56,  161./ 

OAT*  DATAPTI  56 , 1) , OATAPT(  56 , 2) , 0 AT APT ( 56 ,3) .OATAPT ( 56,6), 
lOATAPTI  56,5)/  995.7,  9.6,  .8f,  .112,  882.0/ 

OATA  LAT(  55),L0N«  55),ELEV1I  55)/  68.50,  9.12,  315./ 

OATA  DATAPTI  55 , 1 ), OATAPT ( 55 ,2) , OATAPT ( 55 ,3) , OATAPT ( 55,6), 
104T*PT(  55,5)/  979.4,  9.5,  .75,  .113,  687.0/ 

OAT*  L*T<  56),LON(  56),ELPV1(  56)/  56.11,  12.05,  10./ 

OATA  OATAPT!  56 , 1 ) , 0 AT APT ( 56 ,2 ), OATAPT ( 56 , 3) , OATAPT  ( 56,6), 
lOATAPTI  56,5)/1013. 3,  8.5,  .83,  .092,  595.0/ 

OAT*  L*T(  57),L0N(  5T),clpvi(  c7) / 43.2?,  -6.25,  67./ 

OATA  OATAPT!  57 , 1 ) , OAT APT ( 57 , 2 ) , 0 AT APT ( 57 , 3) , OAT APT ! 57,6), 
lOATAPT!  57,5)/1010.1,  13.9,  .79,  .267,  986.0/ 

OATA  LAT!  58), LON!  58),ELEVH  58)/  69.61,  18.55,  10./ 

OAT*  OATAPT!  58 , 1 ) , 0 A TAPT ! 58 , 2 ) , 0 AT APT  I 5 8 , 3) .OATAPT ! 58,6), 
lOATAPT!  58,5)/1007,7,  3.3,  .80,  .065,1019.0/ 

OATA  LAT!  59), LON!  59),ELEVH  59)/  67.17,  16.25,  13./ 

OATA  OATAPT!  59,  1) , OATAPT ! 59 , 2) , DAT APT  ! 59,3) .OATAPT ! 59,6), 
lOATAPT!  59,5) /1007.7,  i*.6,  .78,  . 060,1027.0/ 

OATA  LAT!  60), LON!  60),FLEV1!  60)/  68.27,  22.39,  327./ 

OATA  OATAPT!  60 , 1 ) , OA T APT ! 60 ,2 ), OATAPT ! 60 ,3) .OATAPT  ! 60,6), 
lOATAPT!  60,5)/  970.3,  -1.5,  .90,  .027,  380.0/ 

OATA  LAT!  61), LON!  61),ELEV1!  61)/  65.50,  26.09,  7./ 

OATA  OATAPT!  61 ,1) , OATAPT ! 61 ,2 ), OATAPT  ! 61 ,3) .OATAPT  ! 61,6), 
lOATAPT!  61,5)/1010.2,  -1.6,  .88,  .033,  565.0/ 

OATA  LAT!  62), LON!  62),ELEV1!  62)/  65.36,  17.10,  327./ 

OATA  OATAPT!  62 , 1 ) , 0 A T APT ! 62 ,2 ), OAT APT  ! 62 ,3 ) ,0 AT  APT  ! 62,6), 
lOATAPT!  62,5)/  971.2,  .7,  .88,  .036,  693.0/ 

OATA  LATI  63), LON!  63),CLEV1!  63)/  63.11,  16.30,  366./ 

OAT*  OATAPT!  63 , 1 ) , 0 A TAPT ! 63, 2 ), OAT APT ! 63,3) .OATAPT ! 63,6), 
lOATAPT!  63,5)/  967.3,  2.9,  .86,  .031,  532.0/ 

OAT*  LAT!  66), LON!  66)tELEV1!  66)/  62.39,  17.57,  8./ 

OATA  OATAPT!  66 , 1 ) , OA T APT ! 66 ,2 ), OATAPT ! 66,3) .OATaPT  ! b6<6), 
lOATAPT!  66,5)/1010.7,  6.6,  .82,  .066,  667.0/ 

OATA  LAT!  65), LON!  65),ELEV1!  65)/  67.22,  26.39,  180./ 

OATA  OATAPT!  65 , 1 ) , 0 AT APT ! 65 , 2 ) , 0 ATAPT ! 65 ,3 ) , DAT A PT ! 65,6), 
lOATAPT!  69,5)/  899.1,  -.6,  .80,  .061,  507.0/ 

OATA  LAT!  66), LON!  66),9LEV1!  66)/  66.17,  27.61,  136./ 

OATA  OATAPT!  66 , 1 ) , 0 AT APT ! 66 ,2) .OATAPT ! 66 ,3) .OATAPT I 66,6), 
lOATAPT!  66,5)/  995.0,  1.9,  .75,  .066,  566.0/ 

OATA  LAT!  67), LON!  67),ELFV1!  67)/  62.26,  25.60,  165./ 

OATA  OATAPT!  67, 1) , OATAPT ! 67 , 2 ) , 0 AT  APT ! 67 , 3 ) , OAT APT ! 67,6), 
lOATAPT!  67,5)/  996.0,  2.8,  .82,  .067,  619.0/ 

OATA  LAT!  68), LON!  69),ELFV1!  68)/  63.03,  21.66,  8./ 

OATA  OATAPT!  68 , 1 ) , OA TAPT ! 68 , 2 ) , 0 A T APT ! 68 ,3) .OATAPT  ! 69,6), 
lOATAPT!  68,5)/1010.7,  3.5,  .83,  .066,  682.0/ 

OATA  LAT!  69), LON!  69),FLEV1!  69)/  66.35,  60.30,  13./ 

OATA  OATAPT!  69, 1)  , OATAPT ! 69 , 2 ) , 0 AT  APT ! 69 , 3) , OAT APT ! 69,6), 
lOATAPT!  69,5) /lOlO. 9,  1.6,  .81,  .066,  521.3/ 

OATA  LAT!  70), LON!  70),9LEV1!  70)/  63.69,  30.69,  181./ 

OATA  OATAPT!  70 , 1 ) , 0* T APT ! 7C , 2 ) , D AT APT ! 7 0 ,3) , 0 AT  APT ! 70,6), 
lOATAPT!  70,5)/  989.7,  -1.0,  .79,  .030,  500.0/ 

OATA  LAT!  71), LON!  ’D.^LEVl!  71)/  61.01,  36.27,  59./ 

OATA  OATAPT!  71 , 1 ) , 0 A T APT ! 7 1 , 2 ) , 0 AT APT ! 71 , 3)  , OA T APT ! 71,6), 
lOATAPT!  71,5) /1006, 3,  2.5,  .79,  .072,  'lO.O/ 

OATA  LAT!  72), LON!  7?),pLEV1!  72)/  59.56,  10.66,  96./ 

OATA  OATAPT!  72 , 1 ) , OA TAPT ! 72 ,2 ) , 0 AT APT ! 72 , 3)  , OA T APT ! 72,6), 
lOATAPT!  7?,5)/1000,8,  5.9,  .77,  .115,  760.0/ 
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□»T»  L»T<  73I,L0N(  r3),eLFWl(  73)/  59.31,  17,57,  33./ 

OAT*  OATAPTI  73 , 1 ) , OAT APT  I 73 ,3) ,OATAPT « 73,3) ,OATAPT ( 73, A), 
lOATAPTI  73,5)/10tl.3,  6.6,  .S3,  .046,  545.0/ 

OATA  LAT(  74),L0N(  74),FLEVlt  74)/  59.33,  13.38,  55./ 

OATA  OATAPTI  74, 1 ) , OATAPT C 74,3) ,OATAPT ( 74 ,3) , OATAPT  { 74,4), 
lOATAPTI  74,5)/1006.3,  5.9,  .63,  ,054,  546.0/ 

OAT*  LATI  75),L0NI  75),FLE¥1I  75)/  57.46,  14.05,  333./ 

OAT*  OATAPTI  75,1) .OATAPTI  75 ,3) .OATAPT I 75 ,3) .OATAPT I 75,4), 
lOATAPTI  75,5)/  984.7,  6.3,  .85,  .048,  536.0/ 

OATA  LATI  76),L0NI  76),ELEV1I  76)/  57,43,  11.47,  5./ 

OATA  OATAPTI  76 , 1 ) , OATAPT  I 76 , 3) , DAT APT  I 76,3) .OATAPT  I 76,4), 
lOATAPTI  7e,5)/1013.4,  7.6,  .80,  .053,  565.0/ 

OATA  LATI  77),L0NI  77I,ELEV1I  77)/  57.40,  18,31,  47,/ 

OATA  OATAPTI  77 , 1 ) , DAT  APT  I 77 ,3) , OAT APT  I 77 , 3)  , OAT  APT  I 77,4), 
lOATAPTI  77,5)/1007,4,  7.1,  .65,  .039,  559.0/ 

OATA  LATI  78),LONI  78),ELEV1I  78)/  60.31,  33.16,  54./ 

OAT*  OATAPTI  78 , 1 ) , OATAPT I 78 , 3 ), OATAPT I 78,3) .OATAPT I 78,4), 
lOATAPTI  78,5)/1005.9,  4.6,  .80,  .043,  586.0/ 

OATA  LATI  79),L0NI  79),ELEW1I  79)/  60.19,  34.58,  58./ 

OAT*  OATAPTI  79,1) , OATAPTI  79 ,3) , OAT APT  I 79 ,3)  .OATAPT I 79,4), 
lOATAPTI  79,5)/1005,6,  4.4,  .80,  .057,  638.0/ 

OATA  LATI  80),LON(  eO),ELEVll  60)/  53.38,  13.34,  49./ 

OATA  OATAPTI  80,1) , OATAPTI  8C ,3) , 0 ATAPT I 80 , 3) , OAT APT  I 80,4), 
lOATAPTI  80,5) /1009. 3,  9,5,  .73,  ,109,  557.0/ 

OATA  LATI  OD.LONI  81),ELEV1I  81)/  49.30,  11.05,  313./ 

OATA  OATAPTI  81 ,1) , OATAPT I 81 ,3) , OATAPT I 8 1 ,3) , OAT APT  I 81,4), 
lOATAPTI  81,5)/  976.6,  8.4,  .78,  .185,  633.0/ 

OATA  LATI  83),L0NI  83),FLEV1I  63)/  48.08,  11.43,  539./ 

OATA  OATAPTI  83 , 1 ) , OATAPT I 83 , 3 ), OAT APT  I 83,3) .OATAPT I 83,4), 
lOATAPTI  83,5)/  954.3,  7.9,  .90,  .371,  964.0/ 

OATA  LATI  83),L0NI  03),ELEV1I  83)/  54.06,  13.37,  3./ 

OATA  OATAPTI  83 , 1 ) , 0 AT APT  I 83 , 3 ) , D AT APT  I 83,3) .OATAPT  I 83,4), 
lOATAPTI  83,5)/1014.3,  8.3,  .63,  .060,  553.0/ 

OATA  LATI  84),L0NI  84),ELEV1I  84)/  53.31,  13.05,  70./ 

OATA  OATAPTI  84 , 1 ) , 0 AT APT  I 84,3) .OATAPT I 84,3) , OATAPT  I 84,4), 
lOATAPTI  84,5)/  951.6,  8.0,  . 8C , .090,  571.0/ 

OATA  LATI  85),L0NI  85),ELEV1I  85)/  53.13,  14.07,  99./ 

OATA  OATAPTI  85 , 1) , OATAPT I 85 ,3 ), DAT APT  I 85,3) .OATAPT I 85,4), 
lOATAPTI  65,5)/1003. 7,  8.6,  .78,  .116,  543.0/ 

OATA  LATI  86),L0NI  86),ELEV1I  86)/  53.06,  11.35,  85./ 

OATA  OATAPTI  86 , 1) , OATAPTI  86 , 3) , OAT  APT  I 86 ,3) .OATAPT I 86,4), 
lOATAPTI  86,5) /1005. 3,  8,7,  .77,  .160,  515.0/ 

OATA  LATI  87),L0N(  67),ELEV1I  87)/  51.51,  10.46,  340./ 

OATA  OATAPTI  87 , 1 ) , 0 AT APT  I 87 , 3 ), OAT APT  I 87 ,3) .OATAPT I 87,4), 
lOATAPTI  87,5)/  986.6,  8.5,  .74,  .089,  637,0/ 

OATA  LATI  8e),LONI  86),ELEV1I  88)/  51.34,  13.34,  133./ 

OATA  OATAPTI  98 , 1) , OATAPT I 9 6 , 3 ) , D AT  APT  I 88,3) , OATAPT I 88,4), 
lOATAPTI  88,5)/  999.5,  8.7,  ,79,  .081,  549.0/ 

'lATA  LATI  89),L0NI  89),rLEVll  89)/  51.10,  14.57,  338./ 

OAT*  OATAPTI  89 , 1 ) , OATAPT I 89, 3) , OATAPT  I 89,3) , OATAPT  I 89,4), 
lOATAPTI  89,5)/  987,9,  8.0,  .78,  .159,  688.0/ 

OATA  LATI  90),LONI  90),pLE¥1I  90)/  50.59,  10.58,  316,/ 

OATA  OATAPTI  90 , 1 ) , 0 AT APT  I 90 , 3) , 0 AT APT  I 90 ,3)  .OATAPT I 90,4), 
lOATAPTI  90,5)/  976.5,  7.8,  .80,  .145,  533.0/ 

OATA  LATI  91),L0NI  91),ELEV1I  91)/  50.14,  10.33,  456./ 

OAT*  OATAPTI  91 , 1 ) , OATAPT I 91 , 3 ), OAT APT  I 91 , 3)  , OA T APT  I 91,4), 
lOATAPTI  91,5)/  963.3,  7.9,  .79,  ,175,  718.0/ 

OATA  LATI  93),L0NI  93),FLEV1I  93)/  50.39,  10.09,  494./ 

OATA  OATAPTI  93 , 1 ) , OATAPT I 93 ,3 ), OAT APT  I 93 ,3) , OAT  APT  I 93,4), 
lOATAPTI  93,5)/  957.7,  7.9,  .79,  .175,  718.0/ 
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DATA  LAT(  93),L0N(  93)«ELrVl(  931/  10.37.1152./ 

DATA  DATAPTt  93. 1 ) . OATAPT ( 93.2) .OATAPT ( 93.3) .OATAPT ( 93. A). 
lOATAPTI  93.5)/  862.1.  2.9.  .88.  .1A1.1A22.0/ 

DATA  LAT(  9A).L0N(  94).FLEV1(  94)/  50.26.  12.57.1215./ 

DATA  OATAPTt  94 . 1 ) . DATAPTl  94 .2) . OATAPT  ( 94.3) .OATAPT ( 94.4). 
lOATAPTI  94.5)/  076.2.  2.0.  .07.  .087.1109.0/ 

DATA  LATI  95).L0N<  95).ELEV1I  95)/  48.12.  15.38.  202./ 

DATA  OATAPT!  95 .1) .OATAPT I 95 .2) .OATAPT I 95.3) .OATAPT I 95.4). 
lOATAPTI  95.5)/  983.7.  8.7.  .77.  .146.  741.0/ 

DATA  LATI  96).L0NI  96).ELEV1I  96)/  48.15.  16.22.  212./ 

OATA  OATAPTI  96.1) .OATAPTI  96.2) .OATAPT  I 96.3) .OATAPT  I 96.4). 
lOATAPTI  96.5)/  991.2.  9.8.  .71.  .171.  615.0/ 

OATA  LATI  97).L0NI  97).ELEV1I  97)/  47.48.  13.00.  450./ 

DATA  OATAPTI  97.1) .OATAPTI  97 .2 ) .OATAPT  I 97.3) .OATAPT  I 97.4). 
lOATAPTI  97.5)/  963.2.  0.1.  .77.  .164.1278.0/ 

OATA  LATI  90).L0NI  98).ELEV1I  98)/  47.00.  15.27.  347./ 

OATA  OATAPTI  90.1) .OATAPTI  98.2 ). OATAPT I 98.3) .OATAPT I 98.4). 
lOATAPTI  98.5)/  975.3.  8.3.  .78.  .158.  840.0/ 

OATA  LATI  99).L0NI  99).ELEyil  99)/  50.05.  12.24.  471./ 

OATA  OATAPTI  99.1) .OATAPTI  99.2) .OATAPT  I 99.3) .OAT APT  I 99.4). 
lOATAPTI  99.5)/  959.7.  7.0.  .81.  .150.  540.0/ 

OATA  LATI100).LONI100).ELEV1I100)/  50.06.  14.17.  381./ 

OATA  OATAPT II 00.1) .OAT APT  1100. 2) .OATAPT 1100 .3) . OATAPT 1100.4) . 
10ATAPTI100.5)/  970.0.  7.9.  .78.  .160.  520.0/ 

OATA  LATIlOl) tLONIlOl) .ELEVIIIOI)/  49.47.  18.16.  253./ 

OATA  OATAPTIlOl.l) . OATAPT 1101 .2) • OATAPT I 101.3) .OATAPT 1101.4) . 
10ATAPTI101.5)/  985.7.  8.1,  .79.  .175.  659.0/ 

OATA  LATI1B2).LONI102).ELEV1I102)/  49.09,  16.42,  242./ 

OATA  OATAPTI102.1) .OATAPTI102,2) ,OATAPTI102,3) .OATAPT 1102. 4) , 
10ATAPTI102,5)/  907.1.  8.8.  .73,  .186,  527.0/ 

OATA  LATI103).LONI103).ELEV1 (103)/  53.55.  14.14.  10./ 

OATA  OATAPT 1103,1) .OATAPT 1 103.2). OATAPT  1103. 3) .OATAPT 1103,4) , 
10ATAPTI103.5)/1014.2.  8.3.  .03.  .060.  625.0/ 

OATA  LAT(104).LONI104),ELEV1I104)/  53.24,  14.37,  7./ 

OATA  OATAPT 1104.1) .OATAPT I104.2).OAT APT  1 104,3) .DAT  APT (104.4) , 
10ATAPT(104.5)/1013.0.  0.3,  .81.  .080.  511.0/ 

OATA  LAT(105) .LONIIOS) .ELEV1(105) / 54.12,  16.11,  34./ 

OATA  OATAPT 1105,1). OAT APT  1 105. 2) .OATAPT 1105.3) .OATAPT (105,4) , 
10ATAPT(105.5)/1010.2.  7.5,  .82.  .065,  697.0/ 

OATA  LATI106).LON(1Q6).ELEV1(106)/  54.31,  10.33.  5./ 

OATA  OATAPTIlOS.l) .OATAPT(106.2) .OATAPT(106.3) .OATAPT (106.4) . 
10ATAPT(106,5)/1014.2.  7.7,  .81,  .092.  750.0/ 

OATA  LATI107),LON(107).ELEV1I107)/  54.23.  18.36,  12./ 

OATA  OATAPT (107.1) .OATAPTI 107,2) .OATAPT (107,3) .OATAPT  1107,4) , 
10ATAPT(107.5)/1013.3.  7.7,  .79.  .116,  499.0/ 

OATA  LAT(109I.LON(108).ELEV1(109)/  52.25,  16.50.  92./ 

OATA  OATAPT 1108.1) .OAT APT (108.2). OATAPT (108,3) .OATAPT (108.4) , 
10ATAPT(108.5)/1004.7,  8.0.  .70.  .120,  480.0/ 

OATA  LAT(109) ,LON(109) .ELEV1(109) / 52.09.  20.59.  107./ 

OATA  OATAPTI109.1) ,OATAPT(109.2) .OATAPT (109.3) .DAT APT (109.4) , 
10ATAPT(109.5)/1002.9.  7.0.  .78.  .175,  471.0/ 

OATA  LAT(110),LON(119).ELEVllli:)/  53.06,  23.10.  151./ 

OATA  OATAPTIllO.l) . DATAPTIllO .2) .OATAPT 1 110 .3)  .OATAPT (110 ,4)  , 
10ATAPT(110.5)/  997.2.  6.0,  .80,  .165.  560.0/ 

OATA  LATIllll.LONIllD.ELEVlIlll)/  51.06.  16.23.  124./ 

OATA  OATAPT (111,1) .OATAPT I 111, 2)  .OATAPTI 111. 3) .OATAPT (111.4)  , 
10ATAPT(111.5)/1001.0.  8.3,  .70.  .168.  557.0/ 

OATA  LAT(112) .L0N(112).ELEV1 (112)/  51.03.  17.21.  134./ 

OATA  0ATAPT(112.1) .0ATAPT(112.2) .OATAPT (112.3) .OAT APT (112 ,4) , 
10ATAPT(112.5)/  999.7.  0.0.  .70.  .175.  525.0/ 
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DATA  LAT(113).L0N(113I«ELEV1(113)/  50. OS,  19.40,  237./ 

DATA  OATAPT 1113,1) ,OATAPT (113,21 ,0 AT APT (113,31 ,OATAPr (113,41  , 
10ATAPT(lt3,S)/  907.9,  7.0,  .00,  .160,  S7S.0/ 

DATA  LAT(114) ,L0N(114I ,fLP¥l (114)/  49.40,  22.46,  280./ 

DATA  0ATAPT(114,1) ,0ATAPT(114,2),DATAPT(114,3) .OATAPT (114,4) , 
iOATAPT(114,S)/  902.7,  0.0.  .77,  .138,  643.0/ 

DATA  LAT(11S),L0N(11S),ELEV1(11S)/  46.47,  23.34,  41S./ 

DATA  OAf APT ,UATAPT(11S,2) ,0ATAPT(11S,3) , OATAPT (lis, 4), 
10ATAPT(115,5)/  966.3,  0.6,  .74,  .132,  507.0/ 

DATA  LAT(116),L0N(116),ELEV1 (116)/  47.09,  24.31,  366./ 

DATA  0ATAPT(116,1) ,nATAPT(lie,2) , OATAPT (116,3) , OAT APT (116,4) , 
10ATAPT(116,5) / 971.3,  8.5,  .75,  .101,  684.0/ 

DATA  LAT(117) ,L0N(117) ,FLEV1 (117)/  59.50,  30.10,  4./ 

DATA  OATAPT (117, 1 ), OATAPT ( 117 ,2) , OATAPT (117,3) , DAT  APT  (117,4) , 
1DATAPT(117,5)/1Q13.2,  4.6,  .79,  .069,  525.0/ 

DATA  LAT(llO)  ,L0N(118)  .FUEVKllO)  / 59.25,  24.40,  44./ 

DATA  OATAPT (118,1 ), OATAPT (110,2) , OATAPT (118,3) , OATAPT (110,4) , 
10ATAPT(lie,5)/1007.6.  5.2,  .81,  .050,  568.0/ 

DATA  LAT(119),L0N(119),ELEV1(119)/  54.53,  23.53,  75./ 

DATA  OAT APT (119,1) , OAT  APT ( 119 , 2) , OAT APT ( 119,3)  , OAT APT (119,4)  , 
inATAPT(119,5)/1005.8,  6.0,  .80,  .150,  625.0/ 

DATA  LAT(120),LON(120),ELEV1(120)/  53.52,  27.32,  234./ 

DATA  OATAPT (120,1 ), OATAPT (120,2 ) .OATAPT (120,3) , OATAPT (120,4) , 
10ATAPT(120,5)/  907.0,  <^.3,  .80,  .172,  506.0/ 

DATA  LAT(121) ,L0N(121),ELEV1 (121)/  56.23,  30.36,  99./ 

DATA  OATAPT (121,1) .OATAPT (121, 2) .OATAPT (121,3) .OATAPT (121,4) , 
IDATAPT (121,5)/1003.  2,  3.0,  .01,  .115,  430.0/ 

DATA  LAT(122) ,LON(122) ,FLEV1(122) / 59.17,  39.52,  118./ 

DATA  OATAPT (122,1) .OATAPT (122,2) .OATAPT (122,3) .OATAPT (122,41 , 
10ATAPT(122,5)/1000.2,  2.6,  .80,  .110,  390.0/ 

DATA  LAT(123),U0N(123),ELEV1(123)/  58.39,  49.37,  164./ 

DATA  0ATAPT(123,1( ,0ATAPT(123,2) ,OATAPT(123,3) .OATAPT (123,4)  , 
1DATAPT(123,5)/  995.7,  1.8,  .78,  .154,  526.0/ 

DATA  LAT(124),L0N(124),ELEV1(124)/  55.45,  37.34,  156./ 

DATA  OATAPT (124,1) .OATAPT (124,2) .OATAPT (124, 3) .OATAPT (124,4) , 
10ATAPT(124,5) / 997.0,  4.4,  .77,  .210,  560.0/ 

OATA  LAT(125) ,LON(125),FLE¥l (125)/  55.47,  49.11,  64./ 

OAT A OATAPT (125,1) , OATAPT (125 , 2) , OAT APT ( 125 ,3) .OATAPT (125,4)  , 
10ATAPT(125,5)/1009.5,  3.1,  .75,  .205,  435,0/ 

DATA  LAT(126) ,LON(126) ,FLEV1 (126)  / 51.42,  39.10,  164./ 

OATA  0ATAPT(126,1) .OATAPT (126,2 >,0 AT APT (126, 3) , OAT APT (126,4) , 
10ATAPT(126,5)/  997.1,  5.5,  .73,  .168,  550.0/ 

OATA  LAT(127) ,L0N(127),5LEV1 (127)/  51.34,  46.02,  156./ 

DATA  OATAPT (127,1 ) ,0ATAPT(127, 2) .OATAPT (127,3) .OATAPT (127,4)  , 
1DATAPT(127,5)/  998.2,  5,9,  .75,  .138,  391.0/ 

OATA  LAT(129),L0N(129),FLEV1(128)/  49.49,  23.57,  325./ 

OATA  OATAPT (128,1 ) .OATAPT (128, 2),0 AT  APT ( 12 8, 3) , OA T APT (128 ,4) , 
10ATAPT(128,5)/  977.2,  7.0,  .80,  .124,  655.0/ 

DATA  LAT(129) ,L0N(129) ,ELFV1 (129) / 50.24,  30,27,  179./ 

OATA  OATAPT (129,1) .OATAPT (129,2) , OATAPT ( 129 ,3) .OATAPT  (129,4) , 
lOATAPT  (129,5)  / <>95,2,  7.4,  .76,  .137,  615,0/ 

DATA  LAT(13Q),LON(130),ELFV1(13Q)/  41.39,  .53,  233./ 

OATA  OATAPT (130, 1) ,DATAPT( 130, 2) , OAT APT ( 130 , 3) .OATAPT (130 ,4) , 
10ATAPT(130,5)/  968.4,  14.8,  .61,  .158,  339.0/ 

OATA  LAT(131),L0N(131),ELEV1(131)/  36.50,  2.28,  7./ 

OATA  OAT APT (131,1) .OAT APT (131, 2) .OATAPT ( 131 ,3) .OATAPT (131 ,4) , 
10ATAPT(131,5)/1015.9,  18.1,  .73,  .139,  227,0/ 

OATA  LAT (132) ,LON (132) ,ELFV1 (132) / 39.57,  32.53,  894./ 

OATA  OAT  APT (132,1) , 0 A T APT (1 32 ,2 ) , 0 AT APT ( 132 , 3) , DAT APT(1 32 ,4) , 
10ATAPT(132,5)/  911.7,  11.7,  .60,  .081,  446.0/ 
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DATA  LATI133)tL0Na33),ELEVl(133)/  38.43,  35.29,1084./ 

9ATA  OATAPTll 33,11 ,0ATAPT(133,2) ,OATAPT (133,3) ,OATAPT (133,4) , 
10ATAPT(133,5) / •91.0,  10.8,  .64,  .285,  352.0/ 

DATA  LA7(134),Li.^(134),ELEVl(134)/  39.45,  37.01,1284./ 

DATA  0ATAPT(134,1) ,0ATAPT(134,2),0ATAPT(134,3) .OATAPT (134,4) , 
1DATAPT(134,5)/  867.7,  8.4,  .65,  .250,  414.0/ 

3ATA  LAT(135) ,L0N(135) ,PLEV1(135) / 39.44,  39.30,1215./ 

DATA  OATAPT (135,1) .OATAPT (135,2), OATAPT (135, 3) .OATAPT (135,4) , 
10ATAPT(135,5)/  881.7,  10.7,  .59,  .117,  363.0/ 

DATA  LAT(136),L0N(136),ELEV1(136)/  38.45,  30.32,1034./ 

DATA  OATAPT (136,1) .OATAPT (136 ,2) ,0 ATAPT ( 136 ,3) .OATAPT (136,4) , 
10ATAPT(136,5)/  986.6,  11.2,  .63,  .066,  621.0/ 

9ATA  LAT(137) ,L0N(137) ,5LrVl(137)/  37.45,  30.33,  997./ 

DATA  OATAPT (137,1) .OATAPT (137, 2), OAT APT (137, 3) .OATAPT (137, 4), 
10ATAPT(137,5)/  900.5,  12.1,  .62,  .244,  617.0/ 

OATA  LAT(138) ,LON(138) ,ELEV1(138)/  37.58,  32.33,1022./ 

OAT A OATAPT (1 38,1 ) , OA TAPT (138 ,2 ) , 0 ATAPT ( 138 ,3) , OAT APT (138,4) , 
10ATAPT(138,5)/  987.7,  11.5,  .60,  .078,  311.0/ 

OATA  LAT(139) ,LON(139) .ELEV1(139)/  37.08,  30.46,  547./ 

OATA  0ATAPT(139,1) , OATAPT( 139,2) , OATAPT (139,3) .OATAPT (139,4) , 
10ATAPT(139,5) / 946.3,  18.1,  .48,  .183,  464.0/ 

OATA  LAT(140),LON(140),rLEVl(140)/  30.26,  38.05,  849./ 

OATA  OATAPT (140,1) .OATAPT (140,2), OATAPT (140,3) .OATAPT (140,4) , 
10ATAPT(140,5)/  915.3,  13.3,  .54,  .058,  370.0/ 

OATA  LAT(141) ,L0N(141) ,ELEV1(141)/  37.53,  40.11,  677./ 

OATA  0ATAPT(141,1) ,0ATAPT(141,2),0ATAPT(141,3) .OATAPT (141,4), 
10ATAPT(141,5)/  932.3,  15.7,  .52,  .104,  489.0/ 

OATA  LAT(142) ,L0N(142) ,ELEV1(142) / 38. 2^,  43.19,1661./ 

OATA  OATAPT (142,1 ) ,0ATAPT(142, 2) .OATAPT (142,3) .OATAPT (142,4) , 
10ATAPT(142,5)/  827.7,  8.9,  .59,  .040,  385.0/ 

OATA  LAT(143) ,L0N(143),ELEV1(143)/  39.55,  41.16,1869./ 

OATA  OATAPT(143,l) , OATAPT ( 143 ,2) , OATAPT ( 143 , 3) , OAT APT ( 143 , 4) , 
10ATAPT(143,5)/  806.7,  5.1,  .61,  .113,  460.0/ 

OATA  LAT(144),L0N(144),ELEV1(144)/  49.56,  36.17,  152./ 

OATA  OATAPT (144, 1 ), OATAPT ( 144 , 2) , OAT APT ( 144 ,3) , OAT APT (144,4)  , 
10ATAPT(144„5)/  998.8,  7.1,  .7*.,  . 110,  573.0/ 

OATA  LAT(145) ,L0N(145),ELEV1(145)/  46.29,  30.08,  64./ 

OATA  0ATAPT(145,1) ,0ATAPT(145,2) ,OATAPT(145,3) .OATAPT (145,4) , 
10ATAPT(145,5)/1006.9,  9.9,  .76,  .100,  389.0/ 

OATA  LAT(146) ,L0N(146),ELEV1(146) / 45.01,  33.59,  205./ 

OATA  OAT APT (146,1) , OATAPT ( 146, 2) , OAT APT ( 146 ,3) , OAT  APT (146 ,4) , 
10ATAPT(146,5)/  990.6,  10.1,  .74,  .168,  519.0/ 

OATA  LAT(147),L0N(147),FLEV1I147)/  47.15,  39.49,  77./ 

OATA  0ATAPT(147,1) ,0ATAPT(147,2) .OATAPT (147,3) , OAT APT ( 147 ,4) , 
10ATAPT(147,5)/1007.6,  9.0,  .72,  .090,  487.0/ 

OATA  LAT(148),L0N(148),ELFV1(148)/  46.16,  48.02,  -18./ 

OATA  OATAPT (148,1 ), OATAPT) 14 8,2) , OAT APT(148,3) .OATAPT (14  8,4)  , 
10ATAPT(148,5)/1019.6,  9.5,  ,79,  .085,  179.0/ 

OATA  LAT(149),L0N(149),ELPV1(149)/  44,03,  43.02,  500./ 

OATA  OAT APT (149,1 ) ,0ATAPT(149, 2) .OATAPT (149, 3) , OAT APT (149 ,4) , 
lOATAPT (149,5)/1017.&,  8.6,  .76,  .108,  482.0/ 

OATA  LAT(150) ,LON(150) ,ELEV1(150) / 45.43,  4.37,  201./ 

OATA  OAT APT (150,1) , OAT APT ( 150 , 2) , 0 AT APT ( 150 , 3) .OATAPT (150 ,4) , 
10ATAPT(150,5)/  993.0,  11.4,  .76,  .178,  513.0/ 

OATA  LAT(151),L0N(151),rLEVl(151)/  41.41,  44,57,  490./ 

OATA  OATAPT (151,1) .OATAPT (151,2) , 0 AT APT ( 15 1 , 3) .OATAPT (151 ,4) , 
nATAPT(151,5) / 959.7,  12.9,  .66,  .145,  496.0/ 

OATA  l.AT(152)  ,LON(152)  ,FLFV1  (152)/  47.15,  9.21,2500./ 

OATA  OATAPT (152,1) .OATAPT (152, 2) .OATAPT (152,3) .OATAPT (152,4), 
10ATAPT(152,5)/  748.7,  -1.9,  ,77,  .187,2478.0/ 
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0»T*  L«T(153I ,L0N(1?3I ,rLEVl(l?3 
OATJk  0ATAPT(1S3«1)  ,OATAPT(lS3.e) 
10ATAPT(153«5)  / 949. 6t  S.5,  .7(>, 
OATA  LAT(194) «L0N(154)«ELeVl (1S4 
DATA  0ATAPT(1S4«1) ,0ATAPT(154«?) 
10ATAPT(154.5I  / 966.0.  10. 3.  .76. 
DATA  LAT<195) .LON ( 1 55 ) .ELEVl ( 155 
DATA  0ATAPTI155.1) .0ATAPT(155.2> 
nATAPTC155,5»/  959.3.  10,0.  .7f, 
DATA  LATI156) .L0NI156).FLEV1(156 
DATA  OATAPT(156.1) .0ATAPT(156.7) 
10ATAPTI156.5)/  982.9.  11.7,  .68. 
DATA  LAT(157) ,L0N(157) .SLEV1I157 
DATA  DATAPT(157,1) ,0ATAPT(157,2) 
10ATAPT(157,5)/  968.7.  8.8.  .82, 

DATA  LAT(158i ,LON(158>,ELeVl<15e 
DATA  0ATAPT(15«,1) ,0ATAPT(158,2) 
10ATAPT<158.5»/  709.8,  -4.7,  .82, 
DATA  LAT(159) .LON ( 159 ) .ELEVl I 159 
DATA  OATAPT(159,ll  , 0ATAPT«159,2) 
1DATAPTJ159.5) / 946.0,  *.6.  .72, 

DATA  LAT(160) .LON(160I.flEV1(160 
DATA  OATAPT(160.1I ,OATAPT(160,2> 
10ATAPT(160.5)/1016.2.  14.2,  .78. 
DATA  LAT(161).LON(161)  .ELEVKiei 
DATA  0ATAPT(161,1I ,0ATAPTfl61.2) 
10ATAPTI161.5)/1010.3,  15.2,  .67, 
DATA  LAT(162) ,L0N I 16? ) . ^LE VI ( 162 
DATA  DATAPT(162.1) .0ATAPT(162.2) 
1DATAPT(162,5)/1008.6.  14.2,  .66, 
DATA  LATI163),L0NI163I.ELEV1I163 
DATA  0ATAPT(163,1) .DATAPT(163.2) 
1DATAPT(163. 51/1014. 6,  14.2,  .71, 
DATA  LATI164),L0NI164),ELEV1I164 
DATA  0ATAPT(164,1I ,QATAPr(164,2l 
10ATAPT(164,5»/1013.3,  14.8,  .76, 
DATA  LATI165) ,L0NC165|,ELEV1(165 
DATA  0ATAPTI165,1) ,0ATAPT(165,2) 
10ATAPTa65,5»/1013.9,  14.7,  .79, 
DATA  LAT(16E) ,L0N(156) ,rLEVl(166 
DATA  0ATAPT(166,1> ,DATAPTI166,2) 
1DATAPTI166,5)/1017.0,  18.2,  .74, 
DATA  LAT(167) ,L0N(167),ELFV1 (167 
DATA  0ATAPT(167,1)  ,0ATAPT(t67,2) 
10ATAPT(167,5)/  936.8,  11.6,  .71, 
DATA  LAT(168),L0N(16<«),ELFV1  (168 
DATA  0ATAPT(168,1) ,DATAPT(t6e,2) 
19ATAPT(168,5) /1008.5,  14.4,  .7*, 
DATA  LAT(169) ,L0N(169),ELEV1 (169 
DATA  0ATAPT(169,1>  ,0ATAPT(16O,2» 
IDATAPT (169,5)/1000.8,  15.9,  .71, 
DATA  LAT(170) ,LON(170>,FLEV1(170 
DATA  OATAPT(170,H ,OATAPT(170,2) 
1DATAPT(170,5)/1005.1,  16.6,  .7C, 
DATA  LAT(171),L0N(17H,FLFV1(171 
DATA  0ATAPT(171,H  ,DATAPT(171,2) 
lDATAPT(iri,5)  / 987.7,  16.2,  .64, 
DATA  LAT(172) ,L0N(172),ELFV1(172 
DATA  0ATAPT(172,1) ,DATAPT(172,2I 
1DATAPT(1’2,5) /1015. 3,  17.8,  .72, 


/ 47.23,  8.34,  569./ 

OATAPT(153,3) ,0ATAPT(153,4I, 
.091,1071.0/ 

/ 46.15,  6.07,  430./ 

0ATAPT(154,3) ,OATAPT (154,4) , 
.126,  831.0/ 

/ 46.49,  6.57,  491./ 

OATAPT (155,3) ,OATAPT  (155,4) , 
.100,  978.0/ 

/ 46.00,  8.58,  276./ 

0ATAPT(156,3) , OAT APT ( 156 ,4) , 
.188,1744.0/ 

/ 47.39,  9.29,  407./ 

OATAPT (157,3) , OATAPT  (157,4) , 
.099,  960.0/ 

/ 47.25,  10.59,2962./ 

OATAPT (158,3) ,OAT APT (158,4) , 
.150,1351.0/ 

/ 47.16,  11.21,  598./ 

OATAPT (159,3) ,OAT APT  (159,4) , 
.083,  899.0/ 

/ 41.02,  40.30,  4./ 

OATAPT(160,3)  ,OAT APT  (160 ,4) , 
.202,2441.0/ 

/ 42.44,  2.52,  48./ 

OATAPT (161,3) , OATAPT (161,4) , 
.359,  639.0/ 

/ 43.52,  4.24,  62./ 

0ATAPT(162,3) .OATAPT  (162 ,4) , 
.224,  743.0/ 

/ 43.27,  5.13,  8./ 

OATAPT(163,3) .OATAPT  (163,4) , 
.227,  550.0/ 

/ 43.39,  7.12,  10./ 

0 ATAPT( 164, 3) .OATAPT (164,4 ) , 
.102,  862.0/ 

/ 41.55,  8.18,  9./ 

OATAPT (165,3) .OATAPT (165,4) , 
.106,  681.0/ 

/ 36.09,  -5.21,  3./ 

OATAPT (166,3) .OATAPT  (166,4) , 
.169,  815.0/ 

/ 41.49,  -6.46,  692./ 
0ATAPT(167,3) .OATAPT ( 167 ,4) , 
.189,  973.0/ 

/ 41.14,  -8.41,  73./ 

OATAPT (168,3) .OATAPT  (168,4), 
. 195,1150.0/ 

/ 40.12,  -8.25,  140./ 

OATAPT (169,3) .OATAPT (169,4) , 
.234,  961.0/ 

/ 38.46,  -9.08,  110./ 

OATAPT  (170,3) , DAT APT  ( 170 ,4) , 
.074,  705,0/ 

/ 38.01,  -7.52,  247./ 

OATAPT (171,3) .OATAPT  (171,4) , 
.125,  549.0/ 

/ 37.01,  -7.58,  9./ 

OATAPT (172,3) .OATAPT (172 ,4) , 
.067,  455.0/ 
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DATA  LAT(173I,L0N(173),ELEV1(173>/  41.39t  -4.43t  715./ 
lATA  OATAPT (173.1 ) .OATAPT (173.71 .OATAPT (173.3) .OATAPT (173.4) . 

1DATAPT(173.S) / 933.6.  12.7.  .64.  .099.  409.0/ 

DATA  LAT(174).L0N(174).7LPV1(174)/  40.74.  *3.41.  657./ 

DATA  OATAPT (174.1) .OATAPT (174.2) .OATAPT (174.3) .OAT APT (174.4) . 

10ATAPT(174.5)/  940.5.  13.9.  .62.  .126.  436.0/ 

OATA  LAT(175) .LON(175) .FLEV1(175) / 36.53.  -6.49.  192./ 

OAT A OATAPT (175.1 ) .0ATAPT(175.2) . OATAPT ( 175.3) .CAT APT (175 .4) . 

10ATAPT(175.5)/  994.6.  16.7.  ,61.  .092.  474.0/ 

OATA  LAT(176).L0N(176).ELEV1(176)/  37.22.  -6.00.  13./ 

OATA  0ATAPT(17E.l) .0ATAPT(17E ,2) , OAT APT ( 176 ,3) .OATAPT (176.4)  . 

10ATAPT(176.5)/1016.0.  18.0.  .66.  .148.  559.0/ 

OATA  LAT(177).L0N(177).ELEW1(177)/  36.22.  -.3C.  82./ 

OATA  OATAPT(177.1) .0ATAPT(177.2).0ATAPT(177.3) .nATAPT(177.4). 

10ATAPT(177.5) /1007.2.  17.8.  .66.  .230.  302.0/ 

OATA  LAT(178).L0N(178).ELEV1(178)/  41.24.  2.09.  95./ 

OATA  OATAPT (176.1) . OAT  APT ( 176 .2 ) .OATAPT (178.3) .OATAPT (176.41 . 

10ATAPT(178*5)/1005.2.  16.5.  .70.  .363.  596.0/ 

OATA  LAT(179) .L0N(179) .FLEV1(179) / 39.36.  2.42.  45./ 

OATA  OATAPT (179. 1) .OATAPT( 179.2) .OATAPT (179.3) .OATAPT (179.4) . 

10ATAPT(179.5)/1010.3.  16.8.  .75.  .157.  447.0/ 

OATA  LAT(100).LON(100).ELEVl(ieO)/  39.52.  4.16.  59./ 

OATA  OATAPT(lOO.l)  .OATAPT(100.2) .OATAPT(100.3) .OATAPT(160.4)  . 

10ATAPT(100.5)/1008.9.  16.6.  .72.  .352.  637.0/  ] 

OATA  LAT(101).LON(181).ELrVl(161)/  39.37.  19.55.  2./ 

OATA  OATAPT (181.1) . OATAPT (181 . 2) . OATAPT ( 181 .3) . OAT APT ( 1 81 . 4)  . 

10ATAPT(181.5)/1014.1.  17.6,  .71.  .125.1311.0/ 

OATA  LAT(1B2) .L0N(182).ELEV1 (102)/  30.15.  21.44.  3./ 

OATA  OATAPT (182.1) .OATAPT (182.2) .OATAPT (102. 3) . OAT  APT (1 02 . 4) . 

10ATAPT(182.5)/1014.3.  17.6.  .69.  .164.  753.0/ 

OATA  LAT(183)  .L0N(183)  .ELEVl  (183)  / 37. <.7.  20.53.  6./ 

OATA  OATAPT (183.1) . OAT APT ( 183 .2) , OAT APT ( 103. 3) .OATAPT (103.4)  . 

10ATAPT(183.5)/1014.2.  18.1.  .69,  .175,  933.0/ 

OATA  LAT(ie4),LON(104),FLEVl(104)/  37.04,  22.01,  0./ 

OATA  0ATAPT(184,1) .0ATAPT(184,2) .0ATAPT(184,3) .0ATAPT(104,4) , 

10ATAPT(184.5)/1014.0,  16.6,  .68.  .210.  813.0/ 

OATA  LAT(185) .LON(195) .ELEVl (165)/  40.31,  22.56.  4./ 

OATA  OATAPT(185.1),OATAPT)165,2),OATAPT(105,3).OATAPT(105.4), 

10ATAPT(185,5)/1015.4,  16.1.  .66.  .101.  465.0/ 

OATA  LAT(186),L0N(186).ELEV1(166)/  39.38,  22.25,  74./ 

OATA  OATAPT(196.1).OATAPT(186.2).OATAPT(186.3),OATAPT(106,4), 
lOATAPT (166.51/1005.9.  16.2.  .67.  .120.  400.0/ 

OATA  LAT(107),LON(187),ELEV1(187I/  37.58,  23.43.  107./ 

OATA  0ATAPT(187.1)  ,0ATAPT(187, 2) .OATAPT (187,3) .OATAPT (187,4)  , 

10ATAPT(187,5)/1002.1.  17.8,  .63.  .167,  403.0/  i 

OATA  LAT(186).LON(188).FLEVl(ie6)/  35.20,  25.11.  39./ 

OATA  OATAPT(180,1).OATAPT(188.2).OATAPT(180,3),OATAPT(188.4),  j 

10ATAPT(188.5) /lOlO.O.  18.9.  .65.  .320.  441.0/  I 

OATA  LAT(189).L0N(189),ELFV1(189)/  35.09,  33.17.  220./ 

OATA  OATAPT (189.1) .OATAPT (189.2) .OATAPT ( 189, 3) , OAT  APT ( 1 89, 4)  , 

10ATAPT(189.5)/  987.0,  19.1,  .61.  .280,  339.0/ 

OATA  LAT(190).LON(190).^LEV1(190)/  34.41.  32.49.  113./ 

OATA  OATAPT (190 .1) , OAT APT (190 . 2 ) .OATAPT (190 .3) . DAT  APT (190 . 4)  , 

10ATAPT(190.5)/  999.0.  20.0.  .61.  .245.  450.0/ 

OATA  LAT(191)  ,L0N(191)  .‘^LEV1(191)  / 41.17,  36.20.  44./ 

OATA  0ATAPT(191,1) ,0ATAPT(19l, 2) , 0ATAPT(191, 3) .OATAPT (191.4)  , 

10ATAPT(191,5)/1010.2.  14.3,  .73,  .045,  731.0/ 

OATA  LAT(192),L0N(19?). ELEVl (192)/  41.40.  26.34.  48./ 

OATA  OAT APT (192,1) .OATAPT (192, 2) .OATAPT (192, 3) , OA T AFT (192 ,4)  , 

10ATAPT(192.5)/1011.0.  13.4,  .71,  .090,  727.0/ 
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OATA  LATI193)«L0NI193I«ELEV1(193>/  AO. 58.  29.05.  40./ 

DATA  0AtAPT(193.1) .OATAPTf 193,2) ,DATAPT(193,3)  .OAT APT (193 .4)  , 
10ATAPT(193,5)/1011.1,  12.6.  .75.  .070,  686.0/ 

DATA  LAT(194).L0N(194),ELEV1I194)/  40.11,  29.04,  100./ 

DATA  0ATAPT(194,1I ,OATAPT(194,2) ,0ATAPT(194,3) .OATAPT (194,4) , 
10ATAPT(194,5)/1002.6,  14.4,  .69,  .075,  726.0/ 

DATA  LAT(195),L0N(195),ELEV1(195)/  40.08,  26.24,  3./ 

DATA  0ATAPT(195,1) ,0ATAPT(195,2) ,0ATAPT(195,3> .OATAPT (195,4) , 
10ATAPr(195,5)/1019.2,  14.7,  .70,  .150,  608.0/ 

OATA  LAT(196),L0N(196),ELEV1(196)/  30.26,  27.10,  25./ 

OATA  0ATAPT(196,1) .OATAPT (196, 2) .OATAPT (196,3) .OATAPT (196,4) , 
10ATAPT(196,5)/1010.3,  17.5,  .62,  .271,  694.0/ 

OATA  LAT(197),L0N(197),FLEV1(197)/  36.42,  30.44,  50./ 

OATA  OATAPT (197,1 ) ,0ATAPT(197, 2) .OATAPT (197,3) .OATAPT (197,4) , 
10ATAPT(197,5)/1006.6,  10.6.  .65,  .453,1038.0/ 

OATA  LAT(190),LON(199),ELEVl(19e)/  37.00,  35.25,  66./ 

OATA  OATAPT (198,1) ,0ATAPT(198, 2) .OATAPT (190,3) .OATAPT (190,4)  , 
10ATAPT(190,5)/ia04.6,  18.6,  .65,  .335,  638.0/ 

OATA  LAT(199) ,L0N(199),ELEV1(199)/  41.22,  33.46,  799./ 

OATA  0ATAPT(199,1) ,0ATAPT(199, 2) .OATAPT  (199,3) .OATAPT (199,4) , 
10ATAPT(199,5)/  921.7,  9.7,  .70,  .206,  440.0/ 

OATA  LAT(200),LON(200),ELEV1(200)/  37.12,  28.21,  646./ 

OATA  OATAPT (200,1) , OATAPT( 200 ,2)  .OATAPT (200,3) .OATAPT (200,4)  , 
10ATAPT(200,5)/  937.3,  14.9,  .60,  .420,1189.0/ 

OATA  LAT(201),LON(201),ELEV1(201)/  45.26,  9.17,  103./ 

OATA  OATAPT (201,1 ) ,OATAPT( 20 1,2) .OATAPT (201, 3) .OATAPT (201,4) , 
10ATAPT(201,5)/1003.1,  12.3,  .79,  .141,  902.0/ 

OATA  LAT(202) ,LON(202),ELEV1(202)/  45.23,  10.52,  60./ 

OATA  OATAPT (202,1) .OATAPT (202, 2) .OATAPT (202,3) .OATAPT (202,4)  , 
10ATAPT(202,5)/1008.0,  12.4,  .78,  .202,  755.0/ 

OATA  LAT(203),LON(203),ELEV1(203)/  43.40,  10.23,  1./ 

OATA  OATAPT (203, 1) ,OATAPT(203 ,2) .OATAPT (203 ,3) .OATAPT (203,4) , 
10ATAPT(203,5)/1014.7,  14.9,  .77,  .209,  935.0/ 

OATA  LAT(204),LON(2Q4),ELEV1(204)/  41.48,  12.14,  3./ 

OATA  OATAPT(204,1) ,OATAPT( 204, 2) .OATAPT (204,3) .OATAPT (204,4)  , 
10ATAPT(204,S)/1014.2,  15.9,  .74,  .162,  749.0/ 

OATA  LAT(205),LON(205),ELEV1(205)/  40.57,  14.18,  72./ 

OATA  OATAPT(205,1) ,OATAPT(205, 2) .OATAPT (205,3) .OATAPT (205,4), 
10ATAPT(205,5)/1006.4,  15.8,  .73,  .109,  922.0/ 

OATA  LAT(206),LON(206),ELEVl(20e)/  45.39,  13.45,  20./ 

OATA  OATAPT (206,1) .OATAPT (206,2) .OATAPT  (206,3) ,OAT APT (206,4) , 
10ATAPT(206,5)/1013.0,  14.5,  .67,  .113,  965.0/ 

OATA  LAT(207) ,LON(20^) ,ELEV1(207) / 45.30,  12.20,  6./ 

OATA  OATAPT(207,1) , OATAPT(207, 2) .OATAPT (207,3) .OATAPT (207,4) , 
10ATAPT(207,5)/1015.1,  12.8,  .76,  .130,  854.0/ 

OATA  LAT(208) ,LON(208) ,ELEV1(208)/  45.37,  13.31,  104./ 

OATA  OATAPT(208,1)  ,OATAPT(208,2) .OATAPT (208 ,3) .OATAPT (208,4) , 
10ATAPT(200,5)/1002. 9,  14.8,  .72,  .064,  674.0/ 

OATA  LAT(209) ,LON(209) .ELEVl (209) / 42.26,  14.12,  16./ 

OATA  OATAPT(209,1) ,OATAPT(209,2) .OATAPT (209,3) .OATAPT (209,4) , 
10ATAPT(209,5)/1013.2,  14.4,  .70,  .141,  678.0/ 

OATA  LAT(210),LON(21!)),PLEV1  (210)/  40.26,  16.53,  12./ 

OATA  OATAPT (210,1) .OATAPT (210,2) , OAT  APT (21 0 ,3) ,0 AT  APT (2 10,4)  , 
10ATAPT(210,5)/1012.8,  17.1,  .70,  .082,  397.0/ 

OATA  LAT(211) ,L0N(211) ,5LEV1(211) / 40.39,  17.57,  10./ 

OATA  OATAPT (211,1) .OATAPT (21 1,2) .OATAPT (211,3) .OATAPT (211,4)  , 
10ATAPT(211,5)/1012.8,  16.3,  .76,  .136,  582.0/ 

OATA  LAT(212),L0N(212),ELEV1(212)/  40.38,  8.17,  40./ 

OATA  OAT APT (212,1) .OATAPT (212, 2) , OAT APT ( 212,3) .OATAPT (212,4) , 
10ATAPT(212,5)/1010.5,  15.9,  .76,  .181,  529.0/ 
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0«T«  LAT(213).L0N(213>.ELEV1(213I/  39.19t  9.03«  18./ 

08T8  0ATAPT(213.1I «0AT*PTI213t2)tD«TAFT(213t3) tOATAPT(213.<»)t  | 

10ATAPT(213«9>/1013.0.  16.5.  .76.  .171.  1132.0/ 

DATA  LATI21%),LONI21(»l«ELEVi(21A)/  37.95,  12.30,  lA./  * 

OATA  0ATAPT(21A,ll,0ATAPT{214,2t,0ATAPT(21A,3>,0ATAPT(21A,A>,  i 

10ATAPTI21A, 51/1013. A,  17.2,  .73,  .093,  516.0/ 

OATA  LATf215) ,L0N(215),ELEV1(215)/  37.28,  15.03,  17./ 

OATA  OATAPT (215,11  ,OATAPT 1215, 21,0 ATAPT 1215,3) ,OATAPT (215,A) 

10ATAPT(215,5)/1013.2,  17.2,  .76,  .159,  6A3.0/ 

OATA  LAT(216),L0N(216),ELEV1(216)/  38.12,  15.33,  51./ 

OATA  0ATAPT(216,1) ,OATAPT(216,2),OATAPT(216,3> .OATAPT (216, A) 

10ATAPT(216,5)/1008.9,  17.9,  .70,  .092,  973.0/ 

OATA  LAT(217),L0N(217),ELEV1<217)/  35.51,  1A,29,  72./ 

OATA  0ATAPT(217,1) ,0ATAPT(217, 2) .OATAPT  (217,3) .OATAPT (217, A) 

10ATAPT(217,5)/1006.1,  19.7,  .61,  .220,  6A0.0/ 

OATA  LAT(210),L0N(219),ELEV1(218)/  A3. 31,  16.26,  128./ 

DATA  0ATAPT(218,1) ,OATAPT( 218,2) t OATAPT (218, 3) .OATAPT (218, A) 

10ATAPT(218,5)/  999.8,  16.1,  .60,  .125,  816.0/ 

OATA  LAT(219) ,L0N(219) ,ELEV1(219)/  A2.22,  19.15,  33./ 

OATA  OATAPT(219,1),OATAPT(219,2),OATAPT(219,3) .OATAPT  (219, A) 

10ATAPT(219,5)/1011.0,  l^.A,  .6A,  .165,1632.0/ 

OATA  LAT(220),LON(220),ELEV1(220)/  A1.57,  21.38,  233./ 

OATA  OATAPT (220,1) ,OATAPT(220.2),OATAPT (220, 3) .OATAPT (220, A) 

10ATAPT(220,5)/  989.0,  12. A,  .72,  .160,  5AA.0/ 

OATA  LAT(221),L0N(221),FLEV1(221)/  35.37,  -.36,  9A./ 

OATA  0ATAPT(221,1) ,0ATAPT(221,2),0ATAPT(221,3) .OATAPT (221, A) 

10ATAPT(221,5)/1005.2,  17.2,  .77,  .239,  39A.0/ 

OATA  LAT(222),L0N(222),ELEV1(222)/  36. A3,  3.15,  25./ 

OATA  OATAPT (222,1), OATAPT (222. 2) .OATAPT (222,3) .OATAPT (222, A) 

10ATAPT(222,5)/1013.3,  17.3,  .77,  .310,  691.0/ 

OATA  LAT(223) ,L0N(223),ELEV1(223)/  36.50,  10. lA,  A./ 

OATA  OATAPT (223, 1) , OATAPT (223,2) .OATAPT (223,3) .OATAPT (223, A) 

10ATAPT(223,5)/101A.A,  18.3,  .72,  .195,  A70.0/ 

OATA  LAT(22A) ,L0N(22A),ELEV1(22A)/  35.30,  35.A8,  9./ 

OATA  0ATAPT(22A,1) ,0ATAPT(22A,2),0ATAPT(22A,3) .OATAPT (22A, A) 

10ATAPT(22A,5)/1010.1,  18.2,  .65,  .335,  71A.0/ 

OATA  LAT(225) •L0N(225),fLEVl(225)/  3A.27,  35. A8,  10./ 

OATA  0ATAPT(225,1) ,0ATAPT(225,2) ,0ATAPT(225,3) .OATAPT (225, A) 

10ATAPT(225,5>/1011.9,  19.1,  .67,  .295,  7A5,0/ 

OATA  LAT(226),L0N(226),ELEV1(226)/  33. A9,  35.29,  16./ 

OATA  OATAPT (226, 1 ) ,OATAPT( 226,2), OATAPT (226,3) .OATAPT (226, A) 

10ATAPT(226,5)/1011.2,  20.2,  .68,  .279,  517.0/ 

OATA  LAT(227),L0N(227),ELEV1(227)/  32.00,  3A.5A,  A9./ 

OATA  OATAPT (227,1) .OATAPT (227,2), OATAPT (227, 3) .OATAPT (227, A) 

10ATAPT(227,5)/10C7.7,  19.3,  .70,  .269,  529,0/ 

OATA  LAT(220) ,LON(228) ,ELEV1(228)/  31. A7,  35.13,  809./ 

OATA  0ATAPT(228,1) ,OATAPT(220,2) .OATAPT (229,3) .OATAPT  (228, A) 

10ATAPT(22e,5)/  921.1,  16.9,  .56,  .093,  A92.0/ 

OATA  LAT(229),L0N(229),ELEV1(229)/  A6.39,  1A.21,  A52./ 

OATA  0ATAPT(229,1) ,0ATAPT(229,2) .OATAPT (229,3) .OATAPT (229 , A) 

10ATAPT(229,5)/  963. A,  7,7,  .91,  ,197,  925.0/ 

OATA  LAT(230),LON(237),ELEV1 (230)/  A7,03,  12.57,3107./ 

OATA  OATAPT(230,1) ,OATAPT(23C,2) ,OATAPT(230,3) .OATAPT (230, A) 

10ATAPT(230,5)  / 692,7,  -6.0,  .8*»,  . 090,1A95.0/ 

OATA  LAT(231) ,L0N(231) .ELEVl (231)/  A6.02,  13.11,  92./ 

OATA  OATAPT (231,1) .OATAPT (231, 2) .OATAPT (231  ,3) .OATAPT  (231 , A) 

10ATAPT(231,5)/1005.0,  13.2,  .68,  .15A,1AA1.0/ 

OATA  LAT(232)  ,L0N(232)  ,ELEV1(232)  / A5.<.9,  15.59,  163./ 

OATA  OATAPT (232,1) .OATAPT (232, 2) .OATAPT (232,3) .OATAPT (2 32, A) 

10ATAPT(232,5)/  997.3,  11.6,  .72,  .163,  96A.0/ 
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OftTA  LAT(?33> ,LON(231) ,FLEV1(?33) / 43.52,  18.26,  637./ 

OAT  a DATAPT (233,1 ) .OATAPT 1233,2 ) .OATAPT  <233,31  , OAT APT (233,4) , 
10ATAPT(233,5)/  941.7,  9.8,  .72,  .149,  889.0/ 

OATA  LAT(234),L0N(234I,FL'^V1(234)/  43.12,  27.55,  41,/ 

OAT  A OAT  APT  (2  34,1)  , 0 ATAPT  ( 234 , 2)  , D AT  APT  (234 ,3)  , OAT  APT  (2  34 , 4)  , 
10ATAPT(234,5)/1011.6,  10.9,  .75,  .123,  475.0/ 

OATA  LAT(235),L0N(235),rLEWl(235)/  45.09,  29.40,  5,/ 

OATA  0ATAPT(2  35,1) , OATAPT ( 235 , 2 ) , 0 AT AP T ( 235 , 3 ) , OAT  APT (235,4), 
10ATAPT(235,5) /1015.7,  11.0,  .78,  .110,  352.0/ 

OATA  LAT(236),L0N(236),FLEV1(236)/  44.13,  28,38,  32./ 

OATA  OAT APT (2 36,1) , OAT APT ( 236 ,2 ), OAT APT (236 ,3 1 , OAT APT (236,4) , 
10AT4PT(236,5)/1012.4,  10.7,  .76,  .090,  375.0/ 

OATA  LAT(237) ,LON(237) ,FLEW1(237) / 49.04,  20,15,  709./ 

OATA  OAT APT (2 37,1) , OAT APT ( 237 ,2 ), OATAPT (237,3) , OAT  APT (237 ,4 ) , 
IOATAPT  (237,5)/  931.8,  *^.7,  .77,  .153,  608.0/ 

OATA  LAT  (238)  ,LON  (238)  ,'^LFVl  (238)  / 48.08,  20.48,  120./ 

OATA  OATAPT (2 38 , 1 ) , 0 AT APT ( 238 ,2 ) , D AT APT ( 238 , 3) , OA T A PT ( 238 , 4)  , 
10ATAPT(23e,5) /1002. 0,  9.7,  .77,  .165,  600.0/ 

OATA  LAT(239) ,L0N(239) ,FUEV1 (239) / 47.71,  19.02,  130./ 

OATA  OATAPT (2 39, 1 ) , OAT APT ( 239 ,2 ) .OATAPT ( 239, 3 ) , OA T APT (239,4)  , 
10ATAPT(239,5)/1000.6,  11.2,  .69,  .180,  630.0/ 

OATA  LAT(240)  ,LON(240),«^LEV1  (240)/  47.27,  21.38,  114./ 

OATA  OATAPT(240 ,1) , OA TAPT ( 240 , 2 ) , 0 AT A PT ( 24 0 , 3 ) .OATAPT (240,4)  , 
IOATAPT (240, 5)/lC02. 8,  10.3,  .75,  .160,  558.0/ 

OATA  LAT(241) ,L0N(241) ,FLFV1 (241)/  46.00,  18.14,  202./ 

OATA  DAT APT (241,1) , 0 AT  APT ( 241 , 2) , 0 AT APT ( 24 1 , 3) , OATAPT ( 241 , 4) , 
10ATAPT(241,5)/  992.3,  11.5,  .71,  .162,  683.0/ 

OATA  LAT(242),L0N(242),FLEV1(242)/  46.15,  20.06,  84./ 

OATA  OATAPT (242,1) .OATAPT (242,2) , 0 AT APT ( 242 , 3 ) , OATAPT (242 ,4 ) , 
10ATAPT(242,5) /1006. 1,  11.5,  .73,  .173,  559.0/ 

OATA  LAT(243) ,LON (243) ,ELEV1 (243) / 44.48,  20.28,  132,/ 

OATA  OATAPT (243,1) , 0 A T APT ( 243, 2 ) , DAT  APT ( 243 , 3) .OATAPT (243,4) , 
104TAPT(243,5)/10C1. 0,  11.8,  .70,  .155,  659.0/ 

OATA  LAT(244) ,L0N(244) ,5LEV1 (244) / 43.49,  23,14,  33./ 

OATA  OAT APT (244,1) , OA TAFT ( 244 , 2 ) , OAT  APT ( 244 ,3 ) , OATAPT (244 ,4 ) , 
IOATAPT (244, 5)/1012. 8,  11.7,  .74,  .158,  555.0/ 

OATA  LAT(245) ,LON(245) .FLFVl (245)/  42.49,  23.23,  588./ 

OATA  OAT AFT (245,1)  , 0 A T APT ( 245 , 2) , 0 AT APT ( 245 , 3) .OATAPT (245,4) , 
10ATAPT(245,5)  / 950.7,  10.4,  ,72,  . 149,  622. C/ 

OATA  LAT  (246) ,LCN (246) ,FLPV1 (246) / 47.10,  27,36,  103./ 

OATA  0ATaPT(246,l) .OATAPT (246,2) .OATAPT (246,3) .OATAPT (246,4) , 
10ATaPT(246,5) /1005.7,  9.4,  .71,  .125,  506.0/ 

OATA  LAT(247) , LON (247) .FLFVl (247)/  45.46,  21.15,  91./ 

OATA  OATAPT (247,1) .OATAPT (247, 2), OAT  APT (24 7, 3), OATAPT (247, 4), 
10aTaPT(247,5) /1005. 3,  10.9,  .73,  .175,  625.0/ 

^ATA  LAT(248)  ,L0N(2‘*8)  .FLFVl  (248)  / 45.48,  24.  OP,  452./ 

OATA  DATAPT(248,1  ) , OATAPT ( 248 , 2 ) , OAT APT ( 24 8 , 3 ) .OAT  APT ( 248 , 4)  , 
10aTaPT(24e,5)/  P60.5,  7.7,  ,74,  .155,  643.0/ 

OATA  LAT(24g)  ,L0N(249)  , '■LEVI  (249)/  44.2'",  26.06,  82./ 

TATA  OATAPT (24Q,  1)  .OAT SFT( 240,2)  , 0 AT APT ( 24 9 , 3)  , OAT  AFT ( 24  9 , 4)  , 

IOATAPT (249,5) /1007. 0,  11.1,  .70,  ,095,  579.0/ 

.function  05C1ML  CONV'-'^T':  OFGPEF' -MINUTES  (DO. MM)  INTO  DECIMAL 
.O'^GPEES. 

OFCIML (X)=FLOAT (INT ( X ) ) ♦ ( X-FLT A T ( I NT ( X ) ) ) / , 60 
. TNITIALIZATICN. 

OATA  N,NOF,NFPM, NIP/249, 20, 15,1/ 
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NCP»5 

CONTUR«35.0 


.....ERROR  CHECK  ON  INTERPOLATION  POINT. 

IF  (XL  AT. GT. SZ.  <•.  OR.  XL  AT.  LT.iiO.O.  OR.  XLON.GT.  25.0.  OR.XLON.lt.  1.2) 
1MRITE(6,1000) 

CONVERT  POSITIONS  OF  DATA  POINTS  TO  OECINAL  OEGREES  ANO  REDUCE 

PRESSURE  TO  SEA  LEVEL  FOR  HOPE  ACCURATE  INTERPOLATION. 

00  100  Isl.N 

□ECLAT(I)xOECIHL(LAT(I)) 

OECLON(I)sOECIML(LON(II) 

PRS(I)xOATAPTII,l)*EXPiq.8*ELEVl(I)/(287.0«**(273.16*OATAPT (1,2) ) ) ) 
100  CONTINUE 

. . .CONVERT  POSITION  OF  INTFRPOLATION  POINT  TO  OECINAL  OEGREES. 

OCXLAT«OECIHL(XLAT) 

OCXLONsOECIHL(XLON) 

. . .ROUND  ELEVATION  TO  CLOSEST  CONTOUR  INTERVAL. 

El=ELEV/CONTUR 
E2=AINT (El) 

E3xAINT((El-E2)/.5) 

ELV*(E2*E3)PC0NTUR 

FIND  THE  20  DATA  STATIONS  CLOSEST  TO  THE  INTERPOLATION  POINT  ANO 

STORE  THEIR  SUBSCRIPTS  IN  ARRAY  ISUB. 

00  300  Jxl.NOP 

CLSDIS«1.CE«6 

JJ=J-1 

DO  300  Kxl.N 

TRYOIS*(OCXLAT-OECLAT(K) )**2* (OCXLCh-OECLON (K)  ) **2 
IF(TRYOIS-CLSOIS)i0,10,300 
10  IF(JJ)30,30,20 
20  00  200  LxltJJ 

IF(ISUB(L)-K)  200,300,200 
200  CONTINUE 
30  CLSOISxTRYOIS 
ISUB(J)>K 
300  CONTINUE 

POT  HETEOROLOGICAL  OATA  FRON  THE  20  STATIONS  INTO  ARRAYS  TO  SEND 

.....TO  lOBVIP. 

00  AOO  Kxl.NCP 
KKxISUB(K) 

PRESS(K)aPRS(KK) 

TEMP(K)sOATAPT(KK,2) 

RELHUM(K)sOATAPT (KK,3) 

3ETA(K)«0ATAPT(KK,4) 

H»1(K)«OATAPT(KK,5) 

LATT(K)xOECLAT(KK) 

LONG(K) «OECLON(KK) 

tOBVIP  COES  THE  INTERPOLATIONS  (AKIHA, 1975) . 
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CALL  IDBVIP(1«N0P,LATT,L0NC« PRESS, NIPtOCXL AT. OCXLON.PRES. INK. MK) 
CALL  I0BVIPJ3,N0P.LATT,L0NG.TEHP,NIP,0CXLAT,0CXL0N.TMP,1NK,NKI 
CALL  IOBVIPO.NOP.LATT.LONG.RELHUN.NIP.OCXLAT.OCXLON.RLHUN.INK.WK) 
CALL  lOBVIPO. NOP. L ATT, L0N6, BETA, NIP. OCXLAT.OCXLON, BET. INK. NK) 

CALL  lOBVIPd.NOPH.LATT.LONG.NN.NIP.OCXLAT.OCXLON.H.INK.NK) 

C 

EVALUATE  STATION  PRESSURE  PROM  PSEUOO-SEA-LEVEL  VALUE  AT  THE 

C INTERPOLATION  POINT. 

C 

PRES«PRES*EXP(-9.A*£LV/«2Br.0‘.*(273.16*TMP» 

RETURN 

1000  FORHATI12HO*  NARNING  •/70H  INTERPOLATION  PT  OUT  OF  DESIGNATED  BOUN 
IDS  (SEE  INITIAL  COMMENT  CAROS),/) 

END 


61 


oooooor>oonooonooor>nooonooooooooooooooo 


SUBROUTINE  IDBVIP ( HO ,NOP.XO « VO«ZO .NIP, XI , YI«ZI .IMK«NK) 
THIS  SUBROUTINE  PERFORHS  BIVARIATE  INTERPOLATION  WHEN  THE  PRO- 


NO 


NOP 

XO 


JECTIONS  OF  THE  DATA  POINTS  IN  THE  X-Y  PLANE  ARE  IRREGULARLY 
DISTRIBUTED  IN  THE  PLANE. 

THE  INPUT  PARANETERS  ARE 

NODE  OF  CONPUTATION  tHUST  BE  1,  2.  OR  3). 

1 FOR  NEN  XO-YO, 

2 FOR  OLD  XD-YO,  NEN  XI-YI, 

3 FOR  OLD  XO-YO,  OLD  XI-YI, 

NUMBER  OF  DATA  POINTS  (MUST  BE  k OR  GREATER), 

ARRAY  OF  DIMENSION  NOP  CONTAINING  THE  X 
COORDINATES  OF  THE  DATA  POINTS, 

ARRAY  OF  DIMENSION  NOP  CONTAINING  THE  Y 
COORDINATES  OF  THE  DATA  POINTS, 

ARRAY  OF  DIMENSION  NDP  CONTAINING  THE  Z 
COORDINATES  OF  THE  DATA  POINTS, 

NUMBER  OF  OUTPUT  POINTS  AT  MHICH  INTERPOLATION 
IS  TO  BE  PERFORMED  (MUST  BE  1 OR  GREATER), 

ARRAY  OF  DIMENSION  NIP  CONTAINING  THE  X 
COORDINATES  OF  THE  OUTPUT  POINTS, 

ARRAY  OF  DIMENSION  NIP  CONTAINING  THE  Y 
COORDINATES  OF  THE  OUTPUT  POINTS. 

THE  OUTPUT  PARAMETER  IS 


YO 


ZO 


NIP 


XI  = 


YI 


ZI  = 


INTERPOLATED  Z 


THE 


OTHER 
INK 


ARRAY  OF  DIMENSION  NIP  WHERE 
VALUES  ARE  TO  BE  STORED. 

PARAMETERS  ARE 
INTEGER  ARRAY  OF  DIMENSION 
MAXO<3l,27»NCP)*NOP*NIP 
TO  BE  USED  AS  A WORK  AREA, 

ARRAY  OF  DIMENSION  e*NOP  TO  BE  USED  AS  A 
WORK  AREA, 

WHERE  NCP  IS  THE  NUMBER  OF  ADDITIONAL 


WK  = 


DATA  POINTS  USED  FOR 
ESTIMATING  PARTIAL  DERIVATIVES  AT  EACH  DATA  POINT.  THE  VALUE 
OF  NCP  MUST  BE  GIVEN  THROUGH  THE  IDCM  COMMON.  NCP  MUST  BE  2 
OR  GREATER,  BUT  SMALLER  THAN  NOP. 

THE  LUN  CONSTANT  IN  THE  DATA  INITIALIZATION  STATEMENT  IS  THE 
LOGICAL  UNIT  NUMBER  OF  THE  STANDARD  OUTPUT  UNIT  AMO  IS, 
THEREFORE,  SYSTEM  DEPENDENT. 

DECLARATION  STATEMENTS 

DIMENSION  XO(IOO) ,YO( 100) ,ZO( 100) , XI (1000) ,YI( 1000) , 

1 ZI(IOOO)  ,IWK(i»100)  ,NK(800) 

COMMON/IOCN/NCP 


COMMON/IOLC/NIT 
COMMON/IDPI/NCF 
DATA  LUN/6/ 

SETTING  OF  SOME  INPUT 
(FOR  MD>1,2,3) 

10  MOO»MO 
NOP0*NDP 
NIP0«NIP 
NCP0*NCP 

ERROR  CHECK.  (FOR  M0=l,2,3) 

20  IF(MOO.LT.1.0R.HOO.GT.3) 
IF(NDP0.LT.i.) 

IF(NIPO.LT.l) 

tF(NCP0.LT.2.OR.NCP0.GE,NOP0) 

GD  TO  21 


PARAMETERS  TO  LOCAL  VARIABLES. 


GO 

GO 

GO 

GO 


TO  90 
TO  90 
TO  90 
TO  90 
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21 

NOPPV=IMK(l) 

NCPPV*1NK(2) 

IF(NOPO.NE.NOPPV) 

GO 

TO 

90 

IF(NCP0.NE.NCPPV1 

GO 

TO 

90 

22 

IF(HOO.GE.3) 

GO 

TO 

23 

IMK(3I«NIP 

GO  TO  30 

23 

NIPPV=INK(3) 

IF(NIP0.NE.NIPPV> 

GO 

TO 

90 

C ALLOCATIOH  OF  STORAGE  AREAS  IN  THE  IMK  ARRAY.  <FOR  MD*1»2*3» 
30  JN1PT*16 

JMIHL36«NOPO«1 
JNIMK>JMINL 
JNIPL»2A*NOPO*1 
JMIMP«30*NOP0^1 
JH1PC»27*NOP0^1 
JMITOsHAXO (31.27»NCPQ ) •NOPO 
C TRIANGULATES  THE  X-Y  PLANE.  JFOP  HO=l) 
itO  IF(HOO.GT.l)  GO  TO  SO 

CALL  IOTANG(NOPO,XO.YO,NT,IMICIJMIPT) .NL.IHKIJMIPL)  . 

1 IMKIJMIML)  .IHKUNIMP)  ,MK) 

IMK(5)*NT 

IMK(6)«NL 

IMK(S|sNT 

IMKI6)«NL 

IF(NT.rO.O)  RETURN 

C DETERMINES  NCP  POINTS  CLOSEST  EACH  DATA  FJINT.  (FOR  ND=1I 
SO  IF(HOO.GT.l)  GO  TO  60 

CALL  IOCLOP(NOPO«XOtYO«NCPO«IMK(JHIPC) ) 

IF(IMK(JMIPC)  .EQ.O)  RETURN 

C LOCATES  ALL  POINTS  AT  WHICH  INTERPOLATION  IS  TO  BE  PERFORMED. 
C (FOR  MO»l«2) 

60  IF(HQ0.EQ.3)  go  TO  70 
NITaO 

JMITsJMITO 
DO  61  IIP»1.NIP0 
JMIT«JI(ITM 

CALL  I0LCTN(NCP(!,X0,Y0,NT,IMK(JHIPT)  ,NL.IMK(JMIPL)  . 

1 XI  (I IP  I .YKIIPI  tIHKIJMITI  .INK(JHIHK)  «MK) 

61  CONTINUE 

C ESTIMATES  PARTIAL  DERIVATIVES  AT  ALL  DATA  POINTS. 

C (FOR  MOsl«2t3) 

70  CALL  IOPORV(NOFOtXD,YD«ZOtNCPO«lMK(JMIPC),HK) 

C INTERPOLATES  THE  ZI  VALUES.  (FOR  MD=1,2,3I 
so  NCF*0 

JMIT=JNITC 
00  81  1IP=1,NIP0 
JNIT*JNIT*1 

CALL  TOPTIP(XO«Vn.Zr),NT,IMK(JMIPT)  ,NL.IMK(JMIPL)  .MK, 

1 IMK(JMIT)  ,XI(I1P)  .YKIIP)  .ZKIIP)) 

81  CONTINUE 
RETURN 
C ERROR  EXIT 

90  WRITE  (LUN,2C9r»  MOO » NOPO.NIPO , NCP C 
RETURN 

C FORMAT  statement  FOR  ERROR  MESSAGE 
2090  F0RMAT(1X/41H  •••  IMPROPER  INPUT  PARAMETER  VALUE (S>./ 

1 7H  mo  =, I4,10X tSHNOP  = , 16 ,1 CX ,5 HNIP  =,I6, 

2 lOX.SHNCP  =,16/ 

3 35H  ERROR  OFTECTEO  IN  ROUTINE  IDBVIP/) 

END 
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SUBROUTINE  IDCLOP(NOP,XO*yO.NCP«IPC) 

C THIS  SUBROUTINE  SELECTS  SEVERAL  OATA  POINTS  THAT  ARE  CLOSEST 
C TO  EACH  OF  THE  OATA  POINT. 

C THE  INPUT  PARANETERS  ARE 
C NOP  * NUMBER  OF  OATA  POINTS* 

C XO.VO  B ARRAYS  CONTAINING  THE  X ANO  Y COOROINATES 

C OF  OATA  POINTS, 

C NCP  • NUMBER  OF  DATA  POINTS  CLOSEST  TO  EACH  DATA 

C POINTS. 

C THE  OUTPUT  PARAMETER  IS 

C IPC  B ARRAY  OF  DIMENSION  NCP*NOP,  WHERE  THE  POINT 

C NUMBERS  OF  NCP  OATA  POINTS  CLOSEST  TO  EACH  OF 

C THE  NOP  OATA  POINTS  ARE  TO  BE  STORED, 

C THIS  SUBROUTINE  ARBITRARILY  SETS  A RESTRICTION  THAT  NCP  MUST 
C NOT  EXCEED  25, 

C DECLARATION  STATEMENTS 

DIMENSION  XO(IOO)  ,YOI100),IPC(<»0  0> 

DIMENSION  OSOO (25) ,IPC0(25) 

OATA  NCPMX/25/,  LUN/5/ 

C STATEMENT  FUNCTION 

0SQF(X1,Y1,X2,V2)b(X2-X1)**24(Y2-V1)**2 
C PRELIMINARY  PROCESSING 
10  NOPObNOP 
NCPObNCP 

IF(NDP0,LT.2I  go  to  90 

IFINCPO.LT.l.OR.NCPO.GT.NCPMX.OR.NCPO.GE.NOPO)  GO  TO  90 
C CALCULATION 

20  no  29  IP1=1,NDP0 

X1=X0(IP1) 

ri*ro(iPi) 

JIbO 

OSQMXbO.O 

00  22  IP2=1,NDP0 

IF(IP2,EQ.IP1)  GO  TO  22 
DSQI=0S0F(Xl,Yl,Xn(IP2) ,yO(IP2)) 

OSQO(J1)=OSQI 

IPC0(J1)=IP2 

IFIPSQI.LE.OSQMX)  GO  TO  21 

OSOMXsOSQI 

JMXsJl 

21  IF( Jl.GE.NCPO)  GO  TO  23 

22  CONTINUE 

23  IP2NNbIP2»1 

IF(IP2MN.CT,NOPO)  GO  TO  26 

00  25  IP2=1P2MN,NOPO 

IF(IP2.EQ,IP1)  GO  TO  2S 
0SQI*0SQF(Xl,Yl,Xn(IP2) ,Y0(IP2)) 

IFiOSQI.GE.OSQHXI  GO  TO  25 

OSQO(JMX)bOSOI 

IPC0IJMX)bIP2 

OSQMXbO.O 

00  24  J1=1,NCP0 

IF(OSQO(J1) .LE.OSQMX)  GO  TO  24 
OSQMXsOSQO ( Jl) 

JMXbJI 

24  CONTINUE 

25  CONTINUE 

26  J1*)IP1-1)*NCP0 

00  28  J2*1,NCP0 


,1 

1? 
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IPC(JllsIPC0(J2l 
Za  CONTINUF 
Z9  CONTINUF 
PFTURN 
C ERROR  EXIT 

90  NRITE  (LUN,Z090)  NOPO.NCPO 
IPCIDaO 
RETURN 

f RN»T  STATEMENT 

£U,J  FORNATIIX/AIH  •••  IMPROPER  INPUT  PARAMETER  VALUEtS)./ 

1 8H  NOP  «,I5.5X,9HNCP  =,15/ 

2 35H  ERROR  DETECTED  IN  ROUTINE  lOCLOP/) 

END 
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SUBROUTINE  lOCRIO IXD, Y0«NT« IPT .NL *IPL .NXI .NYI , XI , YI, 

1 NGP.IGP) 

THIS  SUBROUTINE  ORGANIZES  GRID  POINTS  FOR  SURFACE  FITTING  BY 
SORTING  THEN  IN  ASCENDING  ORDER  OF  TRIANGLE  NUHBERS  AND  OF  THE 
BORDER  LINE  SEGMENT  NUMBER. 

THE  INPUT  PARAMETERS  ARE 

XOtVO  s ARRAYS  CONTAINING  THE  X ANC  Y COORDINATES  OF 
DATA  POINTS, 

NT  « NUMBER  OF  TRIANGLES, 

IPT  * ARRAY  CONTAINING  THE  POINT  NUMBERS  OF  THE 
VERTEXES  OF  THE  TRIANGLES, 

NL  « NUMBER  OF  BORDER  LINE  SEGMENTS, 

IPL  » ARRAY  CONTAINING  THE  POINT  NUMBERS  OF  THE  END 
POINTS  OF  THE  BORDER  LINE  SEGMENTS  AND  THEIR 
RESPECTIVE  TRIANGLE  NUMBERS, 

NXI  « NUMBER  OF  GRID  POINTS  IN  THE  X COORDINATE, 

NYI  » NUMBER  OF  GRID  POINTS  IK  THE  Y COORDINATE, 

XI, YI  = ARRAYS  CONTAINING  THE  X AND  Y COORDINATES  OF 
THE  GRID  POINTS. 

THE  OUTPUT  PARAMETERS  ARE 

NGP  = INTEGER  ARRAY  OF  DIMENSION  2*(NTt?*NL»  WHERE  THE 
NUMBER  OF  GRID  POINTS  THAT  BELONG  TO  EACH  OF  THE 
TRIANGLES  OR  OF  THE  BORDFP  LINE  SEGMENTS  ARE  TO 
BE  STORED, 

IGP  = INTEGER  ARRAY  OF  DIMENSION  NXI*NYI  WHERE  THE 

GRID  POINT  NUMBERS  ARE  TO  BE  STORED  IN  ASCENDING 
ORDER  OF  THE  TRIANGLE  NUMBER  ANO  THE  BORDER  LINE 
SEGMENT  NUMBER. 

DECLARATION  STATEMENTS 

DIMENSION  XOCIOOJ , YO ( 1001 , IPT (565 J , IPL <300 » , 

1 XI(lOl) ,YI(101) ,NGP(8aO) ,IGP(10Z01) 

C PRELIMINARY  PROCESSING 
10  NT0»NT 
NL0=NL 
NXIOsNXI 
NYIOsNYI 
NXINYI«NXIO*NYIO 
XIMNzAMINKXI  (1 ) , X I I NXI  0 ) ) 

XIMXsAMAXKXKl)  ,XI(NXI0)) 

YIMNsAMINKYKl)  ,YI(NYI0»> 

YIMX>AMAX1(YI(1I  ,YI (NYIO)) 

C DETERMINES  GRID  POINTS  INSIDE  THE  DATA  AREA. 

?0  JNGP0»C 

JNGP1=2*(NT042*NL0) ♦! 

JIGP0»0 

JIGP1*NXINYI41 
DO  39  IT0=l,fJT0 
NGPOsO 
NGPIbO 
IT0T3=IT0*3 
IPl»IPT(IT0T3-2) 

IP2*IFT(IT0T3-1) 

IP3*IFT(IT0T3» 

Xl=XO  (IPl) 

YlsYDdPl) 

X2«X0(IP2J 

Y?»YO(IP2» 

X3=XO  (IP3I 
Y3*YO(IP3» 

XMN3AMIN1(X1,X2,X3I 
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XHXs«HAXl(Xl.X2<X3) 

YMNsAMINl(Yl«VZ«Y3) 

YMX*AMAXt(YliY2.Y3) 

INSO>0 

00  22  IXI*1«NXI0 

IFtXKIXD.GE.XHN.AND.XKIXI)  .LE.XHX)  GO  TO  21 
IFdNSO.EQ.O)  GO  TO  22 

1XIMX*IXT-1 
GO  TO  23 

21  IFdNSO.EQ.ll  GO  TO  22 
INS0«1 

IXiHNsIXI 

22  CONTINUE 

IFdNSO.EQ.O)  GO  TO  38 

IXINXxNXIO 

23  0X21>X2-X1 
0X32«X3-X2 
0X13*X1-X3 
0Y21*Y2-Y1 
OY32*Y3-Y2 
0Y13«Y1-Y3 

DO  37  lYIsl.NYIO 
YIIxYHIYI) 

IFIYII.LT.YMN.OO.VII.GT.YMX)  GO  TO  37 

Al«  <YII-Y1)*0X21 
A2=  JYII-Y2)*0X32 
A3sIYlI-Y3)*OX13 
DO  36  IXIsIXIHN.IXTMX 
XIIxXKIXI) 


LxO 

IF«A1-(XII-X1)*0Y21)  36,25,26 

25  L = 1 

26  IF(A2-<XII-X2)*DY32)  36,27,28 

27  L = 1 

28  IF(A3-CXII-X3)*0Yn)  36,29,30 

29  L = 1 

30  IZT=NXIO*dYI-l)*IXI 
IF(L.EQ.l)  GO  TO  31 
NGP0=NGP0*1 
JIGP0»JIGP0^1 
IGP)JIGPO)«IZI 

GO  TO  36 

31  IFtJIGPl.GT.NXINYI)  GO  TO  33 

DC  32  JIGP1I=JIGP1,NXINYI 

IFCIZI.EQ.1G“(JIGP1I))  GO  TO  36 

32  CONTINUE 

33  NGP1=NGP1»1 
JIGP1=JIGP1-1 
IGP<JIGPl)sIZI 

36  CONTINUE 

37  CONTINUE 

38  JNGPOxJNGPO^l 
NGPCJNGPOIxNGPO 
JNGP1=JN6P1-1 
NGP  (JNGP1)=NGP1 

39  CONTINUE 

C 0FTE9MINES  GRID  POINTS  0UTS1«F  THE  DATA  AREA. 

90  00  79  ILOxl,NLO 
NGP0>0 
NGP1«C 
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41 

42 

43 


45 

46 

47 
45 

49 

50 


51 


52 

53 


56 

57 
55 


XMNsANINl (XI. X2) 
XNX=«HAX1(X1.X2> 
YHNsAMIMl (Y1,Y2> 
YHXsAMAXl (Y1,Y2» 


IL0T3*IL0*3 
lPl«IPL(lL0T3-2» 

IP2«1PL(IL0T3-1) 

X1>XDIIP1) 

Y1*Y0(IP1> 

X2sX0(IP2) 

Y2*yO(IP2l 
XMN«XIHN 
XHXsXIHX 
YMNsYlHN 
YMX»YIHX 
1F(Y2.GE.Y1> 

IFIY2.LE.Y1) 

IFCXZ.LE.Xl) 

IFtX2.GE.Xl> 

INSO*0 

00  42  1X1=1. NXIO 

IFIXKIXD.GF.XHN.ANO.XI  (IXI)  .LE.XHX)  GO  TO  41 
IFCINSO.EQ.OI  GO  TO  42 

IXIMXslXI-1 
GO  TO  43 

IFIlNSOiEQ.U  GO  TO  42 

INSOsl 
IX1HN«IX1 
CONTINUE 

IF(INSO.EQ.O)  GO  TO  55 

IXINXxNXlO 
0X21*X2-X1 
0Y21=Y2-Y1 
DO  57  IYI=1.NY10 
YI1*YI(IYI) 

IFIYII.LT.YMN.OP.YII.GT.YMX)  GO  TO  57 

A1*CY1I-Y1)*0X21 
Bls(YII-Yl)*0Y21 
B2«(YII-Y2)*OY21 
DO  56  IXI=IXIMN.1X1NX 
XII*XI(IXI) 

L»0 

IF(A1-(XII-X1»*0Y21)  46.45.56 


L = 1 

IF(0X21*(XII-X1»»BH  56.47.48 

L*1 

IF(0X21*(XII-X2»»82)  50.49.56 

L»1 

IZlsNXIO*(lYI-l)*IXI 
IFd.EQ.ll  GO  TO  51 
NGP0«NGP0*1 
JIGP0=JIGP0*1 
IGP(JIGPO)«IZI 
GO  TO  56 

IF(JI6P1.GT.NXINYI»  GO  TO  53 

00  52  JIGP1I=JIGP1.NXINYI 

IFdZI.EQ.lGP  (JIGPII)  ) GO  TO 

CONTINUE 
NGP1*NGP1*1 
JIGP1=J1GP1-1 
1GP(JIGP1)=IZI 
CONTINUE 
CONTINUE 
JNGPO= JNGP0»1 
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*NGPIJNGPO)sNGPO 
JNGP1«JNGP1-1 
NGPIJNGPD^NGPl 
(0  NGP0*0 
NGP1«0 

ILP1«MOOIILO.NLO)»1 

ILP1T3*ILP1*3 

IP3*IPLIILPtT3-l) 

X3>XD(IP3I 

Y3*Y0IIP3) 

XNNsXIMN 

XHXsXIHX 

YMN*YIHN 

YMX»YIMX 

IF(Y3.GE.Y2,»N0.Y2.GE.Y1>  XMN=X2 
IF(Y3«LE.Y2*«N0.Y2.LE.Y1)  XMX«X2 
IFIX3.LF.X2.XNO.X2.LE.X1)  YMN=Y2 
IF(X3.GC.X2.AN0.X2.GE.X1)  YHXsY2 
INSO«0 

00  62  1X1=1, NXIO 

1F(XHIX1».GE.XMN.*ND.XI(IXU  .LE.XMX)  GO  TO  Gi 
IFIINSO.EQ.O)  GO  TO  62 

1XIMX*IXI-1 
GO  TO  63 

61  IFIlNSO.FQ.l)  GO  TO  62 
INS0=1 

IXIHNsIXI 

62  CONTINUE 

IFIINSO.EQ.O)  GO  TO  76 

IXINX=NXI0 

63  0X21>X2-X1 
0Y21*Y2-Y1 
0X32*X3-X2 
0Y32=Y3-Y2 

00  77  IYI*1,NYI0 
YIIsYIlIYI) 

IFITII.LT.YHN.OR.YII.GT.YNX)  GO  TO  77 

B2*IYII-Y2)*0Y21 
B3*  IYII-Y2)*0Y32 
00  76  IXI=IXIMN,IXINX 
XII=XI(IXI) 

L = 0 

IFrDX21*IXII-X2)4P2)  76,65,66 

65  L = 1 

66  IFI0X32*(XII-X2)*B3)  7C,67,76 

67  L = 1 

70  I2I=NXI0*(IYI-1)>IXI 
IF(L.EQ.l)  GO  TO  71 
NGPe=NGPO*l 
JIGP0=JIGP0*1 
IGP|JIGP0)*12I 

GO  TO  76 

71  IFIJIGPl.GT.NXINYI)  CO  TO  73 

DC  72  JIGP1I=JIGP1,NXI‘ YI 

IFIIZI.EQ, IGP (JTGPII) ) GO  TO  76 

72  CONTINUE 

73  NGP1=NGP1*1 
JIGP1=JIGP1-1 
IGPUIGP1)=IZI 

76  CONTINUE 
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77  CONTINUE 

78  JNGPOsJNGPO^l 


NGPIJNGP0)*NGP0 
UNGPl^JNGPl'l 
NGP(JNGP1)=NGP1 
79  CONTINUE 
RETURN 
END 
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SUBROUTINf  It)LCTN(NOP,VOtYO«NTtIPT«NLt IPLtXII.YIItITI, 

1 1NK,MK) 

THIS  SUBROUTINC  LOCATES  A POINT*  I.E.*  DETERMINES  TO  MHAT  TRI> 
ANGLE  A GIVEN  POINT  (XII.YITI  BELONGS.  MHEN  THE  GIVEN  POINT 
DOES  NOT  LIE  INSIDE  THE  DATA  ARPA.  THIS  SUBROUTINE  DETERMINES 
THE  BORDER  LINE  SEGMENT  IN  THE  AREA  ABOVE  WHICH  THE  POINT 
LIES.  OR  TWO  BORDER  LINE  SEGMENTS  BETWEEN  TWO  AREAS  ABOVE 
WHICH  THE  PCINT  LIES, 

the  input  parameters  are 

NOP  « NUMBER  OF  DATA  POINTS, 

XD.YO  « ARRAYS  CONTAINING  THE  X AND  Y COORDINATES 
OF  DATA  POINTS, 

NT  « NUMBER  OF  TRIANGLES, 

IPT  * ARRAY  CONTAINING  THE  POINT  NUMBERS  OF  THE 
VERTEXES  OF  TH^^  TRIANGLES, 

NL  « NUMBER  OF  BORDER  LINE  SEGMENTS, 

IPL  « ARRAY  containing  ThF  POINT  NUMBERS  OF  THE  END 
POINTS  OF  THE  BORPEP  LINE  SEGMENTS  AND  THEIR 
RESPECTIVE  TRIANGLE  NUMBERS, 

XII, YII  = X AND  Y COORDINATES  OF  THE  POINT  TO  BE 
LOCATED. 

THE  OUTPUT  PARAMETER  IS 

ITI  » TRIANGLE  NUMBER,  WHEN  THE  POINT  IS  INSIDE  THE 
DATA  AREA,  OR 

TWO  BORDER  LINE  SEGMENT  NUMBERS,  ILl  AND  IL2, 

CODED  TO  ILIP (NT»NLI »IL2,  WHEN  THE  POINT  IS 
OUTSIDE  THE  DATA  AREA. 

THE  OTHER  PARAMETERS  ARE 

IWK  * INTEGER  ARRAY  OF  DIMENSION  lf*NOP  TO  BE  USED 
INTERNALLY  AS  A WOPK  AREA, 

WK  = ARRAY  OF  CINENSION  fl*NrP  TO  BE  USED  INTERNALLY 
AS  A WORK  AREA. 

DECLARATION  STATEMENTS 

DIMENSION  XOaOC)  ,YOtlOC»  ,IPT  (5  85)  , IPL  <3  3 0) , 

1 IWK(1803) ,MK(80Q) 

DIMENSION  NTSCIB) ,IOSC(9) 

COMHON/IOLC/NIT 

c statement  function 

ST0E(X1,Y1,X2,Y2,X3, Y3) = (Y3-Y1) ’(XS-Xl) -(X3-Xl)*tY2-Yl) 

0 PRELIMINARY  PROCESSING 
13  NOPO>NDP 
NT0*NT 
NLOsNL 
NTLsNTO+NLO 
X0»XII 
YOsYII 

C PROCESSING  FOP  A NEW  SET  OF  DATA  POINTS 
20  IF(NIT,NE,0)  GO  TO  30 
N1T  = 1 

0 - DIVIDES  THE  X-Y  PLANE  INTO  NINE  RECTANGULAR  SECTIONS. 
XMN=XO(l) 

XMXxXMN 

YMNxYDIl) 

YMXsYMN 

00  21  IOPs2,NrPO 
XISXO(IOP) 

YT»YD(IOP) 

IF(XI.LT.XMN)  XMNxXI 

IF(XI.GT.XMX)  XMXrXI 

IF(YI.LT,YMN)  ymn=yi 
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YI1X»YI 


IFIY1.6T.YHX) 

?1  CONTINUE 

XS1=  ( XMH4'XHN4XHX) /3 . 0 
XS2«  (XHN'fXHX^XHX)  /3. 0 
YSl* I YMN*YMN4YHX) /3 . 0 
YS2* ( YHHfYHX^YMX) /3. 0 

C - DETERMINES  AND  STORES  IN  THE  INK  ARRAY  TRIANGLE  NUMBERS  CF 
C - THE  TRIANGLES  ASSOCIATED  NITH  EACH  OF  THE  NINE  SECTIONS. 

DO  22  ISC=1,9 

NTSC(ISC)=0 
lOSCIISOsO 

22  CONTINUE 
IT0T3»0 
JNKsO 

00  27  IT0=1.NT0 

IT0T3»TT0T2*3 
Il*IPTIlT0T3-2) 

12»1PT(IT0T3-1) 

I3*IPTIIT9T3) 

XHNsAMINKXOlIl)  .X3II2)  .XO  (13)  > 

XHXsAMAXKXOCIl)  «X0(I2),X0(I3)) 

YMNxAMINKYOCll)  ,Y0(I2)  ,Y0)I3)  ) 

YMXxAMAXllYOdl)  .Y0(I2).Y0(I3)) 

IF(YMN.GT.YSl) 

IF(XMN.LE.XSl) 

IF(XHX.GE.XS1.AND.XHN.LE.XS2) 

IF(XHX.GE.XS2> 

23  1F( YMX.LT.YS1.0R. YHN.GT.YS2) 

IF(XMN.LE.XSl) 

IF (XHX.GE.XSl.  AN0.XMN.LE.XS2) 

IFIXHX.GE.XS2) 

2<*  IF(YMX.LT.YS2) 

IF(XMN.LE.XSl) 

IF(XMX.GE.XS1.AN0.XMN.LE.XS2) 

IFIXMX.GE.XS2) 

25  00  26  ISCxl«9 

IFCIOSCdSC)  .EO.O)  GO  TO  26 
JINKa9*NTSC<ISC(»ISC 
INK(JINK)>IT0 
NTSC(1SC)*NTSC(ISC)*1 
IOSC(TSC)=0 

26  CONTINUE 

C - STORES  IN  THE  HK  ARRAY  THE  MINIMUM  AND  MAXIMUM  OF  THE  X AND 
C - Y COORDINATE  VALUES  FOR  EACH  OF  THE  TRIANGLE. 

JNKsJNK«<» 

NKIJNK-3)xXMN 
NK( JNK-2) xXMX 
MK( JHK«l)xYMN 
NK(JMK)  xVHX 

27  continue 
GO  TO  60 

C CHFCKS  IF  IN  THE  SANE  TRIANGLE  AS  PREVIOUS. 

30  IT0*1TIPV 

tf(ito.gt.nto)  go  to  LO 

TT0T3*ITn*3 

TPl=IPT(IT0T3-2) 

X1=X0  (IFl) 

YlxYO  (IPl) 

IP2*IPT(IT0T3-t) 

X2»X0(IF2) 


GO  TO  23 

lose  (l)xl 

I0SC(2)=1 

I0SC(3)=1 

GO  TO  2k 

I0SC(4)xl 

I0SC(5)=1 

I0SC(6)xl 

60  TO  25 

IDSC(7)=1 

I0SC(8)xl 

I0SC(9)=1 
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r2»YOClP2l 

IF(SIOE(X1,V1,X2.Y2«XO,YOI .LT.C.O)  GO  TO  60 

IP3*IPT  (IT0T3) 

X3«X0IIP3) 

Y3»YOIIP3» 

IFISIOE(X2.Y2,X3.Y3«XO,YO) .LT.0.0)  GO  TO  60 

tFISIOE(X3«Y3«Xl«Yi,XO,YQI.LT.O.O)  GO  TO  60 

GO  TO  SO 

C CHECKS  IF  ON  THE  SANE  BOROEP  LINE  SEGHENT. 
fcO  IL1«IT0/NTL 
IL2«IT0-IL1*NTL 
IL1T3»IL1*3 
tPl»IPL<ILlT3-2> 

X1«X0IIP1I 

Y1»Y0(IP1» 

IP2=IPHIL1T3-1> 

X2*X0(IP2» 

Y2xrO(IP2l 

OX02xX0-X2 

OY02»Y0-Y2 

0X21»X2-X1 

0Y21*Y2-Y1 

CS0221*CX02*PX21*OY02*OY21 
TF(IL2.NE.ILli  GO  TO  50 

IF(CS0221.GT.0.0)  GO  TO  60 

OX01*X0-Xl 
OYOlsYO-Yl 

if(oyoi*ox21-oxoi*oy?i,gt.o.o)  go  to  50 

TF(OX01*OX21*OY01*OY21.LT.O.O)  GO  TO  60 

GO  TO  8P 

C CHECKS  IF  BETNF«-N  THE  SANE  TWO  50RCER  LINE  SEGMENTS. 

50  rF(CS0221.LT.0.0»  GO  TO  60 

TP3*IPL  C3*IL2-H 
X3»XOIIP3> 

Y3=YOCIP3» 

OX32*X3-X2 

0Y32«Y3-V2 

tFJOX02*OX32^OY02*OY3?.LE.n.C)  GO  TO  90 

C LOCATES  INSIOE  THE  DATA  AREA. 

C - DETERMINES  THE  SECTION  IN  WHICH  THE  POINT  IN  QUESTION  LIES. 
TSC»1 

tF(XO.GE.XSl)  ISC=ISC*1 

IF(X0.GE.XS2)  ISC-TSCTl 

TF(YO.GE.YSl)  ISC=ISC^3 

IF(Y0.GE. YS2)  ISC=ISCT3 

C - SEARCHES  THROUGH  THE  TRIANGLES  ASSOCIATED  WITH  THE  SECTION. 
NTSCI*NTSC(ISC> 

TF(NTSCI.LE.0»  go  TO  70 

JIMK*-9>ISC 
00  61  ITSC=1,NTSCI 
JINK>JIWK»9 
ITOxINKIJIHK) 

JWK«IT0*4 

IF»X0.LT.WKtJWK-3) ) GO  TO  61 

IF(X0.GT.HK(JWK-2) ) GO  TO  61 

IF{ Y0.LT.HK( JHK'l) > GO  TO  61 

IF(Y0.GT.MK(JMK»)  GO  TO  61 

IT0T3*ITn*3 
IPlxIFTaT0T3-2) 

XlxXDdPlI 
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yi«YO(iPi» 

1P2»IPTIIT0T3-1) 

X2«X0aP2) 

V2«YD(IP2) 

lF(S10EtXl.Yl,X2tY2«XO.YO) .LT.0.0)  GO  TO  61 

IP3«lPmT0T3) 

X3sXD(IP3) 

Y3»Y0aP3) 

IF(SIOE(X2«Y2«X3*Y3,XO.YO) .LT.3.01  GO  TO  61 

IF(SICE(X3fY3«Xl,Yl,XC,Y0l .LT.O.O)  GO  TO  61 

GO  TO  30 

61  CONTINUE 

C LOCATES  OUTSIDE  THE  DATA  ARfA. 

70  NL0T3aNL0*3 

IPl*IPL(NL0T3-2» 

XlsXO(IPl) 

Y1«Y0(IP1» 

TP2  = IPL  INL0T3-1I 
X2=XO«IP2> 

Y2=YOaP2» 

OX02«X0-X2 

OY02aY0-Y2 

OXElsXE-Xl 

0Y21=Y2-Y1 

CS0221«OX02*DX21^OY02*OV21 
00  72  IL2sltNL0 
XlrX2 
Y1»Y2 
DX01SOX02 
OY01*OY02 
CSPV=CS0221 
1P2*IPL(3*IL2-1) 

X2=XD«IP2I 

Y2=Y0(IP2» 

nX02»X0-X2 

OY02=Y0-Y2 

0X21*X2-X1 

0Y21=Y2-Y1 

CS0221=PX0?*CX21*OY02*OY21 
IFJCS0221.GT.0.0)  GO  TO  72 

IF(OX01*OX21*OY01*OY21.LT.O.C»  GO  TO  71 

IF(OY01*nX21-OX01*OY2l.LF.O.C»  GO  TO  74 

GO  TO  72 

•^1  IF<CSFV.GT.O.O)  GO  TO  73 

72  CONTINUE 
IL2  = 1 

73  tLlsIL2-l 
IF(lLl.EQ.C)  ILlrNLO 
GO  TO  75 

74  IL1*IL2 

75  IT0»IL1*NTL*IL? 

C NORMAL  EXIT 

30  ITI=IT0 
rTIPV«IT3 
RETURN 
'•NO 
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SUBROUTINE  lOPORV (NOP.XOtVOtZD.NCPtIPC .PO) 

c this  subroutine  estinates  partial  derivatives  op  the  first  ano 

C SECOND  order  AT  THE  DATA  POINTS. 

C THE  INPUT  PARANETEPS  APE 
C NOP  « NUMBER  OF  DATA  POINTS. 

C XO.YO.ZD  s ARRAYS  CONTAINING  THE  X,  Y,  ANO  Z COORQI- 

C NATES  OF  OATA  POINTS. 

C NCP  s NUMBER  OF  DATA  POINTS  TO  BE  USEO  FOR  ESTIMATION 

C OF  PARTIAL  DERIVATIVES  AT  EACH  OATA  POINT. 

C IPN  = ARRAY  CONTAINING  THE  POINT  NUMBERS  OF  NCP  DATA 

C POINTS  CLOSEST  TO  EACH  OF  THE  NCP  DATA  POINTS. 

C THE  OUTPUT  PARAMETER  IS 

C PD  = ARRAY  OF  DIMENSION  5*N0P.  WHERE  THE  ESTIMATED 

C ZX.  ZY.  ZXX.  ZXY.  ANO  ZYV  VALUES  AT  THE  DATA 

C POINTS  ARE  TO  BE  STORED. 

C DECLARATION  STATEMENTS 

DIMENSION  XD(IOO)  .YDdOOl  .ZO(100).IPC(400).PO<S00> 

REAL  NMX.NHY.NMZ.NHXX.NMXY.NMYX.NHYY 

C PRELIMINARY  PROCESSING 
10  NOPOxNOF 
NCPOsNCF 
NCPM1»NCP0-1 

C ESTIMATION  OF  ZX  ANO  ZY 
20  00  2h  IPQxl.NDPQ 

XOXXOIIPO) 

YO»YO(IPO» 

ZOxZDtIPO) 

NMXxO.O 

NMYxO.O 

NMZ«0.0 

JIPC0»NCP0*«IP0-1I 
DO  23  IClsl.NCPMl 
JIPCxJIPCO*ICl 
IPI*IPC<JIPC> 

OX1*XO(IPI)-XO 
OY1*YD«IPI»-YO 
OZlxZD«IPI>-ZO 
IC2MN«IC1*1 
DO  22  IC»=IC2MN.NCP0 
JIPC*JIPC0*IC2 
IFI=IPC(JIPC) 

OXZsXOCIPIl-XO 

OY2=YO(IPI)-YO 

OZ2*ZOCIPI)-ZO 

DNNX=DYl*CZ2-0Zl*rY2 

0NMY*DZ1*0X2-0X1*DZ2 

DNMZ*r)Xl*0Y2-DYl*CX2 

IF(ONMZ.GF.O.O)  GO  TO  21 

ONMXx-DNHX 

ONMYs-ONMY 

ONMZ*-ONMZ 

21  NMX*NMX»DNMX 
NMY*NMY»rNMY 
NMZxNHZtONHZ 

22  CONTINUE 

23  CONTINUE 
JPOO=E*IPO 
PO<JPOO-«»»  = -NMX/NMZ 
P0< JPOO-3)=-NMY/NMZ 

24  CONTINUE 


75 


C ESTINATION  OF  ZXX,  ZXY,  4ND  ZYY 

30  00  34  IPOsl,NOPO 

JPOO*JPPO*5 

XO>XOIIPO) 

JPOO*5*1PO 

YO*YOIIPO) 

ZXO«PC(JPOO-«.J 

ZYO*PC<JPOO-3) 

NMXX«0. 0 
NMXYsO.O 
NHYXsQ.O 
NMYY=0.0 
NMZ  sO.O 

JIPC0=NCP0*(IPC-1I 
00  33  IC1=1,NCPH1 
JIPCsJIPCO^ICl 
IP1=IPC(JIPC» 

DX1=XD(1PI)-X0 

OY1=YO(IPI)-YO 

JPD=5*IPI 

ozxi=PO(jPO-4)-zxn 
0ZYl=PD«JPD-3)-ZYD 
ICZHN=IC1+1 
DO  32  1C2=IC2MN«NCP0 
jrPC=JIPC0*IC2 
IFI=IPCIJIPC) 

OX2=XO(IP1»-XO 

DY2=Y0(IPI)-YD 

JP0»5*IPI 

OZX2sPO(JPD-i»)-ZXO 

CZY2sPn<JPO-3>-7YO 

DNMXX=DY1*0ZX2-0ZXI*0Y2 

DNIXYsDZXl'OXZ-OXl’OZXa 

ONMYXsDYl*OZY2-OZYl*OY2 

0NMYYs['ZYl*DX2-0Xl*0ZY2 

DNMZ  sOXl*OY2  -"Y1*0X2 

IFCDNMZ.GE.O.OI  GO  TO  31 

CNMXX=-DNMXX 

ONMXY=-ONMXY 

DNMYX»-ONMYX 

ONMYYs-ONMYY 

ONHZ  =-ONHZ 

31  NHXXsNMXX*ONMXX 
NPXY=NMXY»ONMXY 
NHYX=NMYX*ONMYX 
nhyy=nmyy*onmyy 
NMZ  »NMZ  ♦DNMZ 

32  CONTINUE 

) 33  CONTINUF 

i PO  IJPDfJ-2)  = -NMXX/NMZ 

PO(JPCO-1)=-<NMXY^NMYX)/I2.0*NMZ) 
POCJPOO)  =-NMYY/NMZ 


^4  CONTINUE 
i^ETURN 
ENO 


subroutine  IDPTIPlXO,VOfZO,NT,IPT,NLtIPL*POOtITI,XII,YII, 
1 ZII) 

C THIS  SUBROUTINE  PERFORMS  PUNCTUAL  INTERPOLATION  OR  EXTRAPO- 
C LATION,  I.E..  DETERMINES  THE  Z VALUE  AT  A POINT. 

C THE  INPUT  PARAMETERS  ARE 

C XO.YD.ZO  » ARRAYS  CONTAINING  THE  X,  Y,  AND  Z 

C COORDINATES  OF  DATA  POINTS, 

C NT  = NUMBFR  OF  TRIANGLFS, 

C TPT  * ARRAY  CONTAINING  THF  POINT  NUMBERS  OF  THE 

C VERTEXES  OF  THE  TRIANGLES, 

C NL  « NUMBER  OF  BORDER  LINE  SEGMENTS, 

C IPL  « ARRAY  CONTAINING  THE  PCINT  NUMBERS  OF  THF  END 

C POINTS  OF  THE  BORCEP  LINE  SEGMENTS  AND  THEIR 

C RESPECTIVE  TRIANGLE  NUMBERS, 

C POO  = ARRAY  CONTAINING  THE  PARTIAL  DERIVATIVES  AT 

c the  data  points, 

C ITI  = TRIANGLE  NUMBER  OF  Thp  TRIANGLE  IN  HHICH  LIES 

C THE  POINT  FOR  HHICH  INTERPOLATION  IS  TO  BE 

C PERFORMED, 

C XII, YII  = X AND  Y COORDINATES  OF  THE  POINT  FOR  HHICH 

C INTERPOLATION  IS  TO  BE  PERFORMED. 

C THE  OUTPUT  PARAMETER  IS 
C ZII  = INTERPOLATED  Z VALUE. 

C DECLARATION  STATEMENTS 

DIMENSION  XO(IOO) ,YO(100) ,ZO(100) , IPT (585 » , IPL ( 30 0 ) , 

1 P00(500) 

COMMON/IOPI/NCF 

DIMENSION  X(3I ,Y(3> ,Z(3),PO(15> , 

1 ZU(3) ,ZV(3)  ,ZUU(3) ,ZUV(3) ,ZVV(3) 

REAL  LU,LV 

EQUIVALENCE  (P5,P05) 

DATA  NCFMX/EO/ 

C SETTING  OF  SOME  LOCAL  VARIABLES. 

10  IT0*ITI 
XI0=XII 
YIOaYII 
NTL=NT*NL 

C DETERMINES  IF  SIMPLE  INTERPOLATION  IS  APPLICABLE. 

?0  IFdTO.LE.NTL)  GO  TO  30 

IL1=IT0/NTL 
ILZ=IT0-IL1*NTL 
IL1T3=IL1*3 
tL2T3=IL2*3 
IT0  = IPL  (IL1T3> 

IF(IL1.NE.IL2)  GO  TO  4Q 

C CALCULATION  OF  ZII  BY  SIMPLE  INTERPOLATION  OR  EXTRAPOLATION. 

70  ASSIGN  31  TO  LPL 
GO  TO  50 
31  ZIIsZIO 
RETURN 

C CALCULATION  OF  ZII  AS  A HEIGHTEO  MEAN  CF  THO  EXTRAPOLATED 
C VALUES. 

40  ASSIGN  41  TO  LBL 
GO  TO  50 

41  ZI1=Z10 
ITOsIPL  (IL2T,7) 

ASSIGN  42  TO  LBL 
GO  TO  50 

42  ZT2=ZI0 

C CALCULATES  THE  HEIGHTING  COEFFICIENTS  FOR  EXTRAPOLATED  VALUES. 
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1.5  IPl»IPL«lLlT3-2> 

IP2»IPL(IL173-1) 

IP3=1PL(IL2T3-1) 

XlxXOIlPl) 

YlsYOClPlI 

X2*X0(IP2) 

Y2*yO(lP2) 

X3«X0«1P3» 

Y3aYOtIP3» 

OX02*XI0-X2 

OYOZaYIO-YS 

0X32aX3-X2 

0Y32SY3-Y2 

0X21=X2-X1 

OY21*Y2-Yl 

Ml*  <OX02*OX32+OYO2*DY32»  **2/ ( 0X32* DX 32 ♦0Y32*DY32 ) 

M2* (DX02*OX21*OY02*OY21I**2/ I0X21*0X21*0Y21*0Y21) 

C CALCULATES  ZII  AS  A METGHTEP  MEAN. 

46  ZII=<M1*ZI1*M2*ZI2) /|M1*M2» 

9ETURN 

C INTERNAL  ROUTINE  FOR  PUNCTUAL  INTERPOLATION. 

C CHECKS  IF  THE  NECESSARY  COEFFICIENTS  HAVE  BEEN  CALCULATED. 

50  IF(NCF,EO,0)  60  to  60 

IFCNCF.EQ.ITO)  go  TO  TO 

C CALCULATION  OF  MEM  COEFFICIENT  VALUES. 

C determines  the  COEFFICIENTS  FOR  THE  COORDINATE  SYSTEM  TRANS- 
C FORMATION  FROM  THE  X-Y  SYSTEM  TO  THE  U-V  SYSTEM,  AND  CALCU- 
C LATES  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  FOR  INTERPOLATION. 

C LOADS  COORDINATE  AND  PARTIAL  DERIVATIVE  VALUES  AT  THE 
C VERTEXES. 

60  NCF*ITO 
JIPT=3*(IT0-1) 

JPO*0 

00  62  1=1,3 

JIPT=JIPT*1 
IOP»IPT(JIPT) 

X(I)*XO  (lOP) 

Y(I>*Y0<10P) 

Z(I)=ZO(IOP) 

JPOO*5*<IDP-1) 

00  61  KP0=1,5 

JPD*JP0*1 
JPOO=JPOD*1 
PDJJPDlaPODC JPCDI 

61  CONTINUE 

62  CONTINUE 

C OFTFRMINING  THE  COEFFICIENTS  FOR  THE  COORDINATE  SYSTEM 
C TRANSFORMATION  FROM  THE  X-Y  SYSTEM  TO  THE  U-V  SYSTEM 
C AND  VICE  VERSA 

63  X0*X(1) 

Y3»Y  ru 
AaXIZl-X" 

R=X<3)-X0 
C*Y (2J-Y0 
0=V»3)-YC 
A0»A*0 
RC»B*C 
1LT«A0-BC 
!•»  0/CLT 
*»«-f/OLT 
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CP*-C/OLT 
iJP*  A/OLT 

C C0NV5PSI0N  OF  THE  PARTIAL  DERIVATIVES  AT  THE  VERTEXES  OF  THE 

C TRIANGLE  FOR  THE  U-V  COORDINATE  SYSTEM 

FA  AA3A*A 

ACT2*2.0*A*C 

CC=C*C 

A9sA*B 

AOBC«AD«BC 

CD=C*0 

BB=B*B 

BOT2=2.0*B*O 

0D*O*O 

DO  65  1=1,3 

JP0*5*I 

ZU(  II=A*PO( JPO-4) ♦C*PP(JPO-3) 

ZV«I)=9*P0(JPD-4» ♦D*PO(JPD-3) 

ZUUII)=AA*P0IJP0-2) fACTZ^PD ( JPD-1) ♦CC’PD (JPD) 
ZUVCI)=A3*P0< JPO-2) ♦AOBC*Pr(JPD-l) ♦CO*PO<JPO> 
ZVVa>=BB*PO(  JP0-2)*B0T2*P0(JPD-1)^0D*P0(JP0» 

65  CONTINUE 

C CALCULATION  OF  THE  COEFFICIENTS  OF  THE  POLYNOMIAL 

66  P00=Z(1> 

pio=zua» 

poi=zva» 

P20=0.5*ZUU(1) 

P11=ZUVI1) 

P02=0.5*ZVV(1) 

H1=Z<2)-P00-P10-P20 

H2»ZUt2)-P10-ZUU«l) 

H3=ZUU(2)-ZUU(1) 

P30=  10.0*H1-4.0*H2^1.5*H3 
P40=-15.0*Hl*7.0*H2  -H3 

050=  6.0*H1-3.0*H2»0.5*H3 

H1=Z(3)-PCO-P01-P02 
H2=ZV(3»-P0l-ZVV(U 
H3=ZVV<3)-ZVV(1I 
“03=  10.0*H1-4,0*H2*0.5*H3 
P04=-15.0*Hlf7.0*H2  -H3 

“05=  6.0*H1-3.0*H2*0.5»H3 

LU=SQRT  (AAFCC) 

LV=SQRT  IBB*OOI 
THXU=ATAN2(C,A) 

THUV»ATAN?tD,B)-THXU 

CSUV=COS(THUV> 

“41=5.0*LV*CSUV/LU*PEC 

“14  = 5.0*LU»CEUV/LV*P0'i 

H1=ZV(2)-P01-P11-P41 

H?*ZUV»2)-Pll-4. 0*P41 

“21=  3.n*Hl-H2 

“31=-2.0*Hl4H2 

H1=7U(3»-P10-P11-P14 

H2=ZUVJ3)-Pll-4.0*P14 

“12=  3.0*M1-H2 

“13=-2.0*H1*H2 

THUS«ATAN?(0-C,B-A) -THXU 

THSV»THUV-THUS 

AA=  SIN<THSVI7LU 

3B»-C0S<THSV)/LU 

CC«  SIN(THUS»/LV 
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00a  COS (THUS) /LV 
«CaA««CC 
A0aAA*00 
eCaBB*CC 

GlaAA^AC* t3.0*BC+2.0*AD) 

G2aCC*AC*(3.0*A0^2.0*BC) 

Hla-AA*AA*AA*(5.0*AA*8B*P50* 14. O*0C*AO) •P41) 

1 -CC*CC*CC* (5.0*CC*00*P05* (4.0*AO*BC) •P14) 
H2a0.5*ZVVI2)-P02-P12 
H3a0.5*2UU(3)-P20-P21 
P22=(G1*H2*G2*H3-H1) /(G1^G2) 

P32=H2-P22 

P23aH3-P22 

C TRAN3F0RMATICN  OF  THE  COORDINATE  SYSTEM  FROM  X-V  TO  U-V 
70  OXaXII-XO 
OY=YII-YO 
UaAP*CX+BP*Cr 
V=CP*OX*OP*DY 

C EVALUATION  OF  THE  POLYNOMIAL 

75  PO  = POO  + U*  (P10»U* (P2  3*U* (P30*U* (P40*U*P50) ) ) ) 
P1=P014U*(P11*U*(P21»U*(P31»U*P41) ) ) 
P2aP02*U*(P12^U* (P22+U*P32) ) 

P3aP03*U*(P13*U*P23) 

P4=P04>U*P14 

ZIO  = PO  + V*(P1*V*(P2*V*(P3*V*(P«»*V*PE)  ) ) ) 

GO  TO  LBL.  (31,41,42) 

END 


SUBKOUTIMr  IOT*NG(NOP,XO«YOtNT,IPT«NLtIPL«IHL,IHP,HK) 

C THIS  SUBROUTINt  PERFORMS  TRI RNGUL AT  ION , IT  DIVIDES  THE  X-Y 
C PLAN?  INTO  A NUMBER  OF  TRIANGLES  ACCORDING  TO  GIVEN  DATA 
C POINTS  IN  TME  PLANE,  DETERMINES  LINE  SEGMENTS  THAT  FORM  THE 
C BORDER  OF  DATA  AREA,  AND  DETERMINES  THE  TRIANGLE  NUMBERS 
C CORRESPONDING  TO  THE  BORDER  LINE  SEGMENTS. 

C AT  COMPLETION,  POINT  NUMBERS  OF  THE  VERTEXES  OF  EACH  TRIANGLE 
C ARE  LISTED  COUNTER-CLOCKNISF.  POINT  NUMBERS  OF  THE  END  POINTS 
C OF  EACH  BORDER  LINE  SEGMENT  APE  LISTED  COUNTER-CLOCKMISE, 

C LISTING  ORDER  OF  THE  LINE  SEGMENTS  BEING  COUNTE R-CLOCKN ISE , 

C THE  INPUT  PARAMETERS  ARE 
C NDP  = NUMBER  OF  DATA  POINTS, 

" C XO  = ARRAY  OF  DIMENSION  NOP  CONTAINING  THE 

C X COORDINATES  OF  THE  DATA  POINTS, 

C YD  = ARRAY  OF  DIMENSION  NOP  CONTAINING  THE 

C Y COORDINATES  OF  THE  DATA  POINTS. 

C THE  OUTPUT  PARAMETERS  AR'^ 

C NT  = NUMBER  OE  TRIANGLES, 

C IPT  r ARRAY  OF  DIMENSION  F*NCP-15,  WHERE  THE  POINT 

C NUMBERS  OF  THE  VERTEXES  OF  THE  tIT)TH  TRIANGLE 

C ARE  TO  BE  STDREO  AS  THE  (3*IT-2)ND,  (3*IT-1)ST, 

C AND  I3*IT)Th  el<^mFNTS,  IT  = 1 ,2,  . . . ,NT  , 

C NL  = NUMBER  OF  BORDER  LINE  SEGMENTS, 

C IPL  = AFRAY  OF  DIMENSION  6*N0P,  WHERE  THE  POINT 

C NUMBERS  OF  THE  ENO  POINTS  OF  THE  (IL)TH  BORDER 

C LINE  segment  and  ITS  RESPECTIVE  TRIANGLE  NUMBER 

C APE  TO  BF  STORTO  AS  THE  (3*IL-2»N0,  (3*IL-1)ST, 

C AND  (3*IL)TH  ElEMFNTS,  IL=1,2,...,  NL. 

C THE  OTHER  PARAMETERS  ARE 

C TWL  = INTEGER  ARRAY  OF  DIMENSION  18*N0P  USED 

C internally  as  a work  AREA, 

C IWP  = INTEGER  ARRAY  OF  DIMENSION  NCF  USED 

C INTERNALLY  AS  A WORK  AREA, 

C WK  = ARRAY  OE  DIMENSION  NOP  USED  INTERNALLY  AS  A 

C WCRK  AREA. 

C DECLARATION  STATEMENTS 

OIMFNSION  XO(IOO) ,YO(19P> ,IPT(5*S) ,IPL(600) , 

1 THL(ie03)  ,THPJ13C>,WK(100) 

OTMFNSICN  rTE(2) 

DATA  RATIO/1. OE-6/ , NRFp/100/,  LUN/6/ 

c statement  functions 

0SQF(X1,Y1,X2,Y2)  = (X’-X1)**2+  (Y2-Y1) **2 

3IDEIX1,Y1,X2,Y2,X3,Y3)= (Y3-Y1)*(X2-X1) -(X3-X1)*(Y2-Y1) 

C PRELIMINARY  PROCESSING 
10  NOPO*NCF 

NOPM1=NOPO-1 

IF(N0PD.LT.<.)  GO  TO  B" 

C O^TERMINES  THE  CLOSFST  dAIP  OF  DATA  POINTS  AND  TH':iR  MIDPOINT. 
PO  DSOHNsOSOFIXOd)  , YO  ( I ) , XO  ( ? ) , YE  ( 2 ) ) 

IFMNlal 

IPMN2=2 

1 DO  2?  IPl=l,NrPMl 

) X1=XD(IP1) 

Y1=YD  tIPl) 

’ IP1P1=IP1*1 

00  21  IP?= IPIPI , NO09 

OSCI  = OSC'F(Xi,yi  ,yD(IP2)> 

IFIOSOI.EQ.O)  GO  TO  B1 

IF  (CSQI.GE.OSQMN)  GO  TO  21 

OSOMNrOSOI 
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1PHN1«IP1 

IPMN2=1P2 

21  CONTINUE 

22  CONTINUE 
0SQ12s0SQMN 

XOMP«(XD(IPMNl)»XO  (TPNN2n/2.0 
YOMPs IY0(IPNNt)«YO(IPNN2))/2.0 
C SORTS  THE  OTHER  (NDP-21  OXTX  POINTS  IN  ASCENDING  ORDER  OF 
C DISTANCE  FRCH  THE  MIDPOINT  AND  STORES  THE  SORTED  DATA  POINT 
C NUMBERS  IN  THE  IMP  ARRAY. 

30  JPls2 

DO  31  IP1=1,NDP0 

IF(IP1,EQ.IPMN1.0R.I«»1.EQ.IPHN2)  GO  TO  31 

J?i«JPl>l 

IMP(JP1)=IP1 

MK( JPl):OSaF(XOHP«YOHP«XO(IPl) .YDIIPI) ) 

31  CONTINUE 

00  33  JP1=3,N0PM1 

0SQHN*MKCJP1) 

JPHNrJPl 

DO  32  JP2=JP1,N0P1 

1F«NKIJP2).GE.0S0MN)  GO  TO  32 

DSQHN=HK(JP2) 

JPMNrJP2 

32  CONTINUE 
ITSsINPlJPl) 

IMPJJP1)=IMP< JPMN) 

IMP  «JPMN)=ITS 
MK(JPMNI=MK(JP1I 

33  CONTINUE 

C IF  NECESSARY,  MODIFIES  THE  OPf'ERING  IN  SUCH  A MAY  THAT  THE 
C FIRST  THREE  DATA  POINTS  APE  NOT  CCLLINEAR. 

3E  AR»0SQ12*RATI0 
xisxdiifmnh 

YlrYOdPMNlI 
0X2lrX0dPMN2)-Xl 
0Y21=Y0CIPMN2)-Y1 
00  36  JP=3,NOPO 
IP*IMF(JP) 

IF( ABE ( (Y0(IF>-Y1 )*0X21-(X0(IF»-X1)*0Y21) .GT.AR) 

1 GO  TO  17 

36  CONTINUE 
GO  TO  92 

37  IF(JP.EC.3)  GO  TO 
JPMX« JP 

IMP  «3)  = IP 
JP= JPMX+1 
00  38  J0C*4,JPMX 
JP=JP-1 

IMP(JP)rlMP(jP-lJ 

18  CONTINUE 

r EOPHS  THE  FIRST  TRIANGLE.  STORES  POINT  NUMBERS  OF  THE  VE=»- 
C TEXES  OF  THE  TRIANGLE  IN  THE  IRT  ARRAY,  AND  STORES  POINT  NUM 
C BERS  OF  THE  BORDER  LINE  SEGMENTS  AND  THE  TRIANGLE  NUMBER  IN 
C THE  IPL  ARRAY. 

<*f)  IP1  = IPHM 
TP2=IPMN2 
IP1=IMF (3) 

IF(SIDE(xn»IPl)  ,YOfIPl) ,X0»IP2) ,Y0«IP2) ,X0(IP3) ,Y0(IP3n 
1 .GE.0.0)  GO  ’0  41 


IP1*IPMN2 
IP2«IPtiNl 
41  NT0=1 
NTT3»3 
IPT<1)*IP1 
IPT (2)=IP2 
IPT<3»=IP3 
NL0  = 3 
NLT3*9 
IPL I1)=IP1 
IPL(2»=IP2 
IPL(3)=1 
TPL (4)=IP2 
IPL  »5I  = IP3 
IPL(6lsl 
IPL<7)sIP3 
IPL<9»*IP1 
IPL  f9)  = l 

C ADOS  THE  REMAINING  (NDP-3)  TATA  POINTS,  ONE  BY  ONE. 

=0  00  79  JP1=4,NOPO 
IP1  = INP(JPU 
X1=X0IIP1) 

YlsYOCIPl) 

C - DETERMINES  THE  VISIBLE  BORDER  LINE  SEGMENTS. 

IP2»IPL<1) 

JPNN*1 

0XMN=X0(IP2)-X1 
0YHNsY0JIP2)-Yl 
ARMN= (0XMN**2*0YMN**2»*RATI0 
JPMX=1 
OXHX=CXHN 
OYMXsCYMN 
ARMXsARMN 
00  52  JP2=2,NL3 

IP2=IPL(3*JP2-2) 

OXsXO  (IP2)-X1 
0Y  = Y0aP2)-Yl 

IF(CY*OXMN-OX*OYMN.GE.-AFMN)  GO  TO  51 

JPMN=JP2 

OXMM=OX 

DYMNrOY 

ARMN=  «0XMN**2*0YMM**2) 'RATIO 
GO  TO  52 

51  if<oy*dxmx-dx'oyhx.lt.-armx)  go  to  52 

JPMX=JP2 

OXMXsOX 

OYMXsOY 

ARMX=(0XHX**2'CYBX'*2> 'RATIO 
■^2  CONTINUE 

IFCJPMX.LT.  JPHN)  JPMX=JP‘1X'NL0 
NSH=JFMN-1 

IF(NSH.LE.O)  GO  TO  60 

r - SHIFTS  (ROTATES)  THE  IPL  ARRAY  TO  HAVE  THE  INVISIBLE  BORDER 

C - LINE  SEGMENTS  CONTAINEO  IN  THE  FIRST  PART  OF  THE  IPL  ARRAY. 
NSHT3=NSH'3 

00  53  JP2T3=3,NSHT3,3 

JP3T3=JP2T3*NLT3 
IPL  tJP3T3-2)=IPL(JP2T3-2) 

IPL (JP3T3-1) =IPL( JP2T3-1) 
iPL(JP3T3)  =IPL(JP2T3) 
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53  CONTINUE 

00  5t  JP2T3«3tNLT3»3 
JP3T3=JP2T3*NSHT3 
IPL (JP2T3-2»*IPL ( JP3T3-2) 

IPL (JP2T3-1ISIPL (JP3T3-1) 

IPL(JP2T3»  =IPHJP3T3» 

54  CONTINUE 

JPMXxJPMX-NSH 

C - ADOS  TRIANGLES  TO  THE  IPT  ARRAY,  UPDATES  BORDER  LINE 
C - SEGMENTS  IN  THE  IPL  ARRAY,  AND  SETS  FLAGS  FOR  THE  BORDER 
C - LINE  SEGMENTS  TO  BE  REEXAMINED  IN  THE  INL  ARRAY, 

60  JNL-O 

DO  64  JP2=JPMX,NL') 

JP2T3*JP2*3 
IPL1=IPL( JP2T3-2I 
IPL2=IPL(JP2T3-1) 

IT  =IPLIJF2T3) 

C ADDS  A TRIANGLE  TO  THE  IPT  ARRAY. 

NT0=NT0*1 
NTT3=NTT3*3 
IPT (NTT3-2)=IPL2 
IPT«NTT3-1)=IFL1 
IPT<NTT3)  =IP1 

C UPDATES  BORDER  LINE  SEGMENTS  IN  THE  IPL  ARRAY. 

IFIJP2.NE. JPMX)  GO  TO  61 

IPL«JP2T3-1I=IP1 
IPL(JP2T3)  =NTO 

61  IF(JP2.NE.NL0)  GO  TO  62 

NLN=JPMX*1 
NLNT3=NLN*3 
IPLJNLNT3-2»=IP1 
IPL  (NLNT3-1I »IPL(1» 

IPL(NLNT3)  =NT3 

- - DETERMINES  THE  VERTEX  THAT  DOES  NOT  LIE  ON  THE  BORDER 
LIN6  SEGMENTS. 

62  ITT3=IT*3 
IPTI=TPT(ITT3-2» 

IFIIPTI.NE.IPL1.AND.1PTI.NE.IPL2)  GO  TO  63 
IPTI=IPT<ITT3-1) 

IF(IPTI.NE.IPL1.AN0.1FTI.NE.IPL2)  GO  TO  63 
IPTIsIPT(ITT3l 
C - - CHECKS  IF  THE  FXCHANGr  TS  NECESSARY. 

-3  IF(I0XCHG(XD,Y0,IP1,IPTI,IPL1,IPL2) .EO.O)  GO  TO  64 

C - - MODIFIES  THE  IPT  ARRAY  WHEN  NECESSARY. 

IPT<ITT3-2)=IPTI 
IPTIITT3-1)=IPL1 
IPHITT3)  =IP1 


IPT<NTT3-1)=IPTI 

IF«JP2.E0.JPMX»  IPL  (JP2T3)*IT 

IF<  JP2.EO.NLO.  AN0.IPH3)  .EO.ITI  IPL  t3)  sNTO 

c - - sets  flags  in  the  iml  array. 

JML=JML«4 
1ML(JNL-3)*IPL1 
INL  (JNL-2»  = IPTI 


I 
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IF(NLF.EQ.O»  60  TO  79 

C - IHPROWES  TRIANGULATION. 

^0  NTT3P3=NTT3t3 

DO  78  IREPsl.NREP 
DO  76  ILF>1,NLF 
ILFT2=ILF*2 
IPL1=IML(ILFT2-1) 

IPL2=INU(ILFT2» 

C LOCATES  IN  THE  IPT  ARRAY  TNO  TRIANGLES  ON  BOTH  SIDES  OF 

C THE  FLAGGED  LINE  SFGHENT, 

NTF=0 

DO  71  ITT3R=3,NTT3,3 
ITT3=NTT3P3-ITT3R 
IPTl  = IPTa7T3-2> 

IPT2=IPT(ITT3-1) 

IPT3=IFHITT3) 

IFdPLl.NE.  IPTl.  AN0.IPL1.NE.IPT2.AN0. 

1 IPL1.NE.IPT3»  GO  TO  71 

IF(IPL2.NE.I'>T1.5ND.IPL2.NE,IPT2,AND. 

1 IPL2.NE. IPT3)  GO  TO  71 

NTF=NTF«-1 
ITF(NTF)=ITT3/3 
IF(NTF.EQ.2)  GO  TO  72 

71  CONTINUE 

IF(NTF,LT.2)  GO  TO  76 

C OETERHINES  THE  VERTEXES  OF  THE  TRIANGLES  THAT  DO  NOT  LIE 

C - - ON  THE  LINE  SEGNENT. 

72  IT1T3=ITF(1)*3 
IPTI1=iptiIT1T3-2) 

IFTIPTIl.NC.IPLl.ANO.rFTIl.NE.IPL?)  CO  TO  73 

IFTI1=IPT (IT1T3-1) 

IFdPTIl.NE.IPLl.  AN0.IPTI1.NE.IPL2)  GO  TO  73 

IPT11=IPT«IT1T3) 

^3  IT2T3=ITF(2)*3 

IPTI2=IPT (IT2T3-2) 

IFCIPTI2.NE.IPL1.AN0.IPTI2.NE.IPL2)  GO  TO  74 

IFTI2=IPTCIT2T3-1) 

IF(IPTI2.NE.IPLl.ANn.IPTI2.NE.lPL2)  GO  TO  74 

IPTI2=IPT  tIT2T3) 

C - - CHECKS  IF  THE  EXCHANGE  IS  NECESSARY. 

74  IFCIOXCHGCXO. Y0,IFTI1,IPTI2,IPL1.IPL2) .EO.O) 

1 GO  TO  76 

C - - NnOIFIES  THE  IPT  ARRAY  MH^N  NECESSARY. 

IPT (IT1T3-2) =IPTI1 
IPT(IT1T3-1)=IPTI2 
IPTCIT1T3)  =IPL1 
IPT(IT2T3-2)  = I'»TI2 
IPTCIT2T3-1)sIdTI1 
IPT(IT2T3)  =I^L2 
C - - SETS  NEM  FLAGS. 

UHLsUNL^B 

IICL  CJNL-7)  =IPL1 

INL(JNL-6)=IPTI1 

IML<JNL-5) sIPTIl 

IHL(JML-4)=IPL? 

IML  IJNL-3)  = IPL2 
INLIJML-2)=IPTI9 
INL(JNL-l) =IPTI2 
IHL(JML)  sIPLl 
OP  75  JLT3=3.MLT3,3 
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IPLJ1*1PL(  Jt.T3-?) 

IPLJ?=IPL(JLT3-1) 

IF( IIPLJI.EQ. IPL1.AND.IPLJ2.EQ.IPTI2) .OP. 

1 IIPLJ2.E0. IPLl.»Nr.lPLJl.E0.IPTI2)) 

2 IPHJLT3)  =ITF  (1) 

IF C IIPLJI.EQ.  I«>L2.*ND.IPLJ2.EQ.IPTI1)  .OR. 

1 (IPLJ2.EQ.IPL2.«NC.IPLJ1.EQ.IPTI1)) 

? IPLI JLT3I =ITF(2) 

’S  CONTINUE 

76  CONTINUE 
NLFC=NLF 
NLF=JNL/2 

IFINLF.EQ.NLFC)  GO  TO  7q 

C - - RESETS  the  IHL  ARRAY  FO*"  THE  NEXT  ROUND. 

JHL*^0 

J»IL1MN=INLFC*1»*2 

NLFT2»NLF*2 

DO  77  JML1=JML1MN,NLFT2,2 
JNL=JML*2 

I NL (JML-1 ) = I NL I JNL 1- 1 » 

IhLIJHL)  =INL(JML1) 

77  CONTINUE 
HLFsJNL/2 

78  CONTINUr 
70  CONTINUF 

C REARRANGE  THE  IPT  ARRAY  SO  THAT  THE  VERTEXES  OF  EACH  TRIANGLE 
C ARE  LISTED  COUNTER-CLOCKMISE . 

80  00  81  ITT3*3,NTT3,3 

IP1  = IPT  (ITT3-2) 


IP2=IPTIITT3-1> 

IP3SIFTIITT3I 

IF (SIDE  IXOIIPH .YOITPII ,X0(IP2) ,Y0(IP2J ,X0(IP3» ,YDIIP3)) 
1 .GE.0.0)  GO  TO  81 

IPT(ITT3-2)=IP2 
IPTCITT7-1) =IP1 


81  CONTINUE 
NTsNTO 
NLsNLO 
RETURN 
C FRROR  EXIT 

90  WRITE  CLUN,2C9P)  NOPO 
GO  TO  93 

91  WRITE  (LUN,2091)  NOR  0 , IPl . IF  2 .XI , Y1 
GO  TO  93 

92  WRITE  (LUN, 20921  NQRO 

93  WRITE  ILUN,2‘393) 


RETURN 

C FCPWAT  STATEMENTS 

2090  F0RMAT(1X/23H  NO®  LESS  THAN  4./SH  NOP  *,I5) 

2091  F0RMATI1X/29H  •••  TOENTICAL  OATA  POINTS./ 

1 8H  NOP  =,IE,tX,^HIPl  *.I';,EX,9HI02  =.IE. 

2 5X,4HXD  *,E12.4,9X,4HY0  «,E1?.4> 

RT-'?  FORMAT  riX/33H  •••  8LL  COLLINEAR  DATA  POINTS./ 
f 8H  NOP  =.I5) 

20PT  FORMATI’SH  EpRO®  OETFCTEP  IN  ROUTINE  tOTANG/) 

PNO 


86 


1 


FUNCTION  lOXCHGtX, 

C THIS  FUNCTION  OfTERNINCS  MHETHER  OR  NOT  THE  EXCHANGE  OF  TMO 

c triangles  is  necessary  on  the  basis  of  hax-hin-angle  criterion 

C BY  C.  L.  LAMSON. 

C THE  input  PARAHETERS  ARE 

C XtY  « ARRAYS  CONTAINING  THE  COORDINATES  OF  THE  DATA 

C POINTS, 

C 11,12,13,14  = POINT  NUMBERS  OF  FOUR  POINTS  PI,  P2, 

C P3,  ANO  P4  THAT  FORH  A QUADRILATERAL 

C NITH  P3  AND  P4  CONNECTED  OIADONALLY. 

C THIS  FUNCTION  RETURNS  A VALUE  1 (ONE)  WHEN  AN  EXCHANGE  IS 
C NEEDED,  ANO  0 (ZERO)  OTHERNISE. 

C DECLARATION  STATEHENTS 

OIHENSICN  X(IOO) ,Y(100) 

DIMENSION  X0(4),Y0I4) 

EQUIVALENCE  (C2SQ,C1SQI , (A3SQ,B2SQ) , (B3SQ, AISQ) , 

1 tA4SQ,BlSQ) , (B4S0,A2SQ) , (C4SQ,C3SQ) 

C STATEMENT  FUNCTIONS 

SLSQ(J1  ,J2)s(X0(J2)  -X0(  J1H**2*(Y0(J2)-V0(  Jl))**2 
AREA( Jl,J2,J3)s(Y0 ( JI)-YO(Jl) )»(X0 (J2)-X0 ( Jl) ) 

1 -(XO ( J3)-X0 (Jl) )• (YO (J2)-Y0  (Jl) ) 

C CALCULATION 

10  XOdlsXdl) 

Y0(l)>Ydl) 

X0(2)«X(I2) 

Y0(2)=Yd2) 

X0(3)3Xd3) 

Y0(3)sY)I3) 

Y0(4)sX(I4) 

Y0(4)sYd4) 

20  IOX*0 

U3>AREA(3,1,2) 

U4sAREA(4,2,l) 

IF(U3*U4,LE.0.0)  GO  TO  3P 
UlsAREA(l,4,3) 

UZsAREA  (2,3,4) 
tlSQxSLSQ(3,lt 
31S0«SLS0(l,i») 

C1SQ*SLSQ(4,3) 

A?SQ«SLSQ(4,2) 

B?SQxSLSO(2,3) 

C3SQ*SLSO(l,2) 

S1S0*U1*U1/(C1SQ*AMAX1( A1SQ,P1S0) ) 

S2SQ*U2*U2/(C2S0*AMAX1( A2SQ,B2SO) ) 

S3SQ*U3*U3/(C3rQ*AMAXl (A3S0,e3S0) ) 

S4SQ*U4*U4/  (C4SQ*AMAXl(Ai.SQ,e«.S0)  ) 
tF(AHINl(SlSa,S2SO) , LT . AMIMl ( S3SQ , S4SQ) ) IOX=l 

30  lOXCHGxIOX 
RETURN 
END 


CALL  TABLES  TO  INTERPOLATE 
AVERAGE  YEARLY  PRESSURE,  TEMP. 
REL.  HUMIDITY  BETA  AND  M. 


READ 

ANY  PRESS.,  TEMP., 
\ R.H.  OR  M?  / 


USE 

KNOWN 

VALUES 


USE 

INTERPOLATED 

VALUES 


A 


I 

I 


OBTAIN 
VARIANCES 
OF  M AND 
BETA 


CALCULATE 

METHOD 

3 


CALCULATE 

METHOD 

4 


CALCULATE 
ATTENUATION 
FOR  METHODS 
1 AND  2 


CALCULATE  5 AND  95% 
CONFIDENCE  LIMITS  OF 
ATTENUATION  DISTRIBU- 
TION FOR  ALL  METHODS 
(1  - 4) 


OUTPUT 

RESULTS 
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